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A  STUDY  CP  THICK-TARGET  X-RAY  SPECTRA 


USING  PHOTONUCIEAR  REACTIONS 


By;  Karl  L,  Prado 


Megavoltage  bremsstrahlung  spectra  which  are  produced  in  thick 
targets  were  investigated  using  photonuclear  reactions.  Ratios  of 
photonuclear  yields  in  teflon  and  in  KjjsiFgTwere  obtained  that  serve  as 
eeneitive  indexes  of  thick  target  x-ray  quality  and  endpoint  energy  in 
the  energy  range  from  about  14  to  30  MeV.  Because  thick  target 
bremsstrahlung  epectra  have  never  before  been  used  in  the  determination 
of  photonuclear  reaction  croee  eect-ions,  a  numerical  analysis  computer 
calculation  of  thick  target  bremsstrahlung  spectra  was  written  to 
approximate  the  epectra  of  x-rays  produced  in  thick  targets  of  clinical 
electron  linear  accelerators  and  examine  the  feasibility  of  using  such 

t 

spectra  for  cross  eection  determinations.  The  spectra  obtained  by 


computer  eimulatiqaysre  used  to  calculate  the  photonuclear  croee  section 
of  the  »2p)  Ha  reaction  from  25  to  33  MeV.  The  calculated  crose 


section  ie  reported  and  ie  compared  to  a  previously  reported  cross  section 
and  to  a  cross  eection  calculated  using  thin  target  bremsstrahlung  spectra. 
Analyse e  of  the  results  obtained  indicate  that  thick  target  x-ray  spectra 
can  be  used  for  photonuclear  reaction  cross  section  determination.^ 
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A  STUDY  OF  THICK-TARGET  X-RAY  SPECTRA 


USING  PHOTONUCLEAR  REACTIONS 


COPTER  I 


INTRODUCTION 


High -energy  linear  accelerators  ere  becoming  incrtaeingly  sore 
popular  in  radiation  therapy  departments  because  of  the  augmented  treat¬ 
ment  flexibility  then  provide  radiation  oncologists.  The  depth-dose 
characteristics  o£  the  photon  beams  produced  by  these  accelerators  maxi¬ 
mise  the  dose  deliverable  to  treatment  volumes  located  at  substantial 
depths  vhila  minimising  the  dose  to  overlying  healthy  tissues.  In  addi¬ 
tion,  the  rapid  fall-off  of  the  depth  dose  of  the  electron  beams  they 
may  provide  permits  treet.jcnt  of  tumors  located  at  alight  depths  while 
sparing  underlying  structures.  The  nominal  energies  of  many  high-energy 
accelerators  exceed  the  threshold  energy  for  photonudeer  reactions. 
Thus,  the  study  of  photonudeer  reactions  la  becoming  an  Increasingly 
mors  Important  part  of  the  medical  radiological  sciences. 

Photonudeer  reactions  have  been  applied  in  many  different  ways 
in  medical  radiation  physics.  The  radioactivity  induced  in  radiation 
therapy  patients  has  been  used  to  investigate  the  distribution  of  dose 
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in  the  irradiated  volume  (1,2).  Photonuclear  activation  analysis  has 
been  employed  to  detect  trace  amounts  of  certain  elements  in  biological 
materials  (3,4).  Photoactivation  has  been  used  to  determine  megavdtage 
x-ray  spectra  (5) .  Photonuclear  reaction  thresholds  and  cross-section 
resonances  have  been  employed  in  the  energy  calibration  of  accelerators 
(6-8) .  Ratios  of  photoactivation  yields  in  pairs  of  metallic  foils  have 
been  used  to  specify  x-ray  quality  (9,10).  Photoactivation  has  been 
used  also  to  produce  clinically  useful  amounts  of  short-lived  positron- 
emitting  isotopes  used  in  nuclear  medicine  (11,12) . 

The  research  reported  here  utilized  photonuclear  reactions  to 
study  the  bremsstrahlung  spectra  produced  by  a  high-energy  clinical 
linear  accelerator.  The  study  consisted  of  essentially  three  parts. 
First,  ratios  of  photonuclear  yields  that  describe  the  distribution  in 
*energy  of  photons  in  an  x-ray  beam  were  investigated.  Second,  a ^numeri¬ 
cal  analysis  computer  calculation  of  thick-target  bremsstrahlung  spectra 

was  made  to  simulate  the  spectra  of  x  rays  produced  in  a  clinical  linear 

27  25 

accelerator.  Third,  the  yield  of  the  Al(y,2p)  Na  reaction  was 
measured  and  the  cross  section  determined  to  investigate  the  feasibility 
of  using  thick-target  spectra  for  cross-section  determination  and  to 
report  on  resonance  structure  in  the  cross  section  which  previously  has 
been  unreported. 


Overview 

This  report  has  been  divided  into  five  chapters.  In  Chapter  I 
each  of  the  three  parts  of  the  research  are  introduced  and  the  objectives 
of  each  part  are  set  forth.  In  Chapter  II  the  characteristics  of  photo¬ 
nuclear  reactions  are  described,  the  measurement  and  analysis  of  photo- 
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nuclear  yields  is  diacussed,  and  Che  equations  used  in  subsequent  calcu- 
la t ions  are  derived.  Chapter  III  describes  the  materials  used  in  the 
research  and  delineates  the  procedures  followed  in  each  part  of  the  work. 
Chapter  IV  presents  the  results  of  each  of  the  three  parts  of  the  research. 
In  Chapter  V  tha  conclusions  made  from  tha  resulta  obtained  are  aet  forth. 

Photonudear  Ratios 

X-ray  quality  is  a  term  esaentially  used  to  describe  the  distri¬ 
bution  in  energy  of  photons  in  a  bremsstrahlung  beam.  A  statement  of 
quality  is  useful  to  predict  the  p*nstration  of  the  beam,  to  calculate 
differential  absorption  in  various  materials,  and  to  attempt  to  explain 
variations  in  biological  effectiveness  (13) .  The  most  complete  specifi¬ 
cation  of  quality  is  provided  by  the  apectral  distribution  of  fluence, 
or  enargy  fluence,  as  a  function  of  photon  energy.  The  energy  spectrum, 
however,  is  extremely  difficult  to  measure,  thus  quality  is  commonly 
expressed  in  terms  of  certain  parameters  which  are  more  simply  measured. 

For  a  given  waveform,  the  most  obvious  parameter  which  may  be 
used  to  dsscribe  x-ray  quslity  is  the  x-ray  generator's  accslarating 
potential.  This  potential  detarmlnes  tha  final  enargy  of  the  electrons 
producing  the  x  rays.  The  final  electron  energy  is  taken  to  correspond 
to  the  maximum  energy  of  the  photons  in  the  x-ray  spectrum  and  is  thus 
termed  the  endpoint  energy.  The  specif icatiou  of  the  endpoint  energy 
alone  of  a  bremsstrahlung  spectrum  is  not  sufficient  to  characterize  Its 
spectral  quality.  X-ray  beams  having  tha  same  endpoint  energy  may  have 
quite  different  energy  characteristics  depending  upon  the  thickness  of 
the  bremsstrahlung  target  and  upon  the  amount  of  additional  filtration 


in  the  beam. 


Parameters  that  describe  the  penetrating  ability  of  an  x-ray 
spectrum  have  been  used  for  many  yeara  to  further  characterize  apectral 
quality.  The  most  prominent  of  these  parameters  la  the  Half  Value  Layer 
(HVL) ,  tha  thickness  of  a  material  that  will  reduce  tha  intensity  of  an 
x-ray  beam  to  one-half  the  original  Intensity.  The  HVL  is  sensitive  to 
not  only  tha  endpoint  energy  of  tha  spectrum,  but  elso  to  the  amount  of 
filtration  In  the  beam.  For  many  low-energy  (<  1  MeV)  applications,  an 
appropriate  specification  of  quality  Is  provided  by  the  two  parametera 
endpoint  energy  and  HVL  (13) .  The  quality  of  x-ray  beams  used  In  radia¬ 
tion  therapy  are  often  additionally  apeclfied  in  terms  of  depth  dose  in 
an  absorber.  The  specification  edopted  may  take  the  form  of  a  complete 
depth-dose  curve  from  the  surface  of  the  absorber  to  the  depth  at  which 
the  dose  la  reduced  to  a  certain  value,  say  102;  or  It  may  ba  restricted 
to  the  dapth-dose  at  one  or  two  aelected  depths. 

Unfortunately,  for  x-ray  energies  In  tha  megavoltaga  range,  the 
parameters  that  describe  x-ray  penetration  ability  (HVL  and  depth-doae) 
are  relatively  Insensitive  to  changes  In  quality.  Thia  Is  a  reault  of 
the  slow  changa  In  tha  attenuation  coefficients  of  ell  matarials  with 
increasing  photon  enargy  above  about  3  MeV.  Neth  and  Schulz  (14)  have 
calculated  the  percent  change  In  HVL  per  MeV  In  different  materials  for 
thln-target  apactra  of  endpoint  enargiea  from  5-40  MeV.'  The  euthors 
found  that  In  aven  the  mout  sensitive  material  (water),  the  average  sen¬ 
sitivity  (AHVL/AE)  was  found  to  change  less  than  2Z  per  MeV  for  endpoint 
energies  between  10  and  40  MeV.  The  ratio  of  the  percent  depth  dose  In 
tissue  at  10  cm  depth  to  that  at  2  cm  depth  changes  by  leaa  than  21  of 
the  mean  per  MeV  for  betatron  x  rays,  100  cm  SSD,  10  x  10  cm  field,  from 
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15-35  MeV  (15 , 16) .  The  ratio  of  ionization  measurements  made  at  5  cm 
depth  and  15  cm  depth  in  water  changes  by  only  approximately  0.5%  ot  the 
mean  when  Irradiated  with  15  to  30  MeV  x  rays  from  linear  accelerators 
(17). 

It  is  thus  apparent  that  a  more  sensitive  index  of  megavolt,  ge 
radiation  quality  is  required.  Comparison  of  radiation  therapy  treatment 
regimes  requires  a  complete  specification  of  the  characteristics  of  the 
radiation  beams  used.  A  sensitive  index  of  quality  would  aid  in  deciding 
whether  dosimetry  data  (such  as  percent  depth-dose)  for  one  x-ray  genera¬ 
tor  would  be  applicable  to  another.  Similarly,  the  suitability  of  energy- 
dependent  factors  necessary  for  dose  computations,  such  as  instrument 
calibration  factors  and  absorbed  dose  conversion  factors,  could  be  estab¬ 
lished  by  comparing  the  spectral  quality  index  of  the  calibration  machine 
with  that  of  the  user's  machine.  A  sensitive  index  could  also  be  used  as 
a  quality  control  device  to  assess  the  stability  of  radiation  therapy 
machines  and  signal  equipment  faults  before  they  become  obvious  by  other 
means.  Additionally,  if  the  detector  used  to  obtain  such  an  index  ware 
small  enough,  it  would  be  possible  to  measure  the  variations  in  spectral 
quality  at  selected  locations  in  irradiated  media. 

A  sensitive  index  of  x-ray  quality  in  the  megavoltage  energy 
range  is  provided  by  the  ratio  of  photonuclear  yields  in  carefully  aelected 
materials.  Photonuclear  ratios  are  based  on  the  fact  that  photonuclear 
reactions  in  most  elements  result  in  radioactive  nuclei  whose  decay  rate 
can  be  measured  to  determine  photonuclear  yields.  If  the  elements  are 
chosen  such  that  the  photonuclear  cross  sections  peak  at  different  loca¬ 
tions  in  the  photon  energy  spectrum,  the  ratio  of  the  photonuclear  yields 
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produced  in  each  element  will  vary  quite  sensitively  with  energy.  Photo- 
nudaar  yield  ratios  for  thin-target  bremsstrahlung  spectra  from  15  to  40 
MeV  have  bean  measured  and  published  by  Nath  and  Schulz  (9).  Their  method 
Involved  counting  the  activity  induced  in  each  of  a  pair  of  metallic  foils. 
Their  experiment  produced  uniqua  "Photoactivation  Ratios'1  from  about  25 
to  35  MeV.  The  ratio  of  tha  photonuclear  yield  of  tha  ®^Y(y,2n)®^Y  reac¬ 
tion  to  the  photonuclear  yield  of  tha  ®^Cu(Y»n)^Cu  reaction  yielded  the 
most  sensitlva  quality  index  as  a  function  of  photon  energy  in  the  25-35 
MeV  range.  The  maximum  sensitivity  was  approximately  25Z  per  MeV ,  a 
definite  improvement  over  the  indices  of  quality  previously  mentioned. 

The  purpose  of  the  first  part  of  this  work  was  to  seek  photo¬ 
nuclear  yield  ratios  sensitlva  to  x-ray  spactra  produced  by  thick  targets 
in  the  15  to  25  MeV  energy  ranga  whera  such  ratios  are  lacking.  Many 
high-energy  medical  linear  accelerators  hava  nominal  energies  in  this 
range.  The  materials  were  chosen  such  that  the  photonuclear  reactions  in 
the  elements  producing  the  yields  used  to  compute  tha  ratios  resulted  in 
radioactiva  nuclei  that  emitted  the  same  energy  gamma  ray.  This  permits 
the  determination  of  ratios  without  prior  knowledge  of  the  efficiency  of 
the  radioactivity  counter.  The  samplas  consisted  of  chemical  compounds 
which  containad  the  pair  of  elements  from  which  photonuclear  yiald  ratios 
were  computed.  This  furthar  facilitates  ratio  determination  in  that 
self-absorption  factors  need  not  be  determined. 

ihe  ^Al(y.2p)^Na  Cross  Section 
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heve  been  described  in  some  detail  by  Bogdankevich  and  Nikoleev  (19) . 
Frequently,  a  photonudeer  reaction  results  in  the  formation  of  a  radio¬ 
active  nucleus  es  e  result  of  the  nuclear  transformations  Induced  by  the 
high-energy  photons.  The  activity  Induced  in  e  sample  of  materiel  exposed 
to  a  flux  of  high-energy  photons  is  directly  related  to  the  cross  section 
of  the  particular  photonudeer  reaction  producing  the  redioective  species 
and  to  the  number  and  energy  of  the  photons  traversing  the  sample.  By 
measuring  the  activity  of  the  sample  after  irredietion,  it  is  possible  to 
obtain  the  cross  section  of  the  reaction. 

Most  commonly,  the  source  of  high-energy  photons  used  in  photo- 
nudear  research  is  the  bremsstrahlung  beam  of  either  betatrons  or  linear 
accelerators.  Bremsstrahlung  beams  possess  e  spectrum  of  photon  energies 
ranging  from  essentially  zero  to  the  maximum  kinetic  energy,  or  endpoint 
Energy,  of  the  electrons  producing  the  bremsstrahlung;  thus,  photonudeer 
reactions  produced  by  bremsstrahlung  x  rays  are  induced  by  many  photons 
in  the  energy  epectrum.  On  the  other  hand,  the  number  of  bremsstrahlung- 
induced  reactions  et  a  given  endpoint  energy  is  related  to  the  photo- 
nuclear  cross  section  as  e  function  of  energy.  If  the  reactione  produced 
result  in  radioactive  nuclei,  then  ve  may  define  a  photonuclesr  yield  at 
a  given  endpoint  energy,  proportional  to  the  number  of  photonuclear  reac¬ 
tions,  which  is  determined  from  e  measurement  of  the  redioectlvlty  induced. 
If  the  bremsstrahlung  endpoint  energy  is  varied  end  yield  measurements 
are  made  at  each  endpoint  energy,  then  a  yield  curve,  e  graphic  represen¬ 
tation  of  the  photonudeer  yield  as  a  function  of  bremsstrahlung  endpoint 
energy,  may  be  constructed.  If  the  bremsstrahlung  spectrum  is  known, 
the  photonudeer  cross  section  es  e  function  of  energy  may  then  be 


extracted  mathematically  from  the  yield  curve* 

Implicit  in  the  previous  discussion  is  the  dependence  of  the 
cross  section  obtained  upon  the  mathematical  model  chosen  to  describe  the 
bremsstrahlung  spectrum.  Photonuclear  experiments  performed  in  the  pest 
have  employed  the  Schlff  integreted-over-angle  spectrum  (20)  to  describe 
the  photon  energy  distribution  of  the  bremsstrehlung  beam  (18,19).  The 
Schiff  calculations  were  performed  under  the  assumption  thet  the  electron 
interacts  only  once  in  the  target,  i.e.,  the  target  is  infinitely  thin. 
Thus,  in  e  strict  sense,  the  use  of  Schiff  spectre  is  justified  only  in 
experiments  utilizing  thin  bremsstrahlung  targets.  In  en  experiment  where 
photonucleer  reactions  are  produced  by  thick- target  bremsstrahlung  (es  is 
the  case  in  this  experiment),  the  bremsstrahlung  spectra  must  be  corrected 
for  the  spectre!  distortions  introduced  by  electron  collisionel  losses  as 
well  es  by  photon  attenuation  processes.  Only  then  can  e  realistic  cross 
section  be  determined  from  such  an  experiment. 

In  the  second  portion  of  the  research  presented  here,  the 
27  25 

t  (y,2p)  Na  reaction  yield  wes  measured  in  0.25  MeV  increments  from 
threshold  to  33  MeV  and  the  reectlon  cross  section  was  calculated.  This 
cross  section  hes  never  been  determined  with  such  energy  resolution  before 
although  the  cross  section  has  been  reported  previously  (21).  The 
increased  energy  resolution  of  this  experiment  enables  a  more  detailed 
examination  of  resonance  structure  in  the  cross  section  possibly  indi¬ 
cating  the  existence  of  higher  energy  nuclear  stetes. 

The  Bremsstrehlung  Spectrum 

Photonuclear  cross  sections  deduced  from  measurements  of  yield 
depend  to  a  very  large  degree  on  the  bremsstrahlung  spectrum  utilized  in 
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Che  analysis.  The  review  of  phoeonudear  literature  conducted  in  prepa¬ 
ration  for  this  work  conveyed  the  fact  that  thick-target  bremsstrahlung 
spectra  had  never  before  been  used  in  photonudear  reaction  cross  section 
determinations  regardless  of  the  target  thickness  used  in  the  experiment. 
To  examine  the  feasibility  of  using  thick-target  spectra  for  cross  sec¬ 
tion  determinations  t  in  the  third  part  of  this  research  a  numerical  analy¬ 
sis  computer  calculation  of  thick-target  bremsstrahlung  spectra  was  per- 

27  25 

formed  and  the  spectra  obtained  were  used  to  calculate  the  Al(y,2p)  Ha 
reaction  cross  section.  The  cross  section  was  also  calculated  using 
thin-target  spectra  to  determine  the  effect  of  spectral  shape  on  the  cross 


section  obtained. 


CHAPTER  II 


THEORY 

Photonudcar  Reactions 

A  photomidear  reaction,  or  nuclear  photoeffect,  is  on  electro¬ 
magnetic  Interaction  In  which  a  photon  Is  completely  absorbed  by  a  nuc¬ 
leus.  The  result  Is  the  emission  cf  one  or  more  nuclear  constitr ''\*:s . 

A  reaction  of  this  type  may  be  represented  by 

X  +  Y  b  +  Y  or  X(y,b)Y  (1) 

where  X  Is  the  target  nucleus,  y  represents  the  Incoming  photon,  b  Is  the 
emitted  particle,  and  Y  Is  the  residual  nucleus.  Frequently  the  defini¬ 
te  of  a  photonuclear  reaction  Is  restricted  to  Include  only  those  pro¬ 
cesses  In  which  the  emitted  particles  are  either  protons  (y,p),  (y»2p) » 
etc.,  neutrons  (Y,n),  <Yi2n),  etc.,  or  aggregates  of  protons  and  neutrons 
(Y»np),  (Y»2np),  (y»«)»  etc.  Thus  photomason  production  (y»«)  and  nuclear 
resonance  scattering  (y,y')  *te  excluded. 

Photonuclear  reactions  are  endoerglc,  that  Is,  energy  Is  required 
to  Initiate  the  reaction.  The  minimum  amount  of  energy  which  can  initiate 
the  reaction  la  called  the  threshold  energy.  A  photonuclear  reaction  Is 
allowed  energetically  only  when  the  photon  energy  equals  or  exceeds  the 
reaction  threshold  energy.  A  photonuclear  reaction's  threshold  energy 
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can  be  obtained  from  equations  of  energy  and  momentum  conservation.  For 
a  reaction  X(y,b)Y,  the  threshold  energy  will  be  given  by 

k«i0 '  BVl  +  <W2)1  <2> 


where  is  the  mass  of  the  target  nucleus ,  c  is  the  velocity  of  light, 
and  BE^,  the  binding  energy  of  the  emitted  particle,  ia  given  by 


BEb  *  SEb  +  Ey 


where  SE^  ia  the  emitted  particle's  separation  energy  and  ia  tha 
excitation  energy  of  tha  residual  nucleus.  Tha  separation  energy  SE^ 
may  be  calculated  from  the  expression 


S!^  -  +  AEy  -  BEx  (4) 

where  AE^,  ABy,  and  AEX  are  called  the  mass  decrements  of  the  emitted 

particle,  residual  nucleus,  and  target  nucleus,  respectively.  Typically, 

threshold  energies  of  photonuclesr  reactions  exceed  8  MeV. 

Tha  probability  or  relative  frequency  with  which  a  photon  of  a 

determined  energy  will  undergo  a  specific  photoouclear  absorption  process 

in  a  certain  absorber  ia  Indicated  by  that  reaction's  "cross  section," 

2 

o(cm  /nucleus).  Tha  cross  section  of  the  photonuclear  reaction  X(y,b)Y 
may  be  defined  as  the  ratio  of  the  number  par  sec  of  photons  y  of  anergy 
k  absorbed  by  the  nucleus  to  tha  photos  flux  at  that  anergy. 


(5) 
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Thus  a  particular  photonudear  reaction  will  have  a  specific  probability 

or  cross  section  at  a  specific  Incident  photon  energy.  The  cross  section 

has  dimensions  of  length  squared  with  the  ad  hoc  unit  of  barns  (1  barn  ■ 
-24  2 

10  cm  )  being  In  common  usage. 

Above  the  threshold  energy  the  cross  section  Increases  with 
Increasing  Incident  photon  energy  k  for  several  MeV,  reaches  a  maximum 
value,  then  decreases  with  further  increases  of  photon  energy.  The  shape 
of  this  rather  broad  peak  la  very  characteristic  of  reactions  of  this 
type  and  Is  called  the  "giant  resonance."  Photon  absorption  In  this 
region  Is  almost  exclusively  by  electric  dipole  Interaction  with  only  a 
small  (6X)  contribution  by  electric  quadrupole  transitions  (22) .  The*' 
phenomena  are  further  discussed  subsequently  In  this  section.  The  general 
characteristics  of  the  giant  resonance  have  been  set  forth  by  numerous 
— authors,  among  them  Strauch  (18),  Levinger  (23),  Wilkinson  (24,25)  ,  and 
Hayward  (26).  In  general  the  giant  resonance  may  peak  at  about  0.1  barns 
within  a  region  from  10  to  30  MeV  with  a  full  width  at  half  maximum  (I*) 
of  approximately  3-8  MeV  (18,23,26).  The  energy  E^  at  which  the  cross 
section  Is  a  maximum  Is  approximately  20  MeV  for  light  elements  (23,26). 
For  heavy  elements,  E^  varies  inveraely  with  the  element'a  mass  number 
A,  with  dependences  of  A”®*^  (23)  and  A”®*^  (26)  having  been  observed. 

The  area  under  the  o(k)  curve,  called  the  Integrated  cross  section 
end  given  by 

-  J o(k)dk  (6) 

Is  proportional  to  A  with  proportionality  constant  about  0.02  MeV-barns 
(23). 


Goldhaber  and  Teller  have  explained  the  giant  resonance  In  terms 
of  a  collective  model  of  the  nucleus  (27).  According  to  their  explana¬ 
tion,  the  giant  resonance  excitation  Is  one  in  which  all  protons  collec¬ 
tively  execute  harmonic  motion  relative  to  all  neutrons.  The  protons 
and  neutrons  may  be  regarded  as  two  Incompressible,  Interpenetrating 
fluids  which  oscillate  back  and  forth  relative  to  each  other.  The  reeo- 
nant  frequency,  and  hence  the  quantum  energy  -fiio,  will  depend  on  the  nuc¬ 
lear  size.  Their  calculations  yield  a  resonance,  the  energy  of  which 
-0.17 

varies  as  A  ,  In  rough  eccordance  with  experimental  findings. 

Although  the  collective  model  exactly  exhausts  the  dipole  sum 
rule  (see  section  on  Theories  of  Photonudear  Reections)  and  has  been 
moderately  successful  In  predicting  the  experimentally  observed  trend  of 
Ea  is  a  function  of  A,  It  has  Its  shortcomings.  It  hss  nothing  to  eay 
about  T,  the  width  of  the  resonance,  and  provides  no  explonetion  for  the 
anomalous  emission  of  protons  from  heevy  nuclei.  The  probability  of  this 

4 

process  Is  a  factor  of  xO  times  more  corancn  than  one  would  expect  when 
the  euergy  of  exdtstlon  Is  shared  among  all  of  tbe  nucleons  end  when 
all  have  equal  stetlsticel  probability  of  emission.  In  view  of  theee 
difficulties  It  has  been  tempting  to  seek  an  alternative  explanation. 

Such  an  explanetlon  Is  efforded  by  the  shell  model. 

Detailed  inspection  of  the  giant  resonance  of  photonudear 
reections  In  light  elements  reveals  crose  section  maxima  and  minima,  or 
"structure which  Imply  increased  photonucleer  ebeorptlon  et  photon 
energies  corresponding  to  excitation  energies  of  pertlcular  nuclear 
states.  This  structure  In  the  glent  resonance  Is  best  explained  In  terms 
of  ellowed  transitions  between  nucleon  energy  states  which  ere  obtained 
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through  the  Independent  particle  or  shell  model  of  the  nucleus.  The 
shell  model  of  nuclear  structure  Is  discussed  in  some  detail  in  Appendix 
A.  Allowed  transitions  between  energy  states  ere  treated  in  Appendix  B. 
In  the  next  section  photonuclear  reaction  theories  will  be  discussed 
with  perticular  emphasis  on  the  Independent  particle  model  of  nuclear 
etructure.  Throughout,  existing  theories  arc  compared  to  evallable 
experimental  findings. 

Theories  of  Photonuclear  Reactions 

As  described  in  Appendix  B,  transitions  between  states  are  pro¬ 
portional  to  the  square  of  the  transition  matrix  element.  In  the  inter¬ 
action  of  photons  with  bound  electrons,  es  in  the  atomic  photoeffect, 
the  perturbation  operator  in  the  matrix  element  overlap  lntegrel  is  the 
electric  dipole  operator.  The  matrix  clement  for  electric  dipole  transi¬ 
tions  Mg^  Is  given  by  (28) 

“ei  ■  “"V  /’*  \  *iYT  <7) 

where  z  is  the  component  of  the  displecement  along  the  electric  vector 
of  the  electromagnetic  field,  is  the  energy  difference 

batveen  the  excited  and  ground  states  divided  by  and^^  z^  la  the 
dipole  moment  of  the  atom  found  by  summing  the  component  of  dieplaceaent 
over  all  electrons  of  the  atom. 

In  the  nuclear  photo«ffert,  where  phnton  interactions  with  pro¬ 
tons  and  neutrons  are  treated,  two  major  changes  must  be  introduced: 

1.  An  effective  charge  must  be  used  of  aN/A  for  protene,  and 
-eZ/A  for  neutrons,  where  A  is  the  mivss  number,  Z  is  the  etoralc  number. 


and  N  =  A  -  Z  Is  the  neutron  number.  The  effective  charge  Is  Introduced 
since  the  displacements  of  the  nucleons  should  be  meesured  from  the 
nuclear  center  of  mass. 

2.  The  partly  exchange  (Major ana)  character  of  the  neutron- 
proton  potential  must  be  considered.  Sieger;.' s  theorem  states  that 
exchange  forces  between  neutrons  and  protons  imply  the  existence  of  a 
current  of  chergad  particles  (ir  mesons)  between  the  nucleons.  Thus,  the 
complete  calculation  of  the  dipole  matrix  element  should  Involve  the  sum 
of  nucleon  and  * -meson  currants. 

Two  methods  have  been  used  for  theoretical  calculations  of  the 
photon  absorption  cross  section  for  nuclei.  First,  the  dipole  matrix 
element  can  be  calculated  for  absorption  of  phccons  If  explicit  assump¬ 
tions  are  made  as  to  the  ground  etate  and  excited  state  wave  functions 
— (for  various  energies).  Second,  use  can  be  made  of  "sum- rules”  In  which 
transitions  from  the  ground  state  are  summed  over  all  possible  final 
states.  This  method  involves  knowledge  only  of  the  transition  operators 
end  of  the  wavs  function  for  the  ground  etate.  Sines  the  ground  stats 
wave  function  is  a  solution  of  Schrodlnger's  equation  and  thus  determines 
the  potential,  Schrodlnger's  equation  could  then  In  principle  be  solved 
for  this  potential  for  the  wave  functions  of  the  various  excited  states. 
The  first  method  la  then  used  for  calculations. 

If  the  classical  Thoaas-Rslehe-Kuhn  (TRit)  sum-rule  for  electric 
dipole  absorption  of  photons  by  Z  electrons  In  an  atom  Is  modified  for 
the  nuclear  pbotoefiect  by  the  two  effects  discussed  above  (the  use  of 
effective  charges  and  the  significance  of  Majorana  exchange  forces),  an 
expression  Is  obtained  for  the  Integrated  cross  section  for  electric 


dipole  absorption  of  photons  by  a  nucleus  of  atomic  number  Z,  mass  number 
A,  and  neutron  number  N.  This  dipole  absorption  "sum- rule"  is  given  by 
(22) 

olftt  -  Jo( k)dk  -  0.058  ~  (i  +  0.8X)MeV-bams  (8) 

where  X  is  the  fraction  of  the  neutron-proton  force  that  has  an  exchange 
character.  The  upper  limit  of  integration  must  be  set  at  about  the  meson 
threshold  since  Siegert's  theorem,  on  which  the  sum-rule  is  based,  will 
be  invalid  beyond  that  (24).  Values  of  X  from  0.5  to  0.7  have  been  found 
from  neutron-proton  scattering  experiments.  The  number  0.8  was  calcu¬ 
lated  using  a  nuclear  model  of  a  degenerate  Fermi  gas;  however,  the  value 
of  o^nt  is  relatively  independent  of  the  nuclear  model  except  for  the 
value  of  the  coefficient  0.8  and  possible  small  terms  (23). 

The  Independent  particle  shell  model  (IFM)  of  the  nucleus  has 
been  very  successful  in  describing  dipole  transitions  in  light  elements. 
Nuclear  photoeffect  calculations  have  been  performed  by  various  authors 
for  harmonic  oscillator,  infinite  square  well  and  finite  square  well 
potentials  C 3).  Wilkinson  (24,25)  hss  shown  that  essentially  all  of 
the  dipole  sum-rule  is  exhausted  in  transitions  between  the  last  filled 
shell  levels  to  virtusl  shell  levels  in  the  continuum.  Transitions  from 
unfilled  shells  do  not  contribute  appreciably.  The  most  important  tran¬ 
sitions  ware  found  to  be  those  between  states  whose  principal  quantum 
number  and  increase  in  orbital  angular  momentum  is  unity.  In  addition, 
his  calculations  showed  that  the  strongest  transitions  are  those  between 
states  of  larger  orbital  angular  momentum.  The  well-known  angular  momen¬ 
tum,  parity,  and  isotopic  spin  selection  rules  (see  Appendix  B)  are 
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assumed  to  hold.  Since  all  shells  below  the  uppermost  filled  shell  are 
also  filled,  transitions  between  these  states  are  forbidden  by  the  Fauli 
exclusion  principle  and  excitations  between  these  shells  can  contribute 
no  oscillator  strength  (transition  probability)  to  the  total  sum;  in  the 
same  way,  downward  transitions  which  contribute  to  the  sum  negatively 
are  also  forbidden.  Thus,  oscillator  strength  is  "transferred  upward" 
through  the  shells  and  the  oscillator  strength  from  the  uppermost  filled 
shell  results  from  a  collective  absorption  by  all  the  nucleons  in  the 
shell  with  the  excitation  of  one  nucleon  to  the  next  allowable  unfilled 
level  (29).  Wilkinson  showed  that  the  Kl  transitions  between  last  filled 
and  virtual  levels  which  are  allowed  by  the  selection  rules  are  energeti¬ 
cally  clustered  sufficiently  close  to  blend  into  an  apparent  giant  reso¬ 
nance.  In  addition,  Wilkinson's  I PM  calculations  have  been  able  to 
account  for  a  roughly  correct  oscillator  strength  and  for  reasonable 
resonance  widths. 

In  spite  of  the  success  of  Wilkinson's  IPM  calculations  in 
explaining  the  nuclear  photoeffect  phenomena  described  above,  the  giant 
resonance  energy  E^  resulting  from  the  calculations  was  too  low.  The 
problem  was  in  his  choice  of  a  static  potential  whose  depth  was  deter¬ 
mined  by  nucleon  binding  energies  of  the  nuclei  treated  (29) .  He  sug¬ 
gested  that  use  of  a  velocity-dependent  shell  model  potential  that  does 
not  coonute  with  nucleon  position  and  that  considers  a  nucleon-reduced 
mass  might  result  in  increased  estimates  of  E  in  better  accord  with 
experiments . 

Further  IFM  calculations  were  performed  by  Elliott  and  Flowers 
(30),  who  calculated  the  energies  and  radiative  widths  for  transitions 
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from  the  odd  parity  levels  of  ■"  0,  and  by  Brown  and  Bolsterli  (31)  and 
Brown,  Castillejo  and  Evans  (32)  who  considered  the  effect  of  particle- 
hole  Interactions  on  the  El  energy  levels  in  the  nuclear  photoeffect. 

The  calculations  took  Into  consideration  spin-orbit  coupling.  When  a 
nucleon  Is  elevated  by  El  photon  absorption  to  a  virtual  shell  level,  a 
hole  Is  left  In  the  filled  shell.  Such  holes  can  be  treated  in  many  ways 
like  single  particles  In  the  shell.  Since  the  process  of  hole  formation 
Is  dipole;  the  excited  particle  and  hole  are  strongly  correlated  In 
angle;  that  Is,  *:heir  angular  momentum  must  be  coupled  to  form  a  1- 
state,  assuming  the  original  to  be  In  a  0+  state  (see  Appendices  A  and  B) . 
Because  many  particle-hole  states  can  be  formed  and  because  these  states 
are  almost  degenerate  In  energy,  the  particle-hole  Interaction  can  have 
a  profound  effect  In  redistributing  dipole  transition  strength.  The 
above  authors'  calculations  using  particle-hole  Interactions  resulted  In 
Increased  excitation  energies.  In  better  agreement  with  experimental 
results.  Although  a  specific  mixture  of  exchange  forces  (Majorana  and 
Bartlett)  was  used,  the  general  trend  of  levels  obtained  was  found  to  be 
relatively  Insensitive  to  the  exchange  mixture  chosen.  The  calculations 
Indicated  that  the  entire  dipole  oscillator  strength  will  be  exhausted 
at  photon  energies  below  30  MeV.  Figure  1  Illustrates  particle-hole 
transitions  for  the  ^0  nucleus  (29)  showing  the  simple  harmonic  oscil¬ 
lator  shell  levels  (a)  with  flu  ■  17  MeV,  and  the  same  levels  with  spin- 
orbit  coupling  (b).  The  transitions  predicted  by  Wilkinson  are  shown 
with  vertical  lines  along  with  the  transition  energies. 

Although  the  collective  and  the  IPM  models  may  seem  very  differ¬ 
ent,  they  make  very  similar  predictions  In  many  cases.  Brink  (33)  has 
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demonstrated  that  the  similarities  probably  are  not  accidental  and  that 
the  two  models  are  basically  the  same.  The  particle-hole  or  nucleon- 
nucleon  correlations  necessary  to  correct  the  independent  particle  model 
giant  resonance  predictions  play  an  important  collectivizing  role  in  the 
theory.  The  pure  IPM  account  has  neglected  the  effects  of  interactions 
between  the  many  superposed  transitions.  These  interactions  may  very 
well  be  strong  and  the  transitions  may  be  pulled  together  by  them  to  show 
a  tighter  grouping  than  expected  by  the  raw  IPM.  It  appears  that  intro¬ 
duction  of  nucleon-nucleon  correlations  into  either  model  is  necessary  if 
the  model  is  to  describe  anything  more  than  the  most  gross  feature  of  the 
photonuclear  reaction  cross  saction.  A  decision  regarding  the  validity 
of  one  modal  or  another  is  unimportant.  Use  is  made  of  the  model  that 
yields  the  most  accurate  predictions. 

■*» 

Analysis  and  Solution  of  Photonuclear  Yield  Functions 

Analysis 

The  calculation  of  cross  sections  from  measured  yield  curves  is 
the  subject  of  this  section.  The  analysis  presented  was  formulated  by 
Penfold  and  Leiss  (34) .  In  a  photonuclear  experiment  a  sample  is  exposed 
to  a  bremsst rah lung  beam.  Tha  number  of  photonuclear  reactions  occurring 
within  a  sample  depends  upon  the  availability  of  target  atoms,  on  tha 
cross  section  for  the  particular  reaction,  and  on  the  flux  of  photons 
Incident  upon  the  sample.  The  photon  flux  is  measured  by  a  calibrated 

ionization  chamber  which  is  irradiated  simultaneously  with  the  sample. 

2 

A  sample  that  contains  n^  target  nuclei  per  cm  is  exposed  to  a  beam  of 
maximum  energy  T  where  M(T,k)  photons  of  energy  between  k  and  k  +  dk. 
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per  unit  range  of  k  and  per  unit  monitor  response,  enter  the  sample.  If 
o(k)  is  the  cross  section  in  cm  per  nucleus  as  a  function  of  energy  for 
the  reaction  under  investigation,  then  the  normalized  yield  of  photonuc- 
lear  reactions  at  a  given  breosstrahlung  endpoint  energy  per  unit  monitor 
response,  Yq(T),  is  given  by  the  integral 

Yn(T)  “  nt  I/  M(T,k)o(k)dk  (9) 

The  Integral  equation  of  Eq.  (9)  is  called  the  photonudear  yield  func¬ 
tion.  Repeated  measurements  of  yield  at  different  bremsstrahlung  end¬ 
point  energies  and  a  knowledge  of  the  bremsstrahlung  spectrum  allow 
deduction  of  the  cross  section. 

The  photonudear  yield  function,  or  simply  photonudear  yield, 

may  be  thought  of  as  an  overlap  Integral  as  shown  in  Figure  2.  The 

amount  of  yield  Y(T)  will  depend  upon  the  degree  to  which  the  endpoint 

energy  T  of  the  bremsstrahlung  spectrum  M(T,k)  exceeds  the  threshold 

energy  E  of  the  photonudear  cross  section  o(k).  In  addition,  the 
tn 

amount  of  yield  depends  on  the  relative  shapes  of  the  functions  M(T,k) 
and  o(k). 

The  spectrum  M(T,k)  is  the  bremsstrahlung  photon  flux  that 
enters  the  sample  per  unit  monitor  response.  This  spectrum  is  the  pro¬ 
duct  of  three  factors: 

1)  An  intrinsic  bremsstrahlung  spectrum  function  i(T,k),  which 
describes  the  x-ray  spectrum  emitted  by  a  "thin"  bremsstrahlung  radiator. 
♦(T,k)  is  proportional  to  the  bremsstrahlung  cross  section. 

2)  A  spectrum  modification  function  f  (k),  which  describes  the 

8 


Figure  2.  Photonuclear  yield  as  an  overlap  integral. 
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deformation  of  the  thin-target  spectrum.  The  function  of  f^(k)  accounts 
for  the  spectral  distortion  caused  by  electron  collisional  losses  and  by 
the  attenuation  of  photons  in  the  radiator  and  in  structures  which  may 
be  located  in  the  beam  between  the  radiator  and  the  sample. 

3)  A  monitor  response  function  F(T),  which  normalizes  the  spec¬ 
trum  that  enters  the  sample- to -unit  monitor  response. 

In  view  of  the  above,  the  spectrum  M(T,k)  of  Eq.  (9)  may  be 
expressed  as 


M(T,k) 


*(T,k)fs(k)  N(T,k) 
F(T)  "  F(T) 


(10) 


where  N(T,k)  is  the  actual  bremsstrahlung  spectrum  to  which  the  sample  la 
exposed.  The  spectrum  N(T,k)  can  be  estimated  by  a  numerical  method. 

For  the  purpose  of  the  present  discussion  it  is  assumed  that  N(T»k)  is 
known.  The  characteristics  of  this  spectrum  and  the  details  of  the 
numerical  method  used  for  its  estimation  are  presented  In  a  thorough 
fashion  in  the  next  section  of  this  chspter  (The  Bremsstrahlung  Spectrum) . 

In  an  analogous  fashion  one  may  define  a  spectrum  function 
Ma(T,k)  which  describes  the  photon  flux  entering  the  monitor  chamber  per 
unit  monitor  response.  If  the  geometry  of  the  experimental  setup  is 
such  that  the  photon  spectrum  entering  the  monitor  has  first  traversed 
the  sample,  then  a  spectrum  distortion  function  f^OO  may  be  defined 
which  describes  the  degradation  of  the  bremsstrahlung  spectrum  N(T,k)  by 
photon  absorption  in  the  sample.  Similarly  to  Eq.  (10),  Mm(T,k)  is 


M  (T,k) 

Dl 


N(T,k)fB(k) 

F(T) 


(ID 
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The  monitor  response  function  F(T)  normalizes  the  spectrum  inci¬ 
dent  on  the  monitor,  [N(T,k)fm(k) ] ,  to  unit  response  of  the  monitor.  In 
other  words,  F(T)  gives  the  total  number  of  photons  in  the  bremsst rah lung 
spectrum  that  yields  a  unit  count  in  the  calibrated  ionization  chamber. 
F(T)  is  defined  as 

F(T)  -  R(T)  f  f  (k)kN(T,k)dk  (12) 

I  0 
J  0 

where  R(T)  is  the  calibration  factor  of  the  monitor  in  units  of  response 
per  unit  incident  energy.  If  the  monitor  and  sample  dc  not  subtend  iden¬ 
tical  solid  angles  and  if  the  irradiation  geometry  is  not  perfectly 
reproducible  between  measurements  of  yield  for  different  bremsstrahlung 

endpoint  energies,  an  additional  factor  u>  (T)/u  (T)  must  be  introduced 

m  s 

into  Eq.  (12).  This  factor  is  a  ratio  of  the  angular  distribution  of  the 
radiation  integrated  over  the  monitor  solid  angle  to  that  integrated  over 
the  sample  solid  angle. 

The  spectrum  distortion  function  fffi(k)  describes  the  photon 

attenuation  due  to  interaction  processes  in  the  sample.  At  the  high 

photon  energies  employed  in  photonuclear  measurements,  the  predominant 

interaction  processes  are  pair  production  and  Compton  scattering.  The 

cross  sections  for  these  reactions  at  high  photon  energies  vary  very 

little  as  a  function  of  energy.  Therefore,  the  function  f  (k)  will  not 

m 

distort  the  spectrum's  shape  but  will  only  reduce  its  intensity  by  a 
constant  amount.  One  can  then  remove  the  function  f  (k)  from  the  inte- 

ED 

gral  and  replace  it  by  ft(T),  a  transmission  factor  as  a  function  of 
bremsstrahlung  endpoint  energy  which  represents  the  fraction  of  the  pho- 


tons  incident  on  the  sample  that  are  transmitted  through  the  sample.  In 
view  of  the  above,  the  monitor  response  function  F(T)  is  now 


where 


F(T)  -  R(T) f ^ (T)E(T) 


E(T) 


kN(T,k)dk 


(13) 


(13a) 


is  the  total  energy  in  the  bremsstrahlung  spectrum  of  endpoint  energy  T. 

A  solution  to  the  photonuclear  yield  function,  Eq.  (9),  may  now 
be  considered  since  the  spectrum  M(T,k)  has  been  defined  through  Eqs. 
(10)  and  (13).  Ihe  photonuclear  yield  function  now  has  the  form 


Y„m 


n 


t 


f(tV 


N(T,k)o(k)dk 


(14) 


Here,  the  spectrum  N(T,k)  is  the  actual  bremss trail  lung  spectrum  seen  by 
the  sample.  The  monitor  response  function  F(T) ,  which  is  defined  by  Eq. 
(13),  has  been  removed  from  the  integral  because  it  does  not  vary  as  a 
function  of  the  photon  energy  k.  The  remaining  variables  of  Eq.  (14) 
maintain  their  previous  definitions. 

The  photonuclear  yield  function  as  given  by  Eq.  (14)  is  depen¬ 
dent  upon  experimental  conditions,  in  particular  upon  the  radiation 
monitor  used.  Equation  (14)  may  be  rearranged  into  a  form  which  is  inde 
pendent  of  the  exact  nature  of  the  experimental  arrangement.  This 
rearranged  form  has  been  called  the  reduced  yield  function  by  Penfold 

and  Leiss  (34).  The  reduced  yield  Yf(T)  may  be  defined  as  the  photonuc- 

2 

lear  yield  (number  of  photonuclear  reactions)  per  target  atom  per  cm  . 
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It  is  obtained  by  multiplying  the  normalized  yield  (number  of  photonuclear 
reactions  per  unit  monitor  response)  ^(T)  times  the  monitor  response 

function,  F(T),  and  then  dividing  the  product  by  n^,  the  number  of  target 

„  2 
atoms  per  cm 

Y  (T)F(T) 

Y  (T)  -  -  •  (15) 

r  n 


Using  the  above  definition,  the  reduced  photonuclear  yield  func¬ 
tion  is 


Y  (T) 
r 


N(T,k)o(k)dk  . 


(16) 


It  is  to  the  reduced  yield  function  Y  (T)  that  solutions  are  obtained. 

r 

Note  that  the  reduced  yield  Y^(T)  is  related  to  the  experimentally  deter¬ 
mined  normalized  yield  Yr(T)  through  Eq.  (15)  and  that  the  monitor 
response  function  F(T)  can  be  obtained  with  Eq.  (13). 


Solutions 

The  solutions  of  the  reduced  yield  function  may  be  of  two  types: 
1)  solution  of  the  integral  equation  directly  or  2)  solution  of  a  set 
of  linear  equations  which  approximate  the  integral  equation.  Y^(T)  is 
not  known  as  a  continuous  function  of  endpoint  energy  T;  it  is  known 
only  at  a  finite  set  of  energies,  T^.  In  practice,  the  yields  Y(T^)  are 
measured  at  n  endpoint  energies  spaced  in  equal  intervals  AT.  There¬ 
fore,  a  condition  implicit  in  the  solution  of  the  integral  equation, 
namely  that  the  yield  Y^(T)  be  known  as  a  continuous  function  of  energy, 
is  not  satisfied.  Thus,  the  solutions  sought  here  will  be  of  the  second 
type,  solutions  of  linear  equations  which  approximate  the  integral. 
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The  reduced  yield  function  can  be  replaced  by  a  set  of  linear 
equations  if  a  "step-wise"  approximation  to  the  bremsstrahlung  spectrum 
is  made,  as  shown  in  Figure  3.  Let  be  the  number  of  photons  in  a 

bremsstrahlung  spectrum  N(T^,k)  of  endpoint  energy  that  possess  ener¬ 
gies  between  k^  -  AT  and  k^  where  AT  (also  called  an  energy  bin)  is  a 
photon  energy  interval  equal  to  the  spacing  of  yield  points  on  the  yield 
curve  , 


N’ 


ij 


1 

aT 


/ 


N(T,k)dk  . 


kj-AT 


(17) 


Each  depends  on  the  relative  position  of  the  photon  energy  interval 

within  the  bremsstrahlung  spectrum  k^:  J  -  1,  2,  ...»  i,  and  on  the 
endpoint  energy  of  the  spectrum  T., :  i  -  1,  2-m  ...»  n.  The  N*^  can  be 
calculated  from  functional  approximations  to  the  spectrum  N(T,k)  (see 
next  section) . 

In  viaw  of  the  above  "step-wise"  approximation,  the  reduced 
yield  function  becomes 


tCTj)  - 


1 

j«1 


Vr  1 


- 1, 


(18) 


where 


(18a) 


and 


r 


o (k)dk 


YSI 


(18b) 


Photon  Energy  (V^) 


Figure  3. 


T 


1 


K  "stepwise"  approximation  to  the  bremsstrahlung 


spectrum  of  endpoint  T^. 
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is  an  average  value  of  o(lc)  in  the  photon  energy  bin  between  -  AT  and 
kj.  Note  in  Eq.  (18)  that  the  yield  Y(T)  is  the  reduced  yield  Yr(T)  of 
Eq.  (16). 

In  matrix  notation,  Eq.  (18)  becomes 


Y  -  Na 


(19) 


''here  ?  is  a  column  matrix  with  n  rows  of  elements  Y^,  o  is  a  column 
matrix  with  n  rows  of  elements  o4,  and  N  is  an  nCh  order  triangular 
matrix  with  elements  N^,  N^  ■  0  when  j  >  i.  The  matrix  equation,  E^j. 
(19) ,  can  be  rewritten  as 


N,,  0 


N21  N22 


0 

0 


V  Nu2  «n3 


0 

0 


N 


nnj 


U°n 


(20) 


Three  methods  of  solution  to  the  above  aet  of  equations  have 
been  proposed.  The  first  method  is  the  "Photon  Difference  Method''  of 
Eats  and  Cameron  (35).  This  method  conslsta  of  solving  the  first  equa¬ 
tion  for  o^,  which  when  substituted  ir»  the  second  equation  gives  o^,  and 
so  on.  The  worst  defect  of  the  photon  difference  method  is  the  depen¬ 
dence  of  the  calculated  value  of  the  cross  section  at  a  given  energy  on 
the  results  of  calculations  at  all  preceding  energies.  The  errors 
arising  from  fluctuations  in  the  yield  curves  at  low  energy  are  thus 
cumulated  and  transmitted  to  all  cross  sections  at  higher  energies. 
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The  second  r*thod  is  the  "Matrix  Inversion  Method"  of  Penfold 
and  Lelss  (34) .  This  method  involves  an  inversion  of  the  N  matrix  to 
obtain  an  inverse  operator  N  \  which  when  applied  to  Eq.  (19)  yields 
the  solutions  o. 

Jr 

•  if*?  -  N^No  .  (21) 

The  matrix  N  *  consists  of  an  appropriate  set  of  what  Penfold  and  Lelss 
called  "B-numbers"  that  were  calculated  from  Schiff  spectra  approxima¬ 
tions.  The  matrix  product  N  represents  a  series  of  "weighting  func¬ 
tions"  that  relate  the  solution  N  to  the  actual  cross  section.  Penfold 
and  Lelss  have  calculated  sets  of  B -numbers  and  their  associated  weighting 
functions  for  thin-target  bremsstrahlung  of  endpoint  energies  in  the  range 
from  2  MeV  to  1  BeV  (19.34). 

-  -  The  Penfold  and  Lelss  method  has  been  found  to  be  extre&ely  sen¬ 

sitive  to  the  precision  with  which  yields  sre  measured.  Moderate  uncer¬ 
tainties  in  yield  measurement  may  appear  as  photonuclesr  resonances  often 
resulting  in  cross  sections  which  oscillate  excessively  as  a  function  of 
energy  and  sre  therefore  physically  unacceptable  (36,37).  A  third  method 
of  photonuclesr  yield  function  analysis  was  developed  by  Cook  (36)  which 
minimizes  the  undulation,  or  structure,  of  the  cross  section  function. 

The  method  has  been  called  the  "Least  Structure  Solution’1  of  photonuclesr 
yield  functions. 

The  least  structure  method  la  basically  a  systematic  smoothing 
technique  which  produces  the  most  uniform  set  of  cross  sections  consis¬ 
tent  with  the  experimental  data.  The  method  applies  all  smoothing  to  the 
cross  sections  o  and  not  the  yields  Y^,  and  it  assumes  no  functional 


form  of  the  o^.  Completely  numerical  methods  are  used.  The  amount  of 
smoothing  employed  is  determined  by  the  variability  of  the  input  data 
and  the  distortions  introduced  by  the  smoothing  can  be  estimated. 

Cook's  least  structure  solution  of  Eq.  (19)  selects  those  cal¬ 
culated  values  of  o^  such  that 


I  Nlj0j  =  1±  . 


In  other  words,  a  solution  o^  is  acceptable  if  the  calculated  yield, 

A 

Y^  «  Nij°j»  311  "acceptable  range"  of  the  measured  yield  Y^  at 

each  endpoint  energy  T^.  A  set  of  o^'s  are  considered  to  be  an  accep¬ 
table  solution  if 


0.  i  Ktvv1!  < 


„  (23) 


where  x  (o^)  is  a  function  closely  related  to  the  Chi-square  distribution 
of  statistics,  the  VY^  are  standard  deviations  of  the  measured  yield  at 
each  endpoint  energy,  and  n  is  the  number  of  data  point a  on  the  yield 


curve. 


The  actual  acceptable  solutions,  o^,  are  arrived  at  by  an  itera- 

-1 

tive  procedure.  The  procedure  consists  of  determining  the  matrix  M 


such  that  (see  Eq.  (21)] 


■  M  \o  °  Ro 


yields  the  acceptable  set  of  solutions  as  determined  by  Eq.  (23).  The 
matrix  M  ^  is  determined  by  calculus  of  variations  techniques  (38)  and 
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is  a  function  of:  the  bremsstrahlung  matrix  N;  the  inverse  of  the  brems- 
strahlung  matrix  N  a  smoothing  matrix  S  which  minimizes  the  difference 
between  successive  values  of  and  a  weighting  matrix  W  of  weighting 
factors  inversely  proportional  to  the  variance  of  the  yield  points.  The 
matrix  product  M  'Hi  has  been  called  the  resolution  function  R.  It  can  be 
thought  of  as  an  averaging  of  o  over  a  few  energy  bins  with  an  approxi¬ 
mate  Gaussian  function  weighting  (36). 

A  rigorous  explanation  of  the  Least  Structure  (LS)  method  for 

the  solution  of  photonuclear  yield  functions  is  beyond  the  scope  of  this 

work.  Such  an  explanation  is  available  in  references  (29,36,37).  How- 

27  25 

ever,  becaT  je  LS  was  used  to  determine  the  Al(y,2p)  Na  cross  section, 
which  was  investigated  in  this  work,  a  general  explanation  of  the  LS 
computational  procedure  is  provided.  The  explanation  is  basically  & 
description  of  a  Fortran  computer  program  which  incorporates  LS  and  which 
was  deviaed  by  Cook  (36) . 

Least  Structure  attempts  to  arrive  at  solutions  to  the  photo- 
nuclear  yield  function  [Eqs.  (18)  and  (19)] 


?1  -  l 


This  ia  accomplished  by  calculating  a  matrix  M  (to  be  defined  shortly) 


such  that 


fi '  £  Vj 


yields  solutions  that  do  not  oscillate  as  do  those  obtained  by  simply 
inverting  the  N  matrix  of  Eq.  (25),  The  matrix  H  is  of  the  fora 
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N 


ij 


vv 


(27) 


where  X  is  a  Lagrangian  multiplier,  ^  is  the  inverse  of  the  transpose 

of  the  (bremsstrahlung)  matrix,  is  the  inverse  of  the  weighting 

2 

matrix  •  1/(7Y^)  where  the  VY^  are  the  standard  deviations  of  the 
measured  yields,  and  S  is  a  smoothing  matrix  of  the  form 


n-1 

£2  (°3+l  •  2oJ  +  °3-l)2  • 


(28) 


The  matrix  S  (also  called  the  structure  function)  is  essentially  a 

measure  of  the  rate  of  change  of  the  slope  of  the  cross  section  as  a 

function  of  endpoint  energy.  The  LS  solution  to  Eq.  (25)  is  that  set  of 

o .  which  minimizes  the  structure  function  with  the  restriction  that  [see 
J 

also  Eq.  (23)] 

n 

x2  -  I  VrNij°j  "  V2  -  n  ‘  (29) 


In  the  LS  computer  program,  the  Lagrangian  multiplier  X  is  arbi¬ 
trarily  chosen  as  an  input  parameter.  The  equation 

?i  “  [N1;]  +  <30> 

is  then  solved  for  the  fixed  X.  The  solutions  o^  are  substituted  into 
2  2 

Eq.  (29)  and  x  is  computed.  The  calculated  x  is  compared  to  an  input 
2  2 

X  in  and  an  acceptable  error  Ax  If  the  solution  is  not  acceptable, 

a  new  X  is  selected  and  naw  solutions  are  obtained  until  both  Eqs.  (29) 
and  (30)  are  satisfied. 
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The  effectiveness  of  the  least  structure  solution  of  photonudear 

yield  functions  for  cross  section  determination  has  been  well  proven  and 

documented  in  the  literature  (37,39,40).  Cook  (37)  and  Anderson  (39) 

measured  the  cross  section  for  the  ^0  photoneutron  reaction,  a  reaction 

which  has  been  extensively  studied  both  experimentally  and  theoretically, 

using  the  least  structure  method  and  obtained  results  which  were  in  very 

good  agreement  with  existing  theories  and  with  results  obtained  by  other 

investigators.  The  soundness  of  the  technique  was  further  verified  by 

12 

Cook  (40)  when  he  measured  the  cross  section  for  the  C  photoneutron 
reaction  obtaining  results  consistent  with  those  previously  reported. 

The  least  structure  method  has  been  utilized  since  on  several  occasions 
for  photonudear  reaction  cross  section  determination  (41-43). 

The  Bremsstrahlung  Spectrum 

This  section  describes  the  bremsstrahlung  process  and  brems¬ 
strahlung  epectra.  The  theoretical  calculations  of  the  bremsstrahlung 
crose  section  are  discussed  in  generel  with  particular  emphasis  placed 
upon  celculstlons  performed  in  the  Born  approximation.  Finally,  the 
processes  that  leed  to  differences  in  bremsstrahlung  epectra  from  thin 
and  thick  tergets  are  described. 

Bremsstrahlung  is  electromagnetic  radiation  which  is  produced 
when  e  charged  particle  is  accelerated  in  the  field  of  e  charged  pertlcle. 
Bremset rah lung  as  used  here  describes  the  electromagnetic  radiation 
emitted  by  electrons  decelerated  upon  their  penetration  into  e  brems¬ 
strahlung  radietor  or  target.  The  bremsstrahlung  process  is  depicted 
schematically  in  Figure  4.  An  electron  of  charge  e  and  kinetic  energy 
T  comes  within  ths  Coulomb  field  of  a  nucleus  of  charge  +Ze  where  it  is 
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deflected  producing  a  photon  of  energy  It.  The  energy  k  of  the  photon  Is 
the  difference  between  the  electron's  initial  kinetic  energy  T  and  ita 
final  kinetic  energy  T'. 

k  -  T  -  T'  (31) 

Because  the  deflected  electron  can  have  any  one  of  an  essentially  Infi¬ 
nite  number  of  kinetic  energies  between  0  and  T,  the  photon  may  poetess 
any  energy  up  to  T. 

A  classical  electrodynamics  formulation  may  be  used  to  obtain 
an  approximate  functional  form  for  the  energy  spectrum  of  the  photon 

N(k)  -  Na  (32) 

where  N(k)  is  the  number  of  photons  of  energy  between  k  and  k  +  dk,  1^ 

2 

the  number  of  target  nuclei  per  cm  ,  and  do(k)/dk  Is  the  differential 

2 

crosa  section  for  bremsstrahlung  production  in  cm  .  The  bremsstrahlung 
cross  section,  describes  the  relative  probability  with  which  the  brems- 
atrahlung  process  will  occur.  A  simplified  expression  for  the  bremsatrah- 
lung  crosa  section  has  the  form  (44) 

-  Z2ro2a*(T.Z,k)/k  (33) 

where  Z  la  the  atomic  number  of  the  target  element ,  rQ  Is  the  claasical 
electron  radius,  o  is  the  fine  structure  constant  and  t(T,Z,k)  is  the 
so-called  Intensity  function.  From  Eqs.  (32)  and  (33)  it  Is  evident  that 
for  a  given  target  the  bremsstrahlung  spectrum  la  Inversely  proportional 
to  the  energy  of  the  emitted  photon.  Also  seen  Is  the  direct  relation- 
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ship  between  the  photon* s  energy  spectrum  and  the  intensity  function 
$(T,Z,*). 

An  accurate  expression  for  $(T,Z,k)  at  large  T  and  k  has  not 
been  obtained  classically.  To  describe  $(T,Z,k)  fully,  quantum  electro¬ 
dynamics  is  required.  A  most  general  quantvm-mechanical  analysis  pro¬ 
vides  the  following  approximate  expression  for  the  bremsstrahlung  cross 


section  (19) 


do  (k) 
dk 


(34) 


where  Qf  is  the  density  of  final  states,  pQ  is  the  initial  electron 
momentum,  c  is  the  speed  of  light  in  vacuo,  Eq  is  the  total  initial 
energy  of  the  electron,  and  is  the  matrix  element  for  the  transition 
from  the  Initial  state  i  to  the  final  state  f.  From  Eq.  (34)  it  is  seen 
that  the  bremsstrahlung  cross  section  depends  upon  the  square  of  the 
absolute  value  of  the  transition  matrix  element.  This  dependence  is 
Implicit  in  the  intensity  function  of  Eq.  (33).  The  transition  matrix 
element  in  turn  depends  upon  the  wave  functions  of  the  electron  initial 
and  final  states. 

To  obtain  an  exact  expression  for  the  bremsstrahlung  cross 
section  it  is  essential  to  know  the  exact  wave  functions  that  describe 
the  electron  in  the  Coulomb  field  of  the  nucleus.  However,  the  Dirac 
wave  equation  for  an  electron  in  a  Coulomb  field  does  not  have  finite 
solutions  because  the  wave  function  1c  ^his  case  yields  an  Infinite 
series  (19,45,46).  Thus,  approximate  wave  functions  must  be  used  to 
calculate  the  transition  matrix  elements.  In  addition,  some  simplifying 
assumptions  must  be  made  in  the  calculations  themselves.  Therefore, 
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only  approximations  to  the  bremsstrahlung  cross  sections  are  obtainable. 

Basically ,  the  approximations  made  in  the  calculation  o£  the 
bremsstrahlung  cross  section  fall  into  one  of  two  categories:  1)  cal¬ 
culations  performed  in  the  Bora  approximation ,  and  2)  extreme  relativis¬ 
tic  calculations.  In  the  Born  approximation  calculations,  the  Coulomb 
field  is  considered  a  first-order  perturbation  of  the  electron  initial 
state  *  The  electron  initial  and  final  wave  functions  are  described  aa 
plane  waves  (51,52).  Extreme  relativistic  calculations  (45,49)  essen¬ 
tially  encompass  the  calculations  that  do  not  use  the  Born  approximation. 
In  extreme  relativistic  calculations,  Sommer feld-Maue  wsve  functions  in 
a  screened  (atomic  electrons  "screen"  the  nuclear  charge)  Coulomb  poten¬ 
tial  are  used.  These  wave  functions  are  described  as  plane  and  spherical 
waves  (45).  It  turns  out  that  extreme  relativistic  calculations  are 
valid  only  for  initial  electron  kinetic  energies  exceeding  about  50  MeV 
(19).  For  this  reason,  and  due  to  the  fact  that  energies  of  interest 
here  are  well  below  that  value,  extreme  relativistic  calculations  are 
not  considered  further  In  this  treatment;  Bora  approximation  calculations 
are,  hence,  emphasized. 

The  Born  condition  restricts  the  atomic  number  of  the  target 
element  and  the  electron  initial  and  fiual  velocities  such  thst 

2ttZ/137P  «  1  (35) 

where  $  -  v/c  where  v  is  the  electron's  velocity  and  c  is  the  speed  of 
light.  It  is  evident  that  the  conditions  of  the  Bora  approximation  are 
not  strictly  mit  for  high  Z  tsrgets  such  as  those  made  of  tungsten  (Z  ■ 
74).  In  addition,  the  theory  falls  st  the  high-energy  end  of  the  spectrum 
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(the  so-celled  tip)  when  the  energy  of  the  emitted  photon  approaches  the 
kinetic  energy  of  the  incident  electron  and  the  velocity  of  the  recoil 
electron  is  small  compered  to  the  speed  of  light.  However,  the  use  of 
the  Born  approximation  formulas  predict  bremsstrahlung  spectre  which 
agree  reasonably  well  with  experiment  (19,46). 

The  most  important  contributions  to  bremsstrahlung  cross  section 
formulations  obtained  within  the  framework  of  the  Bora  approximation  ere 
due  to  Bethe  and  Heitler  (47),  Heitler  (48)  and  Schiff  (20).  Bethe  and 
Heitler  obtained  an  expression  for  the  bremsstrahlung  cross  section  which 
is  differentiel  in  electron  angle  and  in  photon  angle  and  energy.  Their 
expression  was  derived  on  the  assumption  that  the  field  of  the  nucleus 
is  a  pure  Coulomb  field,  i.e.,  no  nuclear  screening  by  atomic  electrons 
was  considered.  Nudeer  screening  is  an  important  parameter  which  must 
be  taken  into  consideration,  particularly  in  the  case  of  bremsstrahlung 
spectra  from  high  Z  targets.  Schiff  integrated  the  Bethe-Heitler  differ¬ 
entiel  cross  section  assuming  complete  screening  (20) .  The  expressions 
he  obtained  apparently  show  the  best  agreement  with  available  experimental 
data  (19,46,50).  The  Schiff  expressions  have  been  widely  used  in  the 
analysis  of  experimental  results  obtained  with  bremsst rah lung  spectra 
from  high-energy  accelerators. 

In  his  calculations,  Schiff  started  with  the  Bethe-Heitler  cross 
section  and  integrated  it  assuming  that:  1)  the  initial  and  final  ener¬ 
gies  of  the  electron  are  large  compared  to  its  rest  energy;  2)  the  Cou- 

/a 

load)  field  of  the  atom  may  be  represented  by  the  potential  (Ze/r)e 
where  as^Z1^  (Thcmas-Fermi  model);  3)  terms  of  the  order  of  (Z*^/c)2 
may  be  neglected.  Under  these  assumptions  he  obtained  the  following 
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expression  for  the  bremsstrahlung  cross  section  differential  In  photon 
emission  angle  and  energy 


and 

x  ■  E  Wu  .  (36b) 

o 

II 

In  Eqa.  (36),  Eq  is  the  Initial  total  electron  energy,  E  Is  the  final 
total  electron  energy,  k  (-  Eq  -  E)  la  the  energy  of  the  radiated  photon, 

Z  la  the  atomic  number  of  the  target  element,  rQ  la  the  classical  elec¬ 
tron  radius ,  u  la  the  rett  energy  of  the  electron,  and  x  Is  ths  reduced 
angle  of  emission  of  the  radiated  photon  which  la  a  function  of  the  angle 
e,  In  radians,  at  which  the  photon  la  radiated  with  respect  to  the  Initial 
electron  direction.  The  use  of  the  atomic  Thomas-Fermi  model  in  consid¬ 
ering  the  affects  of  screening  lesd  to  exaggerated  vslues  of  do(k,x),  but 
the  error  Introduced  does  not  exceed  4Z  In  the  worst  cass  when  Z  la  large 
and  screening  is  intermediate  (19,20,46). 

Thick-Target  Bremsstrahlung  Spectra 
The  arguments  previously  presented  considered  only  bremsstrahlung 
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epectra  from  "chin"  targats.  A  chin  CargeC  is  a  bremsscrahlung  radiator 
of  thickness  such  that  only  one  electron  interaction  occurs.  Tha  shape 
of  the  thin-target  spectrum  reflects  the  transition  probabilities  between 
the  electron  initial  and  final  states.  Most  all  linear  accalerators 
possass  "thick"  targets  becnusa  of  the  substantially  larger  photon  flu- 
ences  attainable.  A  thick  target  is  a  radiator  of  thickness  such  that 
many  electron  interactions  occur.  Such  interactions  Include  colllsional 
as  well  as  radiatlonal  interactions .  The  radiated  photons  can  also  inter¬ 
act  within  the  target.  Thus,  in  the  thick-target  situation,  tha  shape  of 
the  bremsst  rah  lung  spectrum  additionally  reflects  the  daformation  of  the 
thin-target  spectrum  due  to  tha  Interactions  in  tha  target  of  tha  alec- 
trons  and  photons.  Tha  difference  in  spectral  shape  between  bremsstrah- 
lung  produced  in  thin  targats  in  contrast  to  that  produced  in  thick  ter- 
gets  is  schematically  represented  in  Figure  5. 

Thick-target  bremsstrahlung  spectra  calculations  in  the  enargy 
ranga  from  about  1  MeV  to  about  50  MeV  heva  been  performed  by  essentially 
two  methods:  1)  by  Monte  Carlo  calculations  such  as  thoaa  of  Berger  and 
Seltzer  (51),  and  2)  by  numerical  analysis  methods  such  as  those  employed 
by  Lent  and  Dickinson  (52)  and  by  Ferdlnande  end  co-workers  (53) .  In  the 
present  work,  a  numerical  analysis  method  similar  in  many  reepects  to 
thet  of  Ferdinands  was  usad  to  eetimate  the  bremsscrahlung  epectra  of  tha 
linear  eccalerator  which  wes  used.  The  approech  is  based  on  Hisdal's 
(54)  method  of  correcting  intermediate  target  spactrs  for  the  effects  dua 
to  electron  multiple  scattering.  Tha  Molldre  (55)  multiple  scetterlng 
distribution  is  aseumed.  The  method  is  described  in  detail  in  tha  Brens- 
strahlung  Spectrum  section  of  Chapter  III  -  Materials  and  Methods. 


Energy  Fluence 


Photon  Energy 


Figure  5.  Breuss trahlung  spectre  from  thin  and 
thick  targets.  Curves  have  been  normalised  to 
equal  areas. 


Measurement  of  Photonuclear  Yield 


The  solution  of  the  photonuclear  yield  function  requires  that 

the  reduced  yield  Yr(T),  photonuclear  reactions  at  a  certain  bremsstrah- 

2 

lung  endpoint  energy  T  per  target  atom  per  cm  ,  be  determined.  This  Is 
accomplished  by  measuring  the  normalized  yield  Yr(T),  the  number  of  pho¬ 
tonuclear  reactions  per  unit  monitor  response,  multiplying  the  normalized 

yield  by  the  monitor  response  function  of  Eq.  (13),  F(T) ,  and  dividing 

2 

this  product  by  n£,  the  number  of  target  atoms  per  cm  In  the  sample. 

The  normalized  yield  Yq(T)  is  obtained  experimentally.  If  photonuclear 

reactions  produced  in  a  sample  result  In  radioactive  nuclei,  then  the 

number  of  such  reactions  is  simply  a  function  of  the  radioactivity  induced 

In  the  sample.  If  a  radiation  monitor,  such  as  an  Ionization  chamber.  Is 

Irradiated  simultaneously  with  the  sample  such  that  Its  response  is 

directly  proportional  to  the  bremsatrahlung  beam  energy  flux,  then  the 

normalized  yield  Yq(T)  can  be  obtained  from  a  measurement  of  the  sample 

activity  and  from  the  monitor  response.  This  section  describes  the 

methodology  by  which  the  normalized  yield  Is  determined. 

The  normalized  yield  Y  (T)  Is  the  quotient  of  two  quantities 

n 

which  are  each  determined  experimentally. 

Yn(T)  -  No(T)/Vq(T)  (37) 

In  Eq.  (37),  Nq(T)  is  the  total  number  of  photonuclear  reactions  Induced 
in  the  sample  at  a  given  bremsatrahlung  endpoint  energy,  and  Vq(T)  is  a 
radiation  monitor  response  which  Is  proportional  to  the  bremsatrahlung 
beam  energy  flux  at  a  given  endpoint  energy.  Nq(T)  Is  determined  from 
a  measurement  of  the  sample  activity.  Vq(T)  Is  determined  from  the 
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response  of  a  monitor  system  which  will  be  described  shortly. 

The  total  number  Nq  of  photouuclear  reactions  Induced  may  be 
expressed  as  the  product  of  the  photonuclear  reaction  production  rate  P^ 
and  the  Irradiation  time  t 

r 


N 

o 


P„t 
N  r 


(38) 


If  the  reactions  result  In  radioactive  nuclei  of  decay  constant  X,  then 
the  number  N(t)  of  radioactive  nuclei  present  In  the  sample  at  any  given 
time  t  Is  related  to  the  constant  production  rate  P^  through  the  expres¬ 
sion 

N(t)  -  j  PN[1  -  exp(-Xt)]  .  (39) 

The  relationship  of  Eq.  (39)  Is  shown  in  Figure  6.  The  number  of  radio¬ 
active  nuclei  ’'grows'*  with  a  1-exp  (-Xt)  dependence  with  PN/X  as  an  upper 
limit.  Solving  Eq.  (39)  for  P^  and  substituting  the  resulting  expression 
in  Eq.  (38)  results  in  the  following  equation  for  Nq(T)  in  terms  of  N(t) 

Nq(T)  -  N(t)Xtr/U  -  exp(-Xt)  ]  .  (40) 

The  number  of  radioactive  nuclei  at  any  given  time,  N(t),  is 
rather  sensitive  to  fluctuations  In  the  production  rate,  which  may  be 
caused  by  variations  In  bremsstrahlung  beam  Intensity.  For  this  reason, 
it  Is  uesirable  that  the  radiation  monitor  respond  In  time  to  the  beam 
Intensity  in  the  same  fashion  as  does  N(t).  To  this  purpose,  the  output 
of  the  bremsstrahlung  monitor,  an  NBS  type  P2  ionization  chamber  (56), 
is  delivered  to  an  RC  electric  circuit  as  shown  in  Figure  7  (57).  The 
response  V(t)  of  a  circuit  of  this  type  as  a  function  of  time  is  shown 


Figure  7.  RC  circuit  of  radiation  monitor  system. 
HV  -  High  Voltage.  IC  -  Ionization  Chamber. 

R  -  Resistor.  C  -  Capacitor.  E  -  Electrometer. 
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In  Figure  8.  Since  both  the  sample  and  the  monitor  are  exposed  sinultan- 
eously  to  the  same  bremsstrahlung  beam  Intensity,  V(t)  varies  with  time 
as  does  N(t) . 

The  total  monitor  response  Vq(T)  is  the  product  of  the  change  in 
potential  as  seen  by  the  electrometer  and  of  the  time  the  potential 
undergoes  change,  the  irradiation  time  t^ 

V  (T)  -  P  t  .  (41) 

o  v  r 

The  rate  of  change  of  the  potential  across  the  RC  circuit  of  Figure  7  can 
be  described  by  the  differential  equation 

srn. .  P  .  m  (42) 

at  v  rc  '  1 

where  V(t),  R,  and  C  are  the  circuit's  voltage  as  a  function  of  jtime  and 
the  resistance  and  capacitance,  respectively.  Integration  of  Eq.  (42) 
over  time  yields  the  following  expression  for  the  voltage  V(t) 

V(t)  -  RCPv£l  “  exp(-t/RC) ]  .  (43) 

Solving  Eq.  (43)  for  P^  and  substituting  in  Eq.  (41)  results  in  the 
following  expression  for  Vq 

Vo  *  V(t)tr/RC[l  “  exp(-t/RC)J.  (44) 

The  normalized  yield  of  Eq.  (37)  is  now  specified  through  Eqs. 
(40)  and  (44).  If  the  RC  constant  of  the  circuit  is  adjustable  by  means 
of  a  variable  resistance  such  that 

(45) 


RC  -  1/A 


Figure  8.  Growth  of  potential  across  the  capacitor  C 
of  the  RC  circuit  of  Fig.  7. 


then  Eq.  (37)  reduces  to 


Yq(T)  =  N(tr)/V(tr)  (46) 

where  N(t^)  is  the  number  of  radioactive  nuclei  In  the  sample  Immediately 
following  an  Irradiation  of  duration  t^,  and  V(tr)  Is  the  integrated  mon¬ 
itor  system  voltage  measured  by  the  electrometer  after  the  ionization 
chamber  has  been  irradiated  for  a  time  period  tf. 

A  sample  containing  N(t^)  radioactive  nuclei  of  decay  constant 
X  decays  with  an  activity 

A(tr)  -  N(tr)X  (A 7) 

disintegrations  per  unit  time.  The  sample  activity  is  decreased  during 
the  time  the  sample  is  transported  from  its  irradiation  location  to  its 
appropriate  position  in  the  radioactivity  counter.  If  this  transit  time 
is  denoted  t^,  the  activity  of  the  sample  at  the  start  of  the  counting 
period  may  be  given  by 

A(tt)  -  A(tr)exp(-Xtt)  .  (A8) 

The  sample  activity  at  time  t^  from  Eq.  (A 7)  is  substituted  in  Eq.  (A8) 
to  yield  an  expression  for  the  sample  activity  at  time  t^  as  a  function 
of  the  number  of  radioactive  atoms  in  the  sample  at  time  tr 

A(tt)  »  N(tr)Xexp(-Xtt) .  (A9) 

The  sample  activity  is  now  measured  using  a  gamma- ray  pulse- 
height  spectrometer  system.  Let  C  be  the  number  of  counts  registered  in 
a  pulse-height  interval  (or  window)  &E  centered  at  the  photopeak  of  an 
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appropriate  gamma-ray  energy  E^  of  the  radioisotope  assayed.  C  is  related 
to  the  activity  of  the  sample  at  the  start  of  the  count  A(t^)  and  to  the 
counting  time  t  through  the  integral  equation 

rtc 

C  -  K  /  A(tt)[exp(-A:.^)]dt  (50) 

The  term  in  brackets  accounts  for  radioactive  decay  during  the  counting 
period.  The  constant  of  proportionality  K  (counts  per  disintegration) 
relates  the  total  number  of  sample  disintegrations  integrated  over  time 
tc  to  the  net  number  of  counts  registered  in  the  pulse-height  window  of 
the  counting  system. 

Integration  of  Eq.  (50)  produces 

C  "  -  A(tt)[l  -  exP(-Atc)]  .  (51) 

If  the  expression  for  A(tt)  of  Eq.  (49)  is  now  substituted  in  Eq.  (51) 
and  the  resulting  equation  is  solved  for  N(tf),  the  number  of  radioactive 
nuclei  in  the  sample  at  time  t^,  one  obtains 

N(tJ  -  C/Kt  (52) 

where 

t  -  exp(-Xtt)U  -  exp(-Xtc)]  .  (52a) 

The  normalized  yield  Y  (T)  of  Eq.  (46)  is  now  completely  specified.  If 

n 

we  now  rename  the  quantity  V(t  )  of  Eq.  (46)  as  simply  V  and  substitute 
Eq.  (52)  into  Eq.  (46),  the  normalized  yield  becomes 
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Yn(T)  -  C/KtV.  (53) 

This  Is  thu  required  normalized  yield  axpresslon  [Eq.  (37)]  as 
a  function  of  experimentally  determined  quantities.  Yn(T)  Is  the  nor¬ 
malized  yield.  In  number  of  photonuclear  reactions  per  unit  monitor 
response,  required  for  the  solution  of  the  photonuclear  yield  function. 

C  Is  the  net  number  of  counts  in  the  gamma-ray  spectrometer  pulse-height 
-  window.  The  proportionality  constant  K  relates  the  total  number  of  radlo- 
actlva  decay  events  of  the  sample  to  the  racioactivity  count  C.  The  con¬ 
stant  K  is  a  function  of  the  counter  efficiency  for  an  emitter  of  the 
specific  gamsa-ray  energy  counted,  of  the  relative  intensity  of  the 
ganana  ray,  and  of  sample  self -absorption.  V  Is  the  output  voltage  signal 
of  the  RC  circuit  of  the  radiation  monitor  system  Integrated  ovar  the 
Irradiation  time  tr.  And  x,  defined  by  Eq.  (52a),  Is  simply  a  function 
of  tha  radioisotope  decay  constant  X  and  of  the  sample  transit  and 
counting  times  t  and  t  ,  respectively. 

Photonuclear  Ratios 

In  this  section  are  derlvad  the  formulae  used  to  calculate  pho- 
tonudaar  ratios.  Threa  assumptions  have  been  made  in  the  derivation. 
First,  photonuclear  reactions  result  In  gamma-emitting  radioactive  nuclei 
such  that  the  photonuclear  yield  may  be  determined  by  a  gan&a-ray  detec¬ 
tion  system.  Second,  the  elements  in  which  photonuclear  reactions  are 
produced  are  contained  in  a  single  chemical  compound.  Third,  the  radio¬ 
activity  produced  in  each  element  can  be  datermlned  separately.  The 
formulae  are  derived  in  a  manner  analogous  to  the  derivation  of  the  equa¬ 
tion  for  the  calculation  of  normalized  yield. 
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The  number  of  photonuclear  reactions  induced  in  a  sample  irradi¬ 
ated  at  constant  intensity  for  a  time  t  may  be  given  as 

N'  -  PNtr  (54) 

where  N*  is  the  total  number  of  photonuclear  reactions  produced  in  the 
aamplet  is  the  constant  photonuclear  reaction  production  rate  (reac¬ 
tions  per  unit  time),  and  t  is  the  irradiation,  time.  If,  as  commonly 
occurs,  the  photonuclear  reactions  result  in  radioactive  nuclei,  these 
decay  as  they  are  produced.  If  the  radioactive  species  decay  with  decay 
constant  X,  .then  the  sample's  activity  immediately  following  irradiation 

Atr  18 

Atr  ‘  V1  -  e"Ul>  (55) 

Solving  Eq.  (55)  for  P^  and  then  substituting  the  resulting  expression 
in  Eq.  (54)  yields 

N*  -  Atrtr/(1  -  e"Xtr)  (56) 

If  the  aample  radioactivity  is  counted  at  a  location  different  from  ita 
irradiation  position,  the  sample  decays  during  the  time  it  ia  transported 
from  its  location  in  the  irradiation  room  to  the  radioactivity  counter. 

If  this  transit  time  is  denoted  t£,  then  the  activity  of  the  sample 
after  transit  and  at  the  onset  of  counting  is 

Att  "  Atr(e"Xtt)  (57) 

where  is  the  activity  of  the  sample  which  remains  after  it  has  been 
removed  from  its  location  in  the  Irradiation  room  and  has  been  placed  in 
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the  counter.  Solving  Eq.  (57)  for  Atr  and  substituting  the  resultant 
expression  in  Eq.  (56)  yields  the  following  expression  for  N': 

N*  “  Atttr/(1  -  e‘Xtr)(e"Ut)  .  (58) 


The 

dow 


C. 


at 


At  this  point,  the  sample  is  counted  in  a  gamma-ray  spectrometer, 
total  number  of  counts  registered  in  a  pulse-height  interval,  or  win- 
AE^,  centered  about  the  photopeak  of  an  appropriate  gamma-ray  energy 
is  integrated.  The  sample  is  counted  for  a  time  t^  and  the  net  number 
counts  (after  background  subtraction)  registered  in  the  photopeak  is 
The  number  of  net  counts  C  is  related  to  the  sample's  activity  A 
the  stare  of  the  count  through  the  integral  equation 


C  -  K 


(59) 


where  the  exponential  within  the  Integral  accounts  for  sample  decay  during 
the  counting  period.  The  constant  of  proportionality  K  relates  the  dis¬ 
integrations  of  the  aample  to  the  counta  registered  in  the  spectrometer 
pulse-height  window.  It  is  a  function  of  the  number  of  gamma  rays  per 
decay  n,  the  counting  system  photopeak  efficiency  £  for  gamma  rays  of 
energy  E  from  the  sample  if  there  were  no  self-absorption,  and  the  sam¬ 
ple  aelf-absorption  factor  a 


K  -  n  £a  .  (60) 

Integration  of  Eq.  (59)  yields 

c '  r  \t(1  -  e'Uc>  • 


(61) 


Solving  Eq.  (61)  for  A  and  substituting  in  Eq.  (58)  yields  an  expression 
for  N'  in  terms  of  the  measured  counts  C 


In  Eq.  (63a),  is  Avogadro's  number,  A  is  the  atomic  mass  of  the  target 

nucleus,  m  is  the  mass  of  the  sample,  y  is  the  fractional  natural  abun- 
8 

dance  of  the  target  Isotope,  and  p  is  the  fraction  of  the  sample  molecular 
weight  made  up  by  the  target  element.  Tne  fraction  p  is 

p  "  xA/M^  (63b) 

where  A  is  the  atomic  weight  of  the  target  element,  is  the  molecular 
weight  of  the  compound,  and  x  is  the  number  of  atoms  of  the  target  ele¬ 
ment  per  molecule  of  compound. 

The  final  expression  for  the  number  of  photonuclear  reactions 


per  target  nucleus  induced  in  the  sample  is 
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CXt  A 

jl  ■  — 

n  €«naVyt 

where 

T  -  (1  -  e_Uc)(l  -  e_Xtr)(e’Xtt)  . 

N  Is  the  number  of  photonuclear  reactions  produced  In  the  sample  per  tar¬ 
get  nucleus,  C  Is  the  number  of  net  photopeak  counts  integrated  over  the 
count  time,  \  is  the  decay  constant  of  the  radlalaotope  produced,  A  is 
the  atomic  mass  of  the  target  element,  n  is  the  number  of  gamma  rays  of 
Interest  emitted  by  the  radioisotope  produced  per  disintegration,  £  is 
the  count lng-ays tern  efficiency  without  self-absorption  for  the  particular 
gamma  ray  in  photopeak  counts  per  disintegration,  a  is  the  self -absorption 

factor  of  the  sample  for  the  gamma  ray,  N  is  Avogadro's  number,  c  is 

A  8 

the  sample  mass,  p  is  the  fraction  of  the  sample  molecular  weight  made 
up  by  the  target  element  [as  given  by  Eq.  (63b)  J,  y  is  the  natural  abun¬ 
dance  of  the  target  isotope,  and  t  ,  t  and  t  are  the  counting,  irradia- 

c  r  t 

tion  and  transit  times,  respectively. 

Equations  (64)  are  similar  to  Eq.  (1)  in  the  paper  of  Nath  and 
Schulz  on  photoactivation  ratios  (9) .  Their  photoactivation  ratios  are 
expressed  in  terms  of  activity  Am  that  would  be  produced  in  the  sample 
by  photonuclear  reactions  if  it  were  Irradiated  indefinitely.  For  irradi¬ 
ation  time a  much  longer  than  the  half-life  of  the  radioisotope  being 
produced,  the  activity  of  the  sample  following  irradiation  [as  given  by 
Eq.  (55)]  becomes 

*-  ■  V1  *  «"Xtr>  •  PK 


(64) 


(64a) 


(65) 
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because  e  r  approaches  zero  when  tf  »  half-life  of  the  radioisotope 
produced.  Substituting  PN  -  in  Eq.  (54)  yields 

N1  -  As#tr  (66) 

where  is  the  activity  that  would  be  induced  in  the  sample  following 

an  indefinitely  long  irradiation  (equilibrium  activity) .  The  value  of 

N'  resulting  from  Eq.  (66)  is  now  substituted  in  Eq.  (62)  and  the  activity 

A  is  divided  by  the  number  of  target  nuclei  n  to  obtain 

0 


where  R  is  the  activity  (in  disintegrations  per  second)  per  target  nucleus 
Induced  in  the  sample  following  an  indefinitely  long  irradiation.  The 
remaining  variables  have  been  previously  defined.  Equation  (67)  is  iden¬ 
tical  to  Nath  and  Schulz's  Eq.  (1)  if  it  is  noted  that  for  very  thin 
foils  of  pure  elements,  which  they  used  as  samples  in  their  experiment, 
the  sample  self-absorption  factor  a  and  the  fraction  p  of  the  sample 
molecular  weight  made  up  by  the  target  element  equal  1.0. 

The  samples  used  in  this  study  consisted  of  chemical  compounds 
containing  the  two  elements  in  which  the  number  of  photonuclear  reactions 
induced  were  of  interest.  The  parameter  used  to  characterize  bremsstrah- 
lung  endpoint  and  spectral  quality  is  the  ratio  of  the  number  of  photo- 
nuclear  reactions  per  target  nucleus  induced  in  element  1  to  the  number 
of  photonuclear  reactions  per  target  nucleus  in  element  2.  The  photo- 
nuclear  ratio  PR  is  defined  as 


PR  - 


(68) 
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where  and  are  the  number  of  photonuclear  reactions  produced  per 
target  nucleus  in  elements  1  and  2,  respectively.  The  photonuclear  ratio 
PR,  in  terms  of  the  numbers  of  counts  and  from  the  radioisotopes 
produced  in  elements  1  and  2,  is  given  by 


PR  - 


C1*1A1/,T11  ^  lalplYlTl 


(69) 


Equation  (69)  follows  from  the  quotient  of  N^/N^  as  B*ven  by  Eqs.  (64, 

64a).  The  terms  t  ,  N  ,  and  m  are  constant  and  drop  out  of  the  ratio, 
r  a  s 

Equation  (69)  is  of  a  very  general  nature  and  applies  to  photonuclear 
ratio  calculations  in  almost  any  experimental  situation. 

Equation  (69)  can  be  simplified  somewhat  if,  as  occurs  in  this 
experimental  work,  the  energias  E^  of  the  gamma  reys  producing  the  counts 
and  are  equal.  Under  this  circumstance,  the  det-ctor  afficiency  £ 
and  the  sample  salf-ebsorption  factor  a,  both  of  which  are  a  function  of 
photon  anergy,  drop  out  of  the  retio.  The  photonuclear  ratios  of  this 
experiment  were  computed  using  the  simplified  formula 


PR 


VlVVlVl 

C2X2A2/n2’12Y2T2 


(70) 


Note  that  the  photonuclear  ratio  given  by  Eq.  (69)  is  algebraically 
equivelent  to  the  retio  R^/R^  8*ven  by  Eq.  (67). 


CHAPTER  III 


MATERIALS  AND  METHODS 

Materials 

Samples 

Samples  used  for  photonudear  ra.  -.u  determination  utilizing  the 
method  of  induced  activity  should  satisfy  certain  criteria: 

1.  Sample  elements  in  which  photonudear  reactions  are  to  be 
studied  must  contain  nuclides  which  produce  radioactive  isotopas  upon 
irradiation.  Tha  radioactive  specimen  should  be  a  gamma-emitter  of  large 
enough  branching  ratio  that  an  activity  measurement  can  be  appropriately 
made  by  counting  that  particular  gamma  ray. 

2.  Tha  half-life  of  the  resultant  radioisotope  should  ba  short 
enough  to  allow  a  fair  amount  of  radionuclide  production  upon  irradiation 
for  a  reasonabla  amount  of  time.  The  radioisotope  half-life,  on  the  other 
hand,  should  be  long  enough  so  that  a  sizable  fraction  of  the  radioacti¬ 
vity  remains  aftar  the  sample  has  been  transferred  from  the  irradiation 
room  to  the  counting  system. 

3.  The  specific  target  isotope  to  be  activated  should  have  a 
high  natural  abundance. 

4.  The  sample  material  should  be  stable  and  preferably  in  solid 
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form  at  room  temperature.  The  sample  material  should  be  readily  avail¬ 
able  and  Inexpensive. 

In  addition  to  the  above  criteria,  the  samples  used  in  this  study 
met  the  following  conditions: 

1.  The  sample  consisted  of  a  chemical  compound  containing  the 
two  elements  for  which  the  photonuclear  yields  were  of  Interest.  Because 
of  this,  both  elements  In  the  compound  were  exposed  to  exactly  the  same 
bremsstrahlung  beam  fluence,  and  measurement  of  the  fluence  to  which  the 
sample  is  exposed  was  unnecessary. 

2.  The  radioisotopes  produced  by  photonuclear  reactions  In  each 
element  emitted  the  same  energy  gamma  ray.  This  condition  made  determina¬ 
tion  of  the  counting-system  efficiency  as  a  function  of  gamma-ray  energy 
unnecessary.  The  condition  further  restricted  the  half-lives  of  the 
radioisotopes  produced.  The  half-lives  of  the  radioisotopes  had  to  be 
sufficiently  different  so  that  the  activity  of  each  could  be  determined 

by  counting  the  sample  at  two  different  times.  The  procedure  for  doing 
this  is  explained  in  detail  in  the  Methods  section  of  this  chapter. 

3.  The  elements  in  the  compound  had  to  possess  photonuclear 
cross  sections  differing  in  threshold  energy  and  giant  resonance  location 
such  that  the  photonuclear  yield  taken  as  a  function  of  bremsstrahlung 
endpoint  energy  in  one  element  remained  fairly  constant  while  the  yield 
in  the  second  element  changed  rapidly  in  the  same  energy  range. 

The  primary  aim  of  determining  photonuclear  ratios  was  to  obtain 
an  index  of  spectral  quality  for  bremsstrahlung  in  the  megavoltage  region 
and  particularly  in  the  energy  range  from  15  to  30  MeV,  where  such  ratios 
are  lacking.  No  one  sample  was  found  that  was  suitable  for  production  of 
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unique  photonuclear  ratios  as  a  function  of  bremsstrahlung  endpoint  energy 
encompassing  the  entire  15  to  30  MeV  range.  Therefore,  the  energy  range 
was  divided  and  two  sets  of  ratios  were  determined  requiring  two  differ¬ 
ent  samples.  In  the  lower  half  of  this  energy  range  (from  approximately 
15  to  about  22  MeV) ,  the  sample  of  choice  was  potassium  hexafluorosilicate 

(^SiFg) .  The  photonuclear  reactions  producing  the  ratios  were  the 
39  38  19  18 

K(y,n)  K  and  the  F(y,n)  F  reactions.  In  the  upper  half  of  the 
range  (approximately  23  to  30  MeV),  the  12C(y,n)11C  and  1®F(y,n)18F 
reactions  in  polytetrafluoroethylene  (FjOCFj) ,  or  "teflon,"  produced  the 
desired  ratios.  The  relevant  characteristics  of  these  reactions  are  sum¬ 
marized  in  Table  1. 

The  KjSiFg  samples  consisted  of  approximately  30.7  g  of  99X 
purity  I^SiFg  powder  contained  in  a  cylindrical  aluminum  can  3.2  cm  in 

•*r 

diameter  and  4.1  cm  high.  The  activity  induced  in  the  aluminum  was 

extremely  short-lived  (6.34  sec)  compared  to  the  half-lives  of  the  radio- 

38  18 

isotopes  of  interest  (  K  -  7.63  min,  F  -  109.8  min).  The  aluminum 
activity  was  allowed  to  decay  before  the  sample  was  counted.  The  mole¬ 
cular  weight  of  I^SiFg  is  220.25.  The  potassium  (atomic  weight  39.098) 
fraction  of  the  sample  molecular  weight  is  0.355.  The  fluorine  (atomic 

weight  18.998)  fraction  of  the  sample  molecular  weight  is  0.517.  The 

39  19 

natural  abundance  of  K  is  93. IX,  that  of  F  is  100Z. 

The  teflon  samples  consisted  of  2.5-ca-long  pieces  of  teflon  rod 

1.3  cm  in  diameter.  The  weight  of  each  sample  was  approximately  7.0  g. 

The  molecular  weight  of  teflon  is  100.02.  The  carbon  (atomic  weight 

12.011)  fraction  of  the  sample  molecular  weight  is  0.24.  The  fluorine 

12 

fraction  is  0.76.  The  natural  abundance  of  C  is  98.89Z. 


Table  1 


Photonuclear  Reactions  Used  for  Ratios 


r-  .  - 1 

-  - 1 

12c<  r,n)nc 

19F(  f,n)l8F 

39K(  y,n)38 K 

Threshold 

(MeV) 

18.7 

10.4 

13.1 

Half-life  of 

Product  (min) 

20.3 

109.8 

7.63 

*.  t 

Decay  const. 

(min  1) 

.034 

.0063 

.091 

Gamma  Energy 

(MeV) 

.511 

.511 

.511 

Photons 

per  decay 

1.98 

1.94 

1.98 
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The  aluminum  samples  used  to  measure  the  Al(y,2p)  Na  yield 

curve  consisted  of  aluminum  shot  ('v  1.0  mm  diameter)  tightly  encased  In 

aluminum  cans  Identical  to  those  used  for  the  K^SIF,  samples.  The  aluml- 

2  o 

27 

num  was  99.95%  pure.  The  natural  abundance  of  the  A1  Isotope  is  100%. 

The  sample  mass  averaged  36.5  grams  with  a  maximum  deviation  from  the 

27  25 

mean  of  less  than  0.1%.  The  characteristics  of  the  Al(y,2p)  Na  reac¬ 
tion  are  summarized  in  Table  2. 

Sample  Holders 

The  precise  determination  of  yields  required  that  samples  be 
irradiated  and  counted  in  a  very  reproducible  fashion.  To  satisfy  this 
requirement  purpose,  sample  holders  were  designed  and  constructed  (57). 

The  sample  holder  used  for  sample  irradiation  allowed  rapid  Insertion  and 
removal.  It  hald  the  sample  close  to  the  bremsstrahlung  target, ..thus 
Increasing  the  photon  beam  fluence  accordingly.  Furthermore,  it  permit¬ 
ted  very  reproducible  sample  positioning.  The  sample  holder  w»ed  for 
sample  counting  also  permitted  rapid  sample  insertion  and  removal.  It 
allowed  reproducible  positioning  of  the  samples  at  any  one  of  several 
precisely  measured  distances  from  tha  detector.  These  holders  are 
descrlbad  in  the  following  paragraphs. 

The  sample  holder  assambly  used  for  irradiation  is  shown  in 
Figura  9a.  It  contained  a  radial  collimator  constructed  of  poured  cerro- 
bend  alloy  (Lipowitz  alloy,  attenuation  properties  similar  to  lead).  The 
radial  collimator  was  7.6  cm  in  diameter  and  13.3  cm  long.  The  collimator 
channel  defined  a  truncated  cone  of  entrance  aperture  1.46  cm  in  diameter 
and  exit  aperture  2.16  cm  in  diameter.  Using  a  small  aluminum  collar, 
the  sample  was  positioned  flush  against  the  exit  aperture  of  the  collimator 


Table  2 


22A1(  f  ,2p)2^Na  Reaction 


Threshold 

(MeV) 

24.2 

25Na  Half- 
life  (sec) 

60.0 

25  _ 

Na  Decay 

Constant 

.01155 

Gamma  Ray 

Energies  (MeV) 
and  Intensities 

.975  (14.5%) 
1.612  (9.5*) 
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DETECTOR  SLOTTED  HOLDER 


SAMPLE  CADDY 


Figure  9.  Sample  holder  assemblies, 
(a)  for  irradiation,  (b)  for  counting. 


with  its  central  exis  coincident  with  that  of  the  truncated  cone  defining 
the  collimator  channel.  The  radial  collimator  and  sample-bearing  collar 
were  supported  by  means  of  an  aluminum  assembly  as  shown  in  Figure  9a. 

The  assembly  was  attached  to  a  0.5-in-thick  aluminum  plate.  The  entire 
atructure  was  then  secured  to  the  linac  collimator  by  fitting  holes  in 
the  plate  over  rods  imbedded  in  the  linac  collimator  housing.  These  roda 
were  tapped  so  that  the  plate  was  held  firmly  in  place  by  nuts  screwed 
onto  the  rods.  When  properly  positioned,  the  entrance  aperture  of  the 
radial  collimator  was  located  27.8  cm  from  the  bremsstrahlung  target. 

When  in  place,  the  entrance  face  of  the  sample  was  41.1  cm  from  the  tar¬ 
get,  and  its  central  axis  was  coincident  with  the  central  axis  of  the 
bremsstrahlung  beam. 

Reproducible  sample  positioning  for  radioactivity  counting  was 
alao  imperative.  This  was  accomplished  by  positioning  the  samples  in 
another  sample  holder  which  was  attached  firmly  to  the  detector  as  shown 
in  Figure  9b.  This  assembly,  made  of  plexiglaas,  was  slotted  so  that  a 
sample  caddy  could  be  inserted  into  these  slots  at  some  well-known  dis¬ 
tance  to  the  detector  face.  The  sample  was  placed  into  the  sample  holder 
so  that  the  cylindrical  sample  sxis  was  vertical  and  parallel  to  the  front 
face  of  the  detector.  Using  this  apparatus  the  reproducibility  of  the 
sample  location  with  respect  to  the  detector  was  within  less  than  1  mm. 

Linear  Accelerator 

The  linear  accelerator  was  the  most  important  piece  of  equipment 
used  in  this  study.  The  linac  was  the  source  of  the  bremsstrahlung  beams 
used  in  this  work.  Operational  principles  of  medical  electron  linear 
acceleratora  in  general  are  discussed  in  Appendix  C.  This  section 


66 


describes  the  particular  accelerator  used.  Later,  the  experimental  use 
of  the  accelerator  will  be  discussed.  The  linac  operating  parameters 
which  influence  the  kinetic  energy  achieved  by  the  accelerated  electrons 
are  described  and  their  appropriate  values  are  defined. 

The  accelerator  used  to  produce  the  bremsstrahlung  beams  of  this 
study  was  a  Sagittaire/Therac  40  electron  linear  accelerator  manufactured 
by  CGR-MaV  of  France  belonging  to  the  Radiation  Therapy  Department  of  the 
University  of  Oklahoma  Health  Sciences  Center.  In  its  clinical  mode,  the 
Sagittaire  can  deliver  electron  beams  at  nominal  energies  of  7,  10,  13, 
16,  19,  22,  25,  28,  32,  and  40  MeV.  A  25  MeV  bremsstrahlung  baam  is  also 
available  in  the  clinical  mode.  Six  dose  rates  at  the  isocenter  can  be 
chosen  for  both  the  bremsstrahlung  and  electron  beams:  .5,  1,  2,  3,  4, 
and  10  Gy  per  min. 

*» 

The  Sagittaire/Therac  40  linac  is  shown  schematically  in  Figure 
Cl  (Appendix  C) .  It  essentially  consists  of  three  major  components:  1) 
the  microwave  power  supply  saction,  2)  the  electron  Injector  and  accel¬ 
erator  system,  and  3)  the  gantry  and  treatment  head  component.  In  tha 
modulator  room,  microwave  power  was  provided  by  a  3,000  MHz  oscillator 
and  a  klystron  microwava  amplifier.  In  the  accelerator  room,  electrons 
ara  injected  into  the  acceleration  guides  by  an  electron  gun.  There, 
the  electrons  ara  accalerated  by  the  amplified  3,000  MHz  microwavas, 
traveling  through  two  iris-loaded  cylindrical  wave  guides.  The  energy 
gained  by  the  electrons  is  controlled  primarily  by  the  degree  of  phase¬ 
coupling  between  the  traveling  waves  in  sections  1  and  2.  This  task  Is 
performed  by  a  phase  shifter  consisting  of  a  quartz  rod  which  is  inserted 
into  the  rectangular  wave  guide  feeding  section  1  (see  Appendix  C) . 
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Accelerated  electrons  exit  from  the  second  accelerator  section 
and  then  enter  the  gantry  of  the  accelerator  which  is  located  in  the 
treatment  room.  The  electron  beam  direction  is  controlled  by  two  magnet 
systems  located  in  the  gantry.  The  first  magnet  system  analyzes  the 
electron  beam  energy  by  means  of  a  bending  electromagnet  and  an  energy- 
defining  slit.  The  second  system  consists  of  two  magnetic  beam-bending 
devices  that  redirect  the  electron  beam  to  a  direction  perpendicular  to 
the  acceleration  wave  guides.  The  beam  then  enters  the  therapy  head 
where  removable  bremsstrahlung  targets  may  be  inserted  into  the  electron 
beam  to  produce  x  rays  or  may  be  removed  to  allow  the  electron  beam  to 
emerge.  The  radiation  beam  is  finally  collimated  by  a  set  of  fixed  and 
variable  collimators,  and  its  intensity  is  monitored  by  a  system  of  trans¬ 
mission  ion  chambers  which  also  serves  to  control  the  beam  centering. 

The  accelerator  could  be  operated  in  what  is  called  the' experi¬ 
mental  mode.  Most  of  the  bremsstrahlung  beams  used  in  this  study  were 
produced  In  this  mode.  The  beams  produced  in  this  mode  are  termed  photon 
II  beams  (in  contrast  to  the  clinical  photon  beam  which  is  called  photon 
I) .  Photon  II  beams  differ  from  photon  I  beams  only  in  the  amount  of 
filtration  after  bremsstrahlung  production.  Both  beams  are  produced  in 
identical  4-wm-thick  tungsten  targets  although  each  has  its  own  target. 

The  clinical  photon  I  beam  must  traverse  an  approximately  2.8-cm-thick 
(  antral  dimension)  lead  flattening  filter.  The  photon  II  beams,  on  the 
other  hand,  are  filtered  by  a  2.7-mm-thick  lead  plate  whose  purpose  is 
solely  to  increase  the  amount  of  signal  available  to  the  ionization  cham¬ 
bers  for  beam  steering. 

A  prominent  parameter  which  describes  the  bremsstrahlung  beam 
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of  a  linear  accelerator  Is  its  endpoint  energy.  The  bremsstrahlung  end¬ 
point  energy  is  the  maximum  photon  energy  in  the  photon  spectrum.  It  is 
usually  taken  to  correspond  to  the  mean  kinetic  energy  of  the  electrons 
incident  on  the  bremsstrahlung  target.  This  section  describes  the  accel¬ 
erator  operating  parameters  which  determine  the  final  kinetic  energy  of 
the  electrons.  Also  discussed  is  the  method  employed  in  the  accelerator 
to  define  and  control  the  kinetic  energy  of  the  electrons  prior  to  radia¬ 
tion. 

The  resultant  energy  of  the  accelerated  electrons  is  determined 
by  a  unique  set  of  values  of  the  following  parameters  (see  Appendix  C) : 

1)  the  magnitude  (amplitude)  of  the  axial  electric  field  of  the  radio¬ 
frequency  traveling  waves,  2)  the  relative  phase-coupling  between  the 
traveling  waves  in  the  first  and  second  acceleration  sections,  and  3) 

it 

the  field  strength  of  the  bending  magnet  which  analyzes  the  energy  of  the 
electron  beam.  The  first  parameter  (axial  field  amplitude)  Is  controlled 
by  the  amount  of  radiofrequency  high  voltage  (RFHV)  applied.  The  second 
parameter  (phase -coup ling)  is  controlled  by  the  relative  position  of  the 
quartz  rod  inserted  in  the  wave  guide  that  feeds  section  1.  The  third 
parameter,  analyzing  magnet  field  strength,  is  controlled  by  the  current 
supplied  to  the  energy-analyzing  electromagnet. 

The  energy  of  the  electron  beam  is  ultimately  determined  by  the 
linac  energy -analysis  system.  This  system  defines  the  electron  energy 
to  ±2Z  of  the  mean  energy  (94) .  The  magnetic  analyzer  deflects  the 
accelerated  electrons  through  angles  determined  by  the  electron  kinetic 
energy.  Because  of  the  major  role  played  by  the  energy -analysis  system, 
the  fundamental  principles  of  this  system  require  discussion. 


A  charged  particle  of  velocity  v  entering  a  uniform  magnetic 
field  of  magnetic  induction  f  experiences  a  force  F  perpendicular  to  the 
direction  of  motion  and  to  the  magnetic  induction.  Because  there  is  no 
component  of  this  force  along  the  particle's  direction  of  motion,  the 
particle  speed  is  unchanged  but  the  particle  is  forced  to  move  in  a  cir¬ 
cular  path.  If  a  particle  is  to  move  in  a  circular  path  at  a  constant 
speed,  it  must  be  subjected  to  a  centripetal  force.  Thus,  we  have  ror 
an  electron  of  charge  e,  mass  m^,  and  velocity  v  in  a  magnetic  induction 
t  (58,59) 

2 

m  v 

Bev  •  -  (71) 

where  r  is  the  radius  of  curvature  of  the  circular  path  into  which  the 
electron  is  forced  to  move.  Rearranging  Eq.  (71)  and  substituting  an 
alternate  expression  for  the  momentum  p  (60)  yields 

Berc  -  (T2  +  2mec2T)*S  -  pc  (72) 

where  T  is  the  electron  kinetic  energy  and  c  is  the  speed  of  light. 

According  to  Eq.  (72),  for  a  giver,  constant  magnetic  Induction 
B,  an  electron  will  move  in  a  circular  path  of  radius  r  determined  by  the 
electron  energy  T.  In  the  energy-analysis  system  of  the  accelerator, 
shown  schematically  in  Figure  10,  the  degree  of  deflection  (which  depends 
upon  r)  of  the  electron  produced  by  the  bending-magnet  Induction  depends 
on  the  electron  energy.  Only  electrons  of  kinetic  energy  within  2Z  of 
an  appropriate  mean  energy  will  be  deflected  into  a  radius  of  curvature 
such  that  they  are  allowed  transmission  by  the  energy -defining  slit. 


figure  10.  Electromagnetic  energy  analysing  system. 
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The  appropriate  mean  energy  T  is  that  which  satisfies  the  relation  of 
Eq,  (72).  RF  power  is  automatically  changed  if  the  electron  energy  devi¬ 
ates  from  the  requested  energy  as  determined  by  the  magnetic  induction 
of  the  bending  magnet. 

Electron  kinetic  energy  selection  is  made  possible  by  virtue  of 
the  fact  that  the  degree  of  deflection  of  an  electron  of  a  certain  energy 
in  the  energy-analysis  system  depends  upon  the  magnitude  of  the  magnetic 
induction  of  the  bending  magnet.  This  also  follows  from  Eq.  (72).  The 
magnetic  induction  is  directly  proportional  to  the  amount  of  current  X 
supplied  to  the  electromagnet.  With  this  in  mind,  Eq.  (72)  thus  can  be 
written  in  the  form 

I  -  K(T2  +  2m  c2!)1*  (73) 

.1 

where  K  is  a  function  of  e,  r,  c  and  of  the  specifics  of  the  magnet  sys¬ 
tem.  It  is  constant  as  long  as  saturation  and  residual  magnetization 
effects  are  negligible.  The  exact  form  of  Eq.  (73)  can  be  found  by 
measuring  t'e  energy  of  electrons  transmitted  by  the  energy-analysis 
system  a  function  of  the  current  applied  to  the  bending  magnet,  (kice 
known  in  its  exact  form,  the  equation  can  then  be  used  to  select  certain 
electron  energies  by  choosing  the  appropriate  current  to  the  magnet  of 
the  energy-defining  system. 

Measurements  of  the  type  described  above  have  recently  been  per¬ 
formed  on  the  accelerator.  The  exact  relationship  between  the  electron 
kinetic  energy  and  the  current  to  the  energy-analyzing  magnet  was  found 
and  published  (6,57).  The  authors  used  three  energy  calibration  polnta 
supplied  by  photonuclear  reactions  to  relate  magnet  current  to  electron 
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energy:  the  10.8  MeV  Cu(y,n)  Cu  threshold,  the  17.3  MeV  break  In  the 
1^0(y,n)1^0  yield  curve  and  the  29.2  MeV  ^2S(y,3p)2^Al  threshold.  They 
found  the  deviation  magnet  current  corresponding  to  each  of  the  above 
electron  energies  and  a  function  of  the  form 

I2  =  KXT2  +  K2T  (74) 

and  &2  constants  was  fit  to  the  experimental  points  by  using  a  least- 
squares  procedure.  Note  that  Eq.  (74)  can  be  obtained  by  squaring  Eq. 

(73)  and  redefining  the  constants .  The  resulting  relation  between  magnet 
current  and  electron  kinetic  energy  was 

I  -  (37. 7T2  +  90.3T)*5  (75) 

From  this  relation,  shown  in  Figure  11,  one  can  obtain  the  magnet  current 
which  is  necessary  to  produce  an  electron  beam  of  a  desired  kinetic  energy. 

The  current  to  the  energy-defining  electromagnet  is  monitored  by 
measuring  the  voltage  drop  across  a  sampling  resistor,  as  shown  in  Figure 
12.  In  this  fashion,  the  reading  of  a  digital  voltmeter  (DVM)  serves  as 
an  index  of  magnet  current.  In  the  accelerator  experimental  mode,  the 
current  to  the  analyzing  magnet  is  continuously  adjustable  by  means  of  a 
ten-turn  helipot  located  at  a  convenient  point  in  an  electronics  rack 
which  in  turn  is  located  in  the  modulator  room.  In  the  clinical  mode, 
the  current  to  the  electromagnet  is  fixed  at  the  value  appropriate  for 
the  clinical  beam  to  be  used.  The  position  of  the  switch  S  of  Figure  12 
determines  whether  the  linac  is  operated  in  the  experimental  mode  or  in 
the  clinical  mode. 

In  the  clinical  mode,  when  one  of  the  ten  available  energies  is 


Figure  11.  Relation  between  mean  electron  kinetic  energy 
and  deviation  current  to  the  bending  magnets  of  the  energy  - 
analysis  system  of  the  Sagittaire  accele.ntor.  From  Ref.  (57). 
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selected,  the  phase -coup ling  and  analyzing  magnet  current  are  set  auto¬ 
matically  and  the  RFHV  is  controlled  by  a  servo  mechanism.  In  the  experi¬ 
mental  mode,  the  energy  selection  switches  are  bypassed  and  the  automatic 
servo  mechanisms  regulating  the  energy  control  parameters  are  disabled  so 
that  the  energy-defining  magnet  current,  phase-shifter  position,  and 
RFHV  can  be  set  and  maintained  manually.  A  bremsstrahlung  beam  of  a 
desired  endpoint  energy  is  obtained  by  first  selecting  the  appropriate 
energy-analyzing  magnet  current  (as  measured  by  the  DVM)  corresponding 
to  the  desired  endpoint  energy.  The  appropriate  phase  and  RFHV  values 
are  then  set.  The  proper  settings  had  been  obtained  earlier  (57)  by 
keeping  to  a  minimum  the  gun  current  required  to  keep  the  selected  end¬ 
point  energy  beam  at  a  constant  given  intensity.  The  appropriate  linear 
accelerator  operating  parameters  as  a  function  of  bremsstrahlung  beam 
endpoint  energy  are  presented  in  Table  3. 

Counting  System 

Sample  activity  was  measured  with  a  high-purity  germanium  (HPGe) 
semiconductor  detector  gamma-ray  spectroscopy  system.  These  systems  are 
bost  known  for  their  excellent  energy  resolution.  This  resolution  enables 
the  separation  of  interference  peaks  from  the  full  energy  (total  absorp¬ 
tion)  peak.  The  system  used  in  the  present  work  for  activity  determina¬ 
tion  possessed  an  energy  resolution  of  0.3X  FWHM  at  1332.5  keV.  A  basic 
description  of  the  principles  and  characteristics  of  semiconductor  spec¬ 
trometer  systems  is  provided  in  Appendix  D.  The  present  section  is  a 
description  of  the  particular  system  used  in  this  study.  Discussed  are 
the  specifications  of  each  component  of  the  system  and  the  system  perfor¬ 


mance  . 
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Table  3 


Linear  Accelerator  Operating  Parameters 


8 


The  detector  was  a  closed-end  coaxial  type  high-purity  germanium 
detector  manufactured  by  Princeton  Gamma-Tech  (PGT) .  It  was  a  cylinder 
4.0  cm  long  and  4.3  cm  diameter  with  an  entrance  window  of  aluminum  0.1 
cm  thick  (61) .  The  detector  was  used  with  a  2100-volt  positive  bias  from 
a  high-voltage  power  supply  and  was  cooled  to  77°K.  A  copper  cold  finger, 
one  end  touching  the  detector  and  the  other  immersed  in  liquid  nitrogen, 
provided  the  conduction  path  for  heat  flow.  Detector  cooling  waa  neces¬ 
sary  to  reduce  thermal  noise.  The  detector  was  shielded  with  10.2  cm 
lead  to  reduce  the  background  radiation  contribution. 

The  detector  output  charge  pulse,  which  is  proportional  to  the 
number  of  electron-hole  pairs  produced  by  a  photon  interaction  in  the 
detector,  was  fed  into  a  PGT  model  RG-11  preamplifier.  The  preamplifier 
was  located  directly  adjacent  to  the  germanium  detector  to  minimize 
pick-up  of  stray  fields.  The  preamplifier  amplified  the  charge  pulae 
amplitude  by  a  factor  of  25,000.  The  output  voltage-pulse  of  the  pream¬ 
plifier  was  shaped  to  have  an  approximately  10  naec  rise  time  and  50  ysec 
decay  time  (61) . 

The  preamplifier  output  pulse  was  processed  by  a  Canberra  Model 
1413  amplifier.  This  amplifier  not  only  increased  the  amplitude  of  the 
signal  pulse  but  also  modified  the  pulse  to  a  Gaussian  shape,  thereby 
increasing  the  signal-to-noise  ratio  and  making  the  output  pulse  compat¬ 
ible  with  the  analog -to-dlgital  converter  (ADC)  of  the  multichannel 
analyzer.  The  spectroscopic  amplifier  also  had  baseline-restorer  circuits 
to  minimize  the  effects  of  baseline  fluctuations. 

The  amplitude  of  the  pulse  from  the  amplifier,  which  was  propor¬ 
tional  to  the  amplitude  of  the  detector  original  charge  pulse,  was  digi- 


Cized  by  the  ADC  of  a  Canberra  Model  4100  multichannel  analyzer  (MCA) . 

The  50  MHz  digitizing  rate  was  under  crystal  clock  control.  The  digitized 
value  of  the  pulse  amplitude  determined  the  location  of  the  analyzer 
memory  to  be  incremented.  Thus,  the  address  of  the  channel  in  the  array 
of  the  MCA  memory  locations  was  directly  proportional  to  the  radiation 
energy  absorbed  in  the  detector  (see  Appendix  D) . 

The  relationship  between  the  channel  number,  which  corresponds 
to  a  given  pulse-height  interval,  and  the  energy  deposited  in  the  detec¬ 
tor  by  a  photon  interaction  event  was  found  by  energy-calibrating  the 
pulse-height  scale  with  11  gamma-ray  emitters  having  a  total  of  24  gamma 
rays  of  known  energies  (Table  4) .  The  energy  calibration  was  performed 
so  that  a  2.5  MeV  energy  range  spanned  the  MCA's  1024  channels.  A 

straight  line  was  fit  to  the  data  of  Table  4  using  linear  regression  to 

■»« 

obtain  the  relation  between  photon  energy  and  channel  number. 

E  -  (2.47) C  +  9.74  (76) 

where  E  is  the  photon  energy  in  keV  ana  C  is  the  channel  number.  The 

energy -channel  relation  was  indeed  linear  with  a  coefficient  of  correla- 
2 

tion  of  r  ■  1.00. 

The  counting-system  efficiency  as  a  function  of  photon  energy 
was  investigated  using  gamma-ray  emitters  calibrated  by  the  National 
Bureau  of  Standards  (NBS) .  The  NBS -calibrated  sources  effectively  were 
point  sources  (minute  quantities  of  isotope  deposited  on  thin  plastic 
tape).  The  characteristics  of  these  NBS -calibrated  sources  are  presented 
in  Table  5.  Intensities  shown  in  the  table  are  source  intensities  at  the 


time  of  calibration. 
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Table  4 


Sources  used  In  Counting  System  Energy  Scale  Calibration 


Soui ces 

Gamma  ray 

Energy  (keV) 

Channel 

Number 

Eu-lf>2 

121.8 

46 

344.3 

135 

778.9 

311 

964.0 

386 

1112.1 

446 

1408.0 

566 

Ba-133 

276.4 

108 

302.8 

118 

356J0 

141 

383.8 

151 

Hg-203 

279.2 

109 

Sb-125 

427.0 

169 

463.0 

184 

599.0 

239 

Cs-137 

661.6 

264 

ltr-85 

514.0 

204 

Bl-207 

569.6 

227 

1063.4 

427 

1769.7 

713 

Co-58 

810.7 

324 

Mn-54 

834.8 

334 

Co-60 

1173.2 

471 

1332.5 

536 

Na-22 

1274.6 

512 
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The  point  sources  of  Tsble  5  were  counted  st  2.1  cm  and  8.2  cm 
from  the  detector  face.  The  photopeak  efficiency  as  a  function  of  gamma- 
ray  energy  was  calculsted  by  dividing  the  net  photopeak  count  rate  (in 
counts  per  aec)  by  the  rate  of  gamma  rays  of  spproprlste  energy  emitted 
by  the  source  (in  gamma  rsys  per  sec) .  Net  photopeak  counts  were  obtained 
by  adding  the  total  number  of  counts  in  the  pulse-height  analyzer  channels 
encompassing  the  entire  photopeak  and  then  subtracting  the  background 
contribution.  Background  counts  were  determined  by  first  averaging  the 
counts  in  the  channels  where  the  photopeak  blends  into  the  background 
continuum  and  then  multiplying  the  resultant  average  background  counts 
per  channel  times  the  total  number  of  channels  in  the  photopeak.  The 
gamma-ray  emission  rates  of  the  sources  were  computed  from  the  source 
activities  as  given  by  the  NBS  certlficstes  of  calibration.  Bsrlum-133 

■ti 

gamma- ray  intensities  were  corrected  for  coincidence-summing  effects  (62)  . 

The  count in g-sys tern  point-source  ef f iciencies ,  at  2.1  and  8.2 
cm,  are  presented  in  Table  6  and  Figure  13  as  a  function  of  photon  energy. 
Efficiency  has  been  defined  as  the  ratio  of  the  net  total  photopeak  count 
rate  to  the  gamma-ray  emission  rate.  The  efficiency,  hence,  is  in  units 
of  counts  per  second  per  gamma  ray  emitted  per  sec. 

The  smooth  curves  of  Figure  13  were  hand-drawn  through  the  exper¬ 
imental  points.  The  curve  shapes  are  in  good  agreement  with  typical  effi¬ 
ciency  curves  for  thia  type  detector  (63) . 

Timing  Control  Apparatus 

The  amount  of  radioactivity  produced  by  photonuclear  reactions 
in  the  sample  and  aubsequently  measured  in  the  counting  system  depended 
very  critically  on  three  separate  time  intervals  which  had  to  be  fixed 
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Energy 

(keV) 

Source 

*  Photopeak  Efficiency  (X10^) 
at  2.1  cm  at  8.2  cm 

276.4 

Ba-133 

23.7 

4.48 

302.9 

Ba-133 

20.6 

3.92 

356.0 

Ba-133 

16.9 

3.22  • 

383.9 

Ba-133 

15.6 

2.82 

■  ^ 

514.0 

Kr-85 

11.2 

2.08 

569.6 

Bl-207 

10.6 

1.89 

661.6 

Ca-137 

9.2 

1.67 

1063.4 

Bi-207 

5.0 

0.98 

1173.2 

Co-60 

4.9 

0.94 

1332.5 

Co-60 

4.3 

0.82 

1769.7 

Bi-207 

3.0 

0.55 

♦  M 
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precisely.  These  were  che  Irradiation  time,  the  time  required  for  the 
transfer  of  the  sample  from  the  linear  accelerator  to  the  counter,  and 
the  counting  time.  The  counting  time  was  controlled  by  the  present  live¬ 
time  counter  in  the  multichannel  analyzer.  The  irradiation  and  transfer 
times  were  controlled  by  a  timing  control  apparatus  described  in  the  next 
paragraph . 

The  timing  control  apparatus  consisted  of  a  module  containing 
two  precision  electrical  clocks  which  could  be  preset.  The  timer  module 
was  Interfaced  with  the  linac  controls  and  with  a  start-count  gate  on  the 
MCA.  The  irradiation  time  was  set  on  timer  1  and  the  sample  transfer 
time  was  set  on  timer  2.  When  the  linac  was  turned  on,  timer  1  was  acti¬ 
vated  and  the  irradiation  comnenced.  When  the  time  set  on  clock  1  ran 
out,  the  linac  was  automatically  shut  off  and  a  relay  started  timer  2 
which  ran  while  the  sample  was  transferred  from  the  linac  to  the  detector. 
When  timer  2  shut  off,  a  signal  pulse  was  sent  to  the  MCA  which  started 
the  counting  period.  This  continued  until  the  MCA  live-time  indicator 
reached  the  preset  value.  With  these  time-controlling  devices,  all  times 
were  set  prior  to  irradiation  and  all  phases  of  the  experiment  advanced 
automatically. 


P-2  Ionization  Chamber 

Normalized  yield  determination  required  measurements  of  the 
energy  fluence  of  the  brenisstrahlung  beam  traversing  the  sample.  This 
was  accomplished  with  an  NBS  type  P-2  ionization  chamber,  described  in 
detail  elsewhere  (56),  which  was  located  approximately  2.4  meters  from 
the  sample.  The  activity  induced  in  the  sample  was  divided  by  the  response 
of  this  ionization  chamber  to  compute  the  normalized  yield. 
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The  experimental  arrangement  cf  the  ionization  chamber  and  its 
associated  components  is  shown  schematically  in  Figure  14.  The  ionization 
chamber  (I)  collecting  potential  was  1200  volts  and  was  supplied  by  a  DC 
high-voltage  power  supply  (HV) .  The  output  current  of  the  chamber  was 
Integrated  on  a  high-quality  10.096  uF  polystyrene  capacitor  (C) .  A  var¬ 
iable  resistor  (R)  was  Introduced  in  the  circuit  so  that  the  RC  constant 
was  Inversely  proportional  to  the  decay  constant  of  the  radioactive 
species  being  produced  (see  Measurement  of  Yield  section  in  Chapter  II) . 
The  potential  across  the  capacitor  was  measured  with  a  Cary  model  401 
vibrating  reed  electrometer  (E) ,  which  was  located  in  the  modulator  room 
of  the  accelerator.  A  Varian  model  401  preamplifier  (P)  provided  initial 
amplification  and  impedance-matching. 

The  response  of  the  ionization  chamber  to  bremsstrahlung  energy 
fluence  was  neerly  independent  of  bremsstrahlung  endpoint  energy.  Celi- 
bration  of  the  chamber  showed  a  response  of  approximately  2.4  micro-cou¬ 
lombs  per  joule  of  incident  energy  of  ±0.9Z  from  15  to  35  MeV  (56).  The 
measurement  system  was  checked  for  leakage  current,  possible  recombination 
losees,  and  contributions  to  the  ionizetion  current  from  scettered  radia¬ 
tion  (57).  All  effects  were  found  to  be  negligible  under  the  experimental 
conditions  of  this  work.  System  reproducibility  was  checked  at  the  begin¬ 
ning  of  each  data-collection  session  by  measuring  the  charge  collected 
226 

when  a  30  mg  Ra  source  was  placed  et  e  specific  location  on  the  P-2 
chamber.  These  measurements  showed  e  verieblllty  of  ebout  1.2Z  (standard 
error)  over  the  course  of  the  entire  experiment  thus  Indicating  good 
stability. 


Figure  14.  Radiation  monitor  system  arrangement. 

I  -  Ionization  Chamber.  HV  -  High  Voltage  Supply. 
C  -  Capacitor.  R  -  Resistor.  P  -  Preamplifier. 

E  -  Electrometer. 
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Methods 

Fhotonudear  Ratios 

This  section  describes  the  procedures  followed  to  acquire  and 
analyze  the  experimental  data  used  to  determine  the  photon uclear  ratios. 

In  general  the  experimental  method  consisted  of  irradiating  the  K^SIF^ 
and  teflon  samples  previously  described  with  bremsstrahlung  of  a  certain 
endpoint  energy,  counting  the  samples,  and  then  computing  the  ratios  as 
a  function  of  endpoint  energy.  This  procedure  was  repeated  several  times 
so  that  a  mean  ratio  with  a  minimum  amount  of  uncertainty  at  each  endpoint 
energy  could  be  reported.  Because  two  sets  of  ratios  were  determined 
using  two  different  sample  types,  the  experimental  method  for  each  dif¬ 
fer*  1  slightly.  Each  method  Is  described  separately  in  the  following 
after  a  general  description  of  the  irradiation  and  tha  radioactivity 
determination  procedures  which  were  common  to  both  seta. 

The  samples  were  irradiated,  in  the  geometry  shown  In  Figure  15, 
with  the  Saglttalre  linear  accelerator  operating  in  the  photon  II,  or 
experimental  mode.  The  operation  and  characteristics  of  the  linac  exper¬ 
imental  mode  have  been  set  forth  in  the  description  of  tha  Sagittaire  in 
tha  Materials  section  of  this  chapter.  The  X^SiF^  were  exposed 

to  bremsstrahlung  of  endpoint  energies  from  13  to  23  MeV  at  integral  MeV 
Intervals.  The  teflon  samples  were  exposed  to  x-ray  beams  of  endpoint 
energies  ranging  from  20  to  32  MeV,  also  at  integral-MeV  intervals. 

Prior  to  sample  irradiation,  the  bremsstrahlung  beam  endpoint 
energy  was  chosen  at  random.  Continuous  selection  of  endpoint  energies 
was  made  possible  by  the  linac  modification  previously  described.  The 
deviation  magnet  current  corresponding  to  the  chosen  endpoint  energy  was 


Figure  IS.  Sample  irradiation  geometry. 
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set*  This  setting  established  the  accelerated  electron  beam  energy.  The 
optimum  phase-shifter  position,  gun  current  and  RF  high  voltage  (RFH.V) 
corresponding  to  the  chosen  endpoint  energy  (Table  3}  were  also  set.  In 
addition,  the  irradiation  and  transfer  times  were  set  on  the  timing  con¬ 
trol  module.  The  sample  was  then  inserted  in  the  llnac  sample  holder. 

»  The  linac  was  started  and  irradiation  commenced. 

At  the  end  of  the  irradiation,  the  monitor  units  utilized  were 
noted  and  recorded.  The  sample  was  then  transferred  to  the  counting  room 
and  was  Inserted  in  the  sample  holder  near  the  germanium  detector.  After 
expiration  of  the  transfer,  or  transit,  time  which  was  set  on  timer  2  of 
the  timing  control  module,  the  sample  was  counted.  The  background-correc¬ 
ted  counts  were  then  used  to  compute  the  photonuclear  ratio  at  the  selac- 
ted  endpoint  energy. 

The  photonuclear  ratio  (FR)  has  been  defined  as  the  number  of 
photonuclear  reactions  per  target  nucleus  Induced  in  element  1  of  the 
eample  compound  divided  by  the  number  of  photonuclear  reactions  per  tar¬ 
get  nucleus  indues d  in  element  2  of  the  compound.  To  determine  tha  number 
of  photonuclear  reactions,  the  activity  of  the  samples  was  measured. 

In  both  target  eleaente  of  both  samplee,  the  radioisotopee  of 

interest  dacayed  by  positron  amiaeioo  with  the  subsequent  emission  of 

0.511  MeV  annihilation  gamma  rays.  It  was  the  0.511  MeV  gamma  raye  that 

were  counted  to  detarmina  the  amount  of  radioisotope  ,  codured.  Since  the 

both  of  the  radioisotopes  produced  in  the  compound  differed 

(7.63  min  38K  vs.  109.8  min  in  KjSiFg  and  20.3  min  11C 
18 

F  in  teflon),  the  sample  was  counted  twice  to  enable 
of  the  activity  of  each  radioisotope. 


half-lives  of 
significantly 
vs.  109.8  min 
determination 
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The  first  count  contained  contributions  from  both  the  short-lived 
and  the  long-lived  radioisotopes.  The  second  count  was  delayed  in  time 
until  tha  shorter-lived  radioisotope  had  decayed  very  nearly  completely 
(ten  half-lives  of  the  shorter-lived  isotope) .  The  second  count ,  there¬ 
fore,  was  due  almost  exclusively  to  the  longer-lived  radioisotope.  Based 
on  tha  second  count,  utilizing  appropriate  decay  factors,  the  contribution 
of  the  longer-lived  isotope  to  the  first  count  was  calculated.  This  valua 
was  then  subtracted  from  the  number  of  counts  measured  in  the  first  count 
to  obtain  the  contribution  due  to  the  decay  of  the  shortar-lived  isotopa. 

The  procedure  described  above  for  determining  the  activity  of 
each  isotopa  of  interest  in  the  sample  was  followed  with  both  the  KjSiFg 
and  the  teflon  samples.  In  the  K^SiFg  experiment  the  first  count  was  due 

3O  1g 

to  the  activity  of  both  x  and  F,  and  tha  second  count,  obtained  after 

3a*  decay,  was  due  almost  entirely  to  the  activity.  In  the  taflon 

experiment,  the  first  measured  count  contained  contributions  from  both 

the  X1C  and  1®F  activities  while  the  second  count,  obtained  after  the 

18 

had  decayed,  was  due  almost  entirely  to  the  F  activity. 

The  counts  mentioned  in  tha  foregoing  were  obtained  from  the 
total  absorption  peak  or  photopeak.  The  gross  counts  in  the  photopeak 
were  corrected  for  background  to  yield  net  photopeak  counts  in  a  manner 
outlined  below.  The  nuaber  of  counts  in  all  channels  encompassing  the 
photopeak  were  summed.  Background  counts  were  determined  by  first  noting 
the  nuaber  of  counts  in  the  channels  bordering  the  photopeak  whare  the 
peak  blended  into  the  background  continuum.  The  counts  in  these  channels 
were  then  averaged  and  multiplied  by  the  nuaber  cf  channels  in  the  photo- 
peak  to  yield  the  background  count.  The  background  count  was  then  sub- 


92 


tree ted  from  the  sum  of  the  totel  photopeak  counts  to  yield  the  net  photo-* 
peak  counts. 

The  K^SiFg  samples  were  irrediated  to  bremsstrahlung  of  endpoint 
energy  ranging  from  13  to  23  MeV.  One  by  one,  eech  sample  wes  positioned 
In  the  llnac  sample  holder  and  wes  irradiated,  as  described  previously, 
for  2  minutes  at  a  certain  endpoint  energy.  The  sample  was  then  trans¬ 
ferred  to  the  counting  system  and  allowed  to  stand  for  10  minutes.  The 
delay  provided  for  the  decay  of  any  Interfering  short-lived  positron  acti¬ 
vity  such  as  6,34  sec  ^Al,  4.14  sec  ^Sl,  .93  sec  ^®inK.  The  10-min  wait 
also  served  to  reduce  the  amount  of  detector  deed  time.  After  the  10-min 

wait  the  sample  was  counted  for  1  minute.  The  counts  -  the  0.511  MeV 

38  18 

photopeak,  due  to  the  activity  of  both  K  and  F,  were  then  recorded. 

The  sample  was  then  recounted  for  1  min  90  minutes  after  the  conclusion 

of  ita  irradiation  and  again  the  0.511  MeV  photopeak  counts  were’ recorded. 

18 

It  was  assumed  that  the  90-min  delayed  count  contained  only  F  counts 

since  11.8  half-lives  had  passed  leaving  only  .03Z  of  the  original 

activity.  The  entire  procedure,  which  constitutes  what  is  termed  a  "no," 

was  then  repeated  for  another  sample  at  another  endpoint  energy. 

The  ^SiFg  photonudear  ratio  at  a  given  endpoint  energy  obtained 

18 

in  a  given  run  is  calculated  as  described  below.  The  net  F  count  10 
minutes  after  irradiation  is  calculated  from  the  90-min  count  through  the 
relation 


FC10 


(77) 


where  FC^g  and  FC^g  are  the  10-min  and  90-min  counts  corresponding  to  the 
18F  activity  at  10  min  and  90  min  respectively.  X  is  the  decay  con- 
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stant  (.00631  min  ).  The  net  10-tain  K  count  is  obtained  by  subtracting 

the  net  *®F  cotint  at  10  min  from  tha  measured  plus  1®F)  10-min  count. 

The  net  38k  and  net  *®F  counts  are  then  substituted  Into  Eq.  (70)  to 

obtain  tha  appropriate  photonuclear  ratio. 

The  teflon  samples  wera  irradiated  with  breosstrahlung  of  andpoint 

energy  ranging  from  20  MeV  to  32  MeV.  Each  sample  was  positioned,  ona  by 

one,  in  the  linac  sample  holder  and  was  irradiated  for  3  min  at  a  selected 

endpoint  energy.  The  sample  was  then  transferred  to  the  counting  system 

and  allowed  to  stand  for  10  min  to  reduce  the  amount  of  detector  dead 

time.  Aftar  the  10-min  wait  the  sampla  was  counted  for  1.67  min  (100 

sec).  The  counts  in  the  0.511  MeV  photopeak,  due  to  the  decay  of  both 
11  18 

C  and  F,  were  then  recorded.  The  sample  was  recounted  for  100  sec 
200  minutes  aftar  the  conclusion  of  its  irradiation  and  again  the  0.511 

ti 

MeV  photopeak  counts  were  recorded.  The  200-mln  count  contained  almost 
18  11 

exclusively  F  counts  since  9.85  C  half-lives  had  transpired  leaving 
only  .11%  of  the  original  activity.  The  entire  procedure,  or  run,  was 
then  repeated  for  another  eaapla  at  another  endpoint  energy. 

For  a  given  run,  the  teflon  photonuclear  ratio  at  a  given  end¬ 
point  energy  was  calculated  in  a  fashion  similar  to  that  used  to  calc’ilate 

18 

tn.  IjSIPj  ratios.  The  nat  P  count  10  min  after  irradiation  was  calcu¬ 
lated  from  the  2 00 -min  count  by 


FC  -  FC 

TC10  FC200* 


(78) 


where  and  FCjqq  are  the  10-min  and  200-mln  counts  due  to  *®F  activity 
at  10  min  end  200  min  respectively.  Agein  X  is  the  decay  constant 
(.00631  min  *).  The  nat  10-min  count  is  obtained  by  subtracting  the 


net 


F  10-min  count  from  the  measured  (^C  plus  10-min  count.  The 

11  IS 

net  C  and  net  F  counts  are  then  substituted  Into  Eq.  (70)  to  obtain 
the  appropriate  photonuclear  ratio. 

The  Bremsstrahlung  Spectrum 

This  section  describes  the  numerical  analysis  method  which  was 
used  to  obtain  a  functional  form  for  the  bremsstrahlung  spectra  produced 
by  the  Saglttaire  accelerator.  Thi  computational  procedure  consisted 
essentially  of  tha  following  (re far  to  Fig.  16):  1)  the  targat  (thick¬ 

ness  D)  was  divided  Into  n  slabs  of  thickness  Ad;  2)  an  intrinsic  spec¬ 
trum  was  assumed  for  each  slat;  3)  the  photons  radiated  were  attenuated 
by  remaining  slabs;  4)  the  electron  energy  in  tha  ittl  slab,  T^,  was 
degraded  by  interactions  in  preceding  slabs;  5)  the  number  of  elactrona 
reaching  subsequent  slabs  was  decreased;  6)  contributions  from  each  slab 
were  added  to  form  a  composite  spectrum. 

2 

The  Saglttaire  accelerator  target,  7.72  g/cm  of  tungsten,  was 

2 

divided  into  200  slabs,  each  of  thickness  0.0386  g/cm  (approximately 
0.006  radiation  lengths).  The  bremsstrahlung  cross  section  used  to  com¬ 
pute  the  intrinsic  spectrum  from  each  slab  was  Schiff 'a  cross  section 
differential  in  photon  emission  tingle  and  energy  [Eq.  (36)].  The  formula 
employed  was  in  the  notation  of  Hisdal  (54)  in  the  rewritten  form  of  Lent 
and  Dickinson  (52).  For  each  slab,  the  intrinsic  spectrum  was  integrated 
over  photon  emission  angles  up  to  that  determined  by  electron  multiple 
scattering  in  the  slab  (54). 

The  reason  for  choosing  this  form  of  intrinsic  spectrin  had  to 
do  with  the  chosen  slab  thickness  in  combination  with  the  solid  angle 
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Figure  16.  Thick  target  parameter*  uaed  in  the 
breasetrahlung  computer  program. 
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subtended  by  the  sample  irradiated.  The  slab  thickness  had  been  chosen 
such  that  the  electron  was  multiple-scattered  prior  to  radiation.  The 
solid  angle  subtended  by  the  sample  was  such  that  the  actual  spectrum 
seen  by  the  sample  contained  photons  emitted  at  angles  up  to  those  deter¬ 
mined  by  multiple-acattaring.  Under  thase  circumstances,  the  differential 
spectrum  integrated  over  photon  emission  angles  up  to  the  multiple- 
scattering  angle  represented  an  approprlata  approximation  to  the  actual 
spectrum  to  which  the  sampla  was  exposed  (19,54,64,65). 

The  radiated  photons  were  attenuated  by  remaining  slabs.  A 
factor  of  the  form 


^(hv)  -  expt-UjftvHD  -  dj)]  (79) 


was  used  for  this  purpose.  In  Eq.  (79),  a^Ow)  is  ths  photon  attenuation 
factor  as  a  function  of  photon  snergy  for  slab  i,  D  is  the  total’ target 
thickness  and  d^  la  the  depth  of  tha  present  slab,  i^fhv),  the  total 
attenuation  coefficient  of  tungsten  as  a  function  of  photon  energy  hv, 
was  obtained  by  fitting  simple  quadratic  expressions  to  tha  total  atten¬ 
uation  coefficient  for  tungsten  given  by  Storm  aud  Israel  (66). 

Tha  average  electron  kinetic  snargy  in  sach  slab  T^  was  deter¬ 
mined  by  energy  lossss  in  previous  slabs 
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In  Eqs.  (80),  T^  is  the  average  electron  kinetic  energy  In  the  present 


slab,  *-s  the  electron  kinetic  energy  In  the  previous  slab,  Tq  Is  the 

Incident  electron  kinetic  energy  (endpoint  energy),  Ad  Is  the  slab  thick¬ 
ness,  end  dT/ds  Is  the  total  mass  stopping  power  of  electrons  in  tungsten. 
The  totel  mass  stopping  power  was  obtained  by  fitting  simple  quadratic 
expreeslons  to  the  mass  stopping  power  data  of  Berger  and  Seltzer  (67). 
The  angle  O-  Is  the  r.m.e.  angle  of  multiple-scattering  according  to  the 
Mcllere  theory  155).  liie  angle  of  multiple-scattering  was  calculated 
using  Blchsel's  expression  (68) 


M57)(Z*  +  Z) (B  -  5. 01) Ad 

+  imT^) 


(81a) 


In  Eqs.  (81),  Z  and  A  ere  the  atomic  number  end  atomic  mass  of  the  terget, 
respectively,  is  the  electron  kinetic  energy  in  the  1  slab  (Eqs.  80), 
in  is  the  reet  energy  of  the  electron  In  MeV,  Ad  Is  the  elab  thickness, 
and  B  is  the  Mollere  B  (68).  The  value  of  5.01  was  determined  In  sn 
experiment  (deecrlbed  leter  In  this  chapter)  In  which  the  engular  distri¬ 
bution  of  the  bremaetrahlung  beam  of  the  Saglttelre  was  lnvestigeted  to 
eetimate  the  electron  multlple-sce tterlng  dletributlon  in  the  target. 

In  a  elab,  the  e.tual  angle  Into  which  an  electron  Is  vultlple- 
ecettered  depends  also  on  the  multlple-ecettering  angle  of  the  previous 
elab  (53).  Thus,  for  e  elab  1,  the  angle  of  milt lple -scattering  C-  was 
obtained  from  the  reletlon 

*2  -  e-2^  +  *1  (8ib) 

2  2 

where  6^  Is  given  by  Eq.  (81a)  and  6^  ie  the  square  of  the  angle  of 
multlple-sce tterlng  of  the  previous  eleb.  The  angles  were  summed  In 
quad return  because  they  essentially  represent  meesuree  of  dispersion. 


In  latter  slabs  of  the  target,  the  electrons  undergo  a  diffusion  process 

2 

In  which  the  electron  angular  distribution  remains  approximately  cos  6- 
2 

for  which  ©■  -  0.5388  radians  (69).  In  these  slabs,  {as  given  by  Eq. 
(81b)]  was  limited  to  0.734  radians  (53). 

The  fraction  of  the  incident  electron  flux  which  reached  the  itb 
slab  was  calculated  with  the  following  formulas  due  to  Ebert  (70). 

r±  -  exp[-a(Ad/R)b]  (82) 

e  -  (1  -  l/b)1-b  (82a) 


b  -  [387TQ/1.000075ZTQV25  (82b) 


In  Eqs.  (82),  is  the  frection  of  the  Incident  electron  flux  that 
reaches  the  i*^  slab  and  U  is  the  extrapoleted  range.  All  remaining 
variables  have  been  defined  previously.  These  equations  resulted  from 
curves  fit  to  experimental  deta  for  monoenergetlc  electrons  up  to  12  MeV 
In  foils  of  select  elements  up  to  uranium  (70) .  The  assumption  was  made 
here  that  these  relations  hold  up  tc  about  30  MeV. 

The  final  photon  spectrum  was  obtained  by  summing  the  photon  con¬ 
tributions  from  each  slab. 

n 

N(hv)  •  I!  N. (hv)i.  (83) 
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In  Eq.  (83),  H(hv)  is  the  total  number  of  photons  of  energy  hv  in  the 


spectrum.  The  summation  is  over  all  slabs,  Is  the  transmission  factor 
of  the  l**1  slab  [Eqs.  (82)],  and  N^(hv)  Is  the  number  of  photons  of  eo-Vgy 
hv  contributed  by  the  l**1  slab 

M^Chv)  "  doChvJn^Chv)  .  (83a) 

In  Eq.  (83a),  do(hv)  is  the  bremsstrahlung  cross  section  In  cm  [Eqs. 
(38)],  o^(hv)  Is  the  photon  attenuation  factor  as  a  function  of  photon 

£t| 

energy  for  the  1  slab  (Eq.  (79)],  and  n^  Is  the  number  of  target  nuclei 

2  th 

per  cm  In  the  1  slab 


!|a  Ad 

A  cos&£ 


(83b) 


where  la  Avogadro's  number,  A  Is  the  atomic  mass  of  the  target,  p  la 
the  density  of  the  target.  Ad  Is  the  slab  thickness,  and  Is  the  r.m.a. 
angle  of  multiple-scattering  (Eqs.  (81)]. 

In  most  linear  accelerators,  the  electrons  producing  bremsstrah¬ 
lung  possess  a  finite  distribution  of  energies  centered  about  a  mean  elec¬ 
tron  energy.  In  the  case  of  the  Sagittaire  accelerator,  electrons  of 
energies  within  ±21  of  a  seen  electron  energy  are  transmitted  by  the 
energy-defining  slit  (94).  The  bremsstrahlung  computer  program  was  modi- 
flad  to  account  for  this  fact.  At  any  given  endpoint  energy  Tq,  the 
spectrum  was  assumed  to  consist  of  photons  produced  by  electrons  of  ini¬ 
tial  energy  T  ,  T  -  2XT  ,  and  T  +  2ZT  .  The  endpoint  energy  dlstrlbu- 
o  o  o  o  o 

tion  was  assumed  to  be  a  Gaussian  of  4ZT  full-width -half -maximum.  Rale- 

o 

tive  weights  of  .25,  .50,  and  .25  thus  were  given  to  the  spectra  produced 
by  electrons  of  energy  Tq  -  21,  T  ,  and  Tq  +  21,  respectively. 


A  computer  program,  written  In  Fortran  and  incorporating  the 
above,  expressions,  was  used  to  produce  the  bremsstrahlung  spectra  approxi¬ 
mations  required  for  this  study.  A  listing  of  this  program  is  provided 
in  Appendix  E. 


Electron  Multiple-Scstterlng 

High-energy  bremsstrahlung  la  emitted  at  small  angles  with 
respect  to  the  direction  of  the  incident  electron.  An  order  of  magnitude 
estimate  of  this  angle  is  (19) 

+  -  u/Eq  radians  (84) 

where  t  Is  the  angle  between  the  direction  of  the  radiation  photon  and  the 
initial  direction  of  the  radiating  electron,  u  is  the  rest  energy  of  the 
electron,  and  Eq  is  the  total  electron  energy.  For  a  bremsstrahlung  beam 
of  25  MeV  endpoint,  t  is  of  the  order  of  0.02  radians.  In  an  experiment 
where  the  angular  distribution  of  high-energy  bremsstrahlung  intensity  is 
measured,  the  electron  multiple-scattering  distribution  csn  be  estimated 
if  the  assumption  is  made  that  bremsstrahlung  is  emitted  st  angles  small 
compared  to  the  angles  of  electron  scattering.  Such  an  experiment  was 
performed  to  examine  the  ability  of  Eqs.  (81)  of  appropriately  estimating 
the  angle  of  multiple-scattering. 

Eleven  teflon  samples  were  placed  in  a  strsight  line  on  a  plane 
parallel  to  and  75  cm  from  the  bremsstrahlung  target.  The  samples  were 
spaced  2  cm  apart  with  the  central  sample  located  on  the  central  axis  of 
the  bremsstrahlung  beam.  The  samples  were  exposed  to  a  25  MeV  photon  II 
beam  for  5  min  and  were  counted  for  60  sec  with  the  germanium  detector. 

The  same  procedure  used  in  the  phoconudesr  ratios  experiment  for  estl- 
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macing  ths  activity  was  followed.  The  decay  factor  was  a  lied 
to  each  count  to  calculate  the  activity  of  each  sample  at  the  time 
the  irradiation  stopped.  The  activity  of  each  sample  was  divided  by  the 
sample  weight  to  yield  a  "specific  activity"  for  each.  This  activity  W83 
then  plotted  as  a  function  of  the  sample  relative  position  during  irradi¬ 
ation  and  a  Gaussian  curve  was  fitted  to  the  plot  using  a  non-linear 

regression  technique  (71).  Thus,  a  distribution  of  the  intensity  of 

12  11 

photons  possessing  energies  between  18.7  MeV  [  C(y,n)  C  reaction  thresh¬ 
old]  and  25  MeV  as  a  function  of  distance  from  central  axis  (and  hence 
angle  from  the  Initial  electron  direction)  was  obtained. 

A  weighted  effective  angle  of  multiple-scattering  was  calculated 
from  the  bremsst rahlung  computer  program  for  the  target  slabs  contributing 
photons  of  energy  greater  than  18.7  MeV.  The  weights  were  obtained  from 
the  electron  energy  in  the  slab  [Eqs.  (80)]  and  from  the  relative 
yield  as  a  function  of  energy  measured  in  the  photonudear  ratio  experi¬ 
ment.  This  weighted  effective  angle  was  then  compared  to  the  standard 
deviation  of  the  Gaussian  distribution  obtained  from  the  regression  fit 
and  a  normalization  factor  F  was  calculated 


eff 


me  as 


(85) 


where  is  the  weighted  effective  angle  of  multiple-scattering  (1/e 

width  of  the  distribution)  and  is  the  standard  deviation  of  the 

me  as 

Gaussian  fit.  The  (B  -  X)  factor  of  Eq.  (81a)  was  then  modified  in  a 
manner  similar  to  that  of  Hanson  (72). 


The  ^Al(y,2p)^Na  Cross  Section 

The  aluminum  samples  described  in  the  Materials  section,  of  this 

27  25 

chapter  were  used  to  obtain  the  A1(y>2p)  Na  yields  as  a  function  of 

bremsstrahlung  endpoint  energy.  Yields  were  measured  at  endpoint  energies 

from  25  to  33  MeV  in  0.25  MeV  increments.  Thus,  a  yield  curve  consisted 

of  33  yield  points.  A  total  of  five  yield  curves  was  obtained  and  mean 

yields  were  computed.  This  section  describes  the  experimental  procedure 
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followed  to  obtain  the  Al(y,2p)  Na  yield  curve  and  the  method  of  deter¬ 
mining  the  reaction  cross  section. 

Prior  to  the  start  of  this  experiment,  the  resistance  value  of 
the  variable  resistor  of  the  RC  circuit  of  the  radiation  monitor  system 

was  set  such  that  the  time  constant  (RC)  was  inversely  proportional  to 
25 

the  Na  decay  constant.  For  a  capacitance  of  10.096  yF  and  a  decay  con- 

** 

stant  of  0.01155  sec~\  the  required  resistance  was  8.57  MO.  The  behavior 
of  the  RC  system  was  checked  initially  and  a^  the  beginning  of  each  exper¬ 
imental  session.  The  capacitor  was  charged  to  a  certain  value  and  the 
decay  of  charge  was  observed  as  a  function  of  time  to  determine  whether 
the  potential  was  decreasing  with  rate  constant  RC  ■  1/X.  The  potential 
drop  half-time  was  measured  with  a  stopwatch.  After  the  initial  setting 
of  the  resistance  value,  no  further  adjustments  were  required  as  the 
measured  potential  drop  half-times  were  consistently  reproducible  to 
within  IX. 

Before  the  irradiation  of  the  aluminum  samples,  the  capacitor  of 
the  radiation  monitor  system  was  shorted  to  remove  any  potential.  The 
samples  were  irradiated  in  the  irradiation  geometry  shown  in  Figure  15. 

The  bremsstrahlung  endpoint  energy  was  chosen  at  random  and  the  approp- 
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riate  Irradiation  parameters  (Table  3)  were  set  on  the  linac.  Sample 
irradiation  time  was  2  min.  Immediately  following  irradiation,  the  vol¬ 
tage  reading  on  the  electrometer  was  noted  and  recorded  and  the  sample 
was  retrieved  from  the  irradiation  room.  After  a  dalay  of  55  sec,  the 
samples  wera  counted  for  195  sec  with  the  germanium  detector  spectrometer 
system.  The  counts  in  the  photopaaks  corresponding  to  the  0.975  and 
1.612  MeV  gamma  rays  were  computed  as  described  in  the  Photonuclear  Ratios 
section  of  this  chapter.  Another  endpoint  energy  was  chosen  and  the  pro¬ 
cedure  was  rapeated  until  a  complete  yield  curve  was  obtained.  Five 
yield  curves  wera  measured  in  this  manner. 

At  each  endpoint  anergy  in  a  yield  curve,  a  normalized  yield 
(number  of  reactions  per  volt)  was  calculated  using  Eq.  (53).  The  con¬ 
stant  K  of  Eq.  (53)  was  obtained  from  the  calibration  of  the  germanium 
detector,  from  the  decay  characteristics  of  ^Na,  and  from  an  e“^r  (r  " 
sample  radius)  estimate  of  sample  self -absorption.  C  Is  the  net  number 
of  photopeak  counts  due  to  a  particular  gamma  ray.  A  mean  Yq(T)  was  cal¬ 
culated  based  on  tha  two  estimates  of  yield  made  corresponding  to  each  of 
the  gamma  rays  counted. 

A  total  yield  in  the  entire  energy  range  of  each  yield  curve  was 
determined  by  summing  the  individual  normalized  yields 

n 

T  -  £  Y  Of  >  (86) 

T  i-1  "  1 

where  Y^  is  the  total  yield,  Y^T^)  is  the  normalized  yield  at  endpoint 
energy  T^,  and  n  is  tha  number  of  points  on  the  yield  curve.  For  each  of 
the  five  yield  curves,  a  total  yield  Y^  was  obtained.  A  mean  total  yield. 
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7^,  was  then  calculated  from  all  the  yield  curves 


m 


(87) 


where  7^  la  total  yield  of  the  yield  curve  and  m  Is  the  total  number 
of  yield  curves.  A  normalization  factor,  Y^,  was  calculated  for  each 
yield  curve  by  dividing  Its  total  yield  Into  the  mean  total  yield 


7T/YTj 


(88) 


where  7f^  and  7^  are  the  normalization  factor  and  total  yield,  respec¬ 
tively,  of  the  yield  curve.  Each  normalized  yield  point  of  a  given 
yield  curve  was  then  multiplied  by  its  corresponding  normalization  factor 

vv  ■  VW  ' <89) 


where  Y  (T^)  is  the  corrected  normalized  yield  point  at  endpoint  energy 
In  the  yield  curve  j,  7f^  is  the  normalization  factor  of  the  yield 

A 

curve  J,  and  Ynj(Tj)  the  uncorrected  yield  point  at  endpoint  energy 
In  the  yield  curve  j .  Note  that  this  normalization  procedure  preserves 
the  shape  of  the  yield  curve  while  correcting  It  for  any  possible  varia¬ 
tions  in  experimental  conditions  from  one  experiment  session  to  the  next. 
The  final  yield  curve  was  obtained  by  averaging  the  corrected  yields  at 
each  endpoint  energy  over  all  yield  curves 


VV 


m 

i  I 

1  3-1 


vv 


(90) 


where  *n(T^)  Is  the  mean  normalized  yield  at  endpoint  energy  T^. 
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In  order  to  determine  the  cross  section  of  the  Al(y,2p)  Na 

reaction,  reduced  yields  were  calculated  from  the  mean  normalized  yields. 

These  reduced  yields  were  calculated  using  Eqs.  (13)  and  (15).  The  total 

energy  in  the  bremsstrahlung  spectrum  E(T)  was  obtained  from  the  brems- 

strahlung  computer  program.  The  sample  transmission  factor  f^(T)  was 

experimentally  determined  and  the  monitor  response  R(T)  was  obtained  from 
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Reference  56.  The  Al(y,2p)  Na  reaction  cross  section  was  obtained  with 
the  Cook  (34)  Least  Structure  computer  program  from  the  calculated  reduced 
yields  and  from  the  bremsstrahlung  spectra  obtained  with  the  bremsstrah¬ 
lung  computer  program  (Appendix  E). 


CHAPTER  IV 


RESULTS 

Photonuclear  Ratios 

The  production  of  unique  photonuclear  yield  ratios  as  a  function 
of  bremsstrehlung  endpoint  energy  in  a  simple  manner  required  that  the 
yield  of  photonuclear  reactions  induced  in  element  1  of  the  a  air- '  -  remain 
fairly  constant  over  a  certain  energy  range  while  the  yield  of  photonuc¬ 
lear  reactions  induced  in  element  2  of  the  sample  Increase  over  the  seme 
energy  range.  This  requirement  was  fulfilled  by  the  K^SiFg  end  teflon 
sanples  as  shown  in  Figure  17.  The  yield  curves  of  the  figure  have  been 
normalized  so  that  all  could  be  presented  on  the  same  graph.  As  seen  in 

IQ  10 

Figure  17,  the  FCy.n)  F  reaction  yield  remains  fairly  constant  while 

the  39K(Y,n)38K  and  l2C(y,n)11C  reaction  yields  increase.  The  shape  of 

the  yield  curves  and  their  relative  placement  on  the  endpoint  energy  axis 

follows  from  shape  and  location  on  the  photon  energy  axis  of  the  cross 

sections  for  the  photonuclear  reactions  contributing  to  the  ratios.  These 

19  18 

are  shown  in  Figure  18  (40,73-77).  The  F(Y,n)  F  reaction  threshold  is 

10.4  MaV  (78);  the  cross  section  peaks  at  approximately  4  mb  at  about  12 

MeV  (73).  The  3\(y,n)3®K  reaction  threshold  is  13.1  MeV  (78);  the  cross 

12  11 

section  peaks  at  about  10  mb  around  20  MeV  (74-76) .  The  C(y  C 
reaction  threshold  is  18.7  MeV  (78);  the  cross  section  peaks  at  about  9  mb 


YIELD  (relative  units) 
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around  24  MeV  (40,77) . 

Photonudeer  yield  ratios  were  computed  using  Eq.  (70).  Potas¬ 
sium  and  fluorine  counts  used  to  compute  the  KjSiF^  ratios  typically 
ranged  from  about  4,000  to  11,000  counts.  Figure  19  shows  the  decsy  over 

time  of  a  K^SiF^  sample  exposed  to  20  MeV  bremsstrahlung.  The  presence 

38  18 

of  the  shorter-lived  K  radioisotope  and  the  longer-lived  ^  F  radioiso¬ 
tope  is  evident.  Carbon  and  fluorine  counts  used  to  compute  the  teflon 
ratios  typically  ranged  from  about  15,000  to  30,000  counts.  The  decay 
over  time  of  a  teflon  sample  exposed  to  25  MeV  bremsstrahlung  is  shown  in 
Figure  20.  The  presence  of  the  shorter-lived  and  the  longer-lived 

18_  j 

F  is  also  apparent. 

Mean  photonudeer  yield  ratios  were  computed  by  averaging  the 
ratios  obtained  in  experimental  runs  performed  at  the  same  endpoint  energy. 
Table  7  shows  a  breakdown  of  the  number  of  runs  contributing  to  the  mean 
ratio  at  each  endpoint  energy.  Runs  were  repeated  a  greater  number  of 
tines  at  endpoint  ent  *gies  in  the  region  where  the  ratios  were  changing 
most  rapidly.  Fever  runs  were  performed  at  energies  where  the  ratios 
were  slowly  changing. 

The  photonuclear  yield  ratios  obtained  as  a  function  of  brems- 
atrahlung  endpoint  energy  are  presented  in  Tables  8  and  9  and  are  shown 
graphically  in  Figures  21  and  22.  Evident  in  Figures  21  and  22  is  the 
sensitivity  (percent  change  of  ratio  per  MeV)  of  the  photonuclear  yield 
ratios  to  endpoint  energy  change.  The  maximum  sensitivity  of  the  l^SiF^ 
ratios  is  43Z  per  MeV  at  19  MeV.  The  maximum  sensitivity  of  the  teflon 
ratios  is  41Z  per  MeV  at  24  MeV.  Thus,  for  a  given  accelerator  target 
and  amount  of  flit ret ion,  the  photonuclear  yield  ratios  are  quite  good 


staple  decay 


Figure  20.  Teflon  eaaple  decay 
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Table  7 


Experimental 

runs  used  to  compute  photonuclear  ratios 

KjSiFg  Ratios 

Teflon  Ratios 

Energy  (MeV) 

#  Runs 

Energy  (MeV) 

#  Runs 

13 

1 

20 

2 

14 

5 

21 

3 

15 

5 

22 

3 

16 

5 

23 

4 

17 

7 

24 

4 

18 

6 

25 

4 

19 

6 

26 

3 

20 

5 

27 

5 

21 

5 

28 

4 

22 

4 

29 

3 

23 

3 

30 

4 

31 

3 

32 

3 

Table  8 


K_S1F.  Ratios 


Endpoint 

(MeV) 

Photon  II 

t 

Ratio 

13 

0 

14 

.029 

(.019) 

15 

.074 

(.006) 

16 

.093 

(.010) 

17 

.151 

(.016) 

18 

.235 

(.016) 

19 

.306 

(.029) 

20 

.438 

(.013) 

21 

.515 

(.016) 

22 

.641 

(.032) 

23 

.641 

(.016) 

■Ratios  are  mean  values,  numbers  In  parentheses  are  standard 


errors  of  the  means. 


*Ratios  are 


errors  cf  the  means. 


' 

i 
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Table  9 


Teflon  Ratios  * 


EndpoiQt 

(MeV) 

Photon  II 

Ratio 

20 

0 

21 

.002  (.002) 

22 

.032  (.004) 

23 

.093  (.006) 

24 

.172  (.006) 

25 

.243  (.004) 

26 

.302  (.002) 

27 

.354  (.006) 

28 

.393  (.008) 

29 

.460  (.004) 

30 

.484  (.006) 

31 

.514  (.008) 

32 

.516  (.010) 

tan  values,  numbers  in  parentheses  are  standard 


.7 
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Figure  21.  K^SIF^  photonuclear  ratios. 
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Indicators  of  bremsst rah lung  endpoint  and  hence  accelerator  energy. 

The  photonuclear  yield  ratios  were  found  to  be  sensitive  to  the 
amount  of  filtration  In  the  bremsstrahlung  beam.  When  the  lead  flattening 
filter  (2.8  cm  thick  at  center)  was  Inserted  In  the  beam,  the  ratios  were 
found  to  decrease  as  shown  In  Table  10.  As  seen  In  this  table,  the  per¬ 
cent  reduction  Increases  with  Increasing  endpoint  energy.  This  phenomenon 
has  to  do  with  the  relative  proportion  of  higher-  and  lower-energy  photons 
In  the  bremsstrahlung  spectrum.  It  Is  a  well-known  fact  that  lesd  filtra¬ 
tion  softens  high-energy  x-ray  beams  (14).  The  Increase  in  the  relative 

18 

number  of  lower-energy  photons  Increases  the  relative  F  yield.  The 

38 

decreased  number  of  higher-energy  photons  decreases  the  relative  K  and 
yields.  These  result  In  decreased  rstlos.  In  the  lower  portion  of 
the  energy  range  of  the  K^SIF^  and  teflon  ratio  curves,  ratios  are  prac- 
tlcally  unaffected  by  the  addition  of  the  flattening  filter.  This  may 
follow  from  the  fact  that  In  the  lower  energy  range  there  Is  very  little 
contribution  from  the  relative  and  yields.  This  would  seem  to 
Indicate  that  the  decrease  of  rctlos  with  Increasing  flltretlon  Is  due 
more  to  the  decrease  of  higher-energy  photon  availability  than  to  the 
relative  Increase  of  lower -energy  photons.  This  can  be  substantleted 
since  the  fluorine  threshold  Is  well  above  the  energy  for  the  minimum  y 
for  Pb. 

The  Influence  of  sample  size  on  the  teflon  photonuclear  retlo  was 
also  Investigated.  Samples  of  teflon  rod  0.S,  0.75,  and  1.0  inch  In  dle- 
meter  were  exposed  to  25  MeV  bremsstrahlung  and  ratios  were  computed.  The 
resultant  retlos  were  not  statistically  different  (standard  deviations 
overlapped)  although  the  larger  (1.0-lnch-diameter)  sample  exhibited  s 


Table  10 


Effect  of  Filtration 
on  Photonuclear  Ratios 


Endpoint 

Energy 

Reduction  of 
K^SiFg  ratio 

Reduction  of 

teflon  ratio 

17 

2Z 

19 

5.6Z 

21 

9.6Z 

23 

10. 7Z 

25 

14.  OZ 

27 

15.3Z 

I  29 

15. 4  Z 
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slightly  lower  (by  4Z)  ratio.  It  Is  very  possible  that  larger  sample 
sizes  may  In  fact  result  In  different  ratios.  Thick-target  bremsstrahlung 
beams  contain  relatively  fewer  high-energy  photons  at  larger  angles  from 
the  Incident  electron  direction  (51) .  The  limited  increased  angle  sampled 
in  this  experiment  by  the  larger  teflon  rod  was  not  sufficient  to  estab¬ 
lish  this  possible  cause,  however. 

The  presence  of  possible  Interference  peaks  in  the  vicinity  of 
the  0.511  MeV  photopeaks  used  for  ratio  determination  was  not  observed. 
Thus,  photonuclear  ratios  could  be  measured  with  confidence  using  a  Nal(Tl) 
detector.  Possible  interference  from  the  B*  decay  of  ®^Cu  and  ^Cu  which 

may  be  produced  in  the  copper  present  in  the  aluminum  can  of  the  K^SiFg 

62 

samples  was  not  noted.  The  1.173  and  1.346  MeV  gamma  rays  of  Cu  and 
64 

Cu,  respectively,  could  not  be  seen  in  the  spectra. 


The  Bremsstrahlung  Spectrum 

The  accuracy  of  the  thick-target  bremsstrahlung  calculation  in 
predicting  spectral  shape  is  best  illustrated  by  comparing  the  results  of 
the  calculation  to  experimentally-measured  spectra.  Unfortunately,  very 
few  experimental  determinations  of  bremsstrahlung  spectra  produced  by 
electrons  of  energies  from  about  20  to  35  MeV  have  been  reported  in  the 
llte7/ature.  The  calculations  performed  in  this  work  were  compared  to  the 
21  MeV  thick-target  spectrum  measured  by  O'Dell  (79).  The  present  calcu¬ 
lations  also  are  compared  to  other  calculations  (51,53)  which  attempted 


to  simulate  the  O'Dell  spectrum.  O'Dell  measured  the  spectra  produced  in 

2  2 

a  composite  0.490  g/cm  tungsten  -  0.245  g/cm  gold  radiator  by  electrons 


with  energies  between  5  and  21  MeV.  His  experimental  method  consisted  of 
measuring  the  time  of  flight  of  neutrons  produced  by  the  ^U(y,n)^H  reaction 
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in  e  heevy-water  target  (79).  The  highest-energy  spectrum  was  chosen  for 
comparison  purposes  because  it  was  closest  to  the  bremsstrahlung  endpoint 
energies  of  interest  in  this  work. 

Figure  23  comperes  the  result  of  the  monoenergetic  endpoint  cal¬ 
culations  of  this  work,  indicated  by  the  solid  line,  to  thet  of  Berger 
and  Seltzer  (51),  indicated  by  the  squares.  The  agreement  is  quite  good 
although  the  computational  procedure  of  Berger  and  Seltzer  is  quite  dif¬ 
ferent  from  the  calculation  performed  here.  Berger  and  Seltzer's  compu¬ 
tation  was  besically  e  Monte  Cerlo  simulation.  They  used  the  intrinsic 
bremsstrahlung  cross  section  of  Olsen  and  Maxlmon  et  small  values  of 
screening  parameter  and  the  cross  section  of  Seuter  et  lerge  values  of 
screening  parameter  (51;  see  also  46).  The  computation  of  this  work  was 
e  numerical  analysis  calculation;  the  Schiff  cross  section,  differential 
in  photon  emission  angle  and  energy  (20) ,  was  used  regardleee  of 'the  value 
of  screening  parameter  (the  Schiff  cross  section  assumes  full  screening). 

Figure  24  comparee  the  monoenergetic  endpoint  celculation  per¬ 
formed  here  (solid  line)  to  the  calculation  of  Ferdinands  (53) ,  ehown  by 
the  circles.  Agreement  is  good  except  at  the  epsctrum  tip.  Ferdinande's 
intrinsic  cross  section  was  the  extreme  relativistic  Bethe-Heitler  croes 
section  integrated  over  photon  emission  angle.  He  approximated  the  angle 
of  photon  emission  with  e  Gaussian  fit  to  the  calculation  of  bremsstrah¬ 
lung  mean  square  angles  performed  by  Steeras  (£0) .  The  increased  number 
of  photons  in  the  tip  of  the  Ferdinande  spectrum  may  be  due  to  his  choice 
of  intrinsic  cross  section.  The  extreme  relativietic  Bethe-Heitler  croee 
eection,  with  screening  and  Coulomb  corrections,  contains  a  lerger  number 
of  photons  in  the  tip  than  does  the  Schiff  cross  section  (81).  In  the 
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Figure  24.  Comparison  of  tha  O'Dell  spectrum 
calculated  in  this  work  (solid  line)  with  the  spectrum 
of  Ferdlnande  (S3).  Spectra  have  been  normalized  at 
mid-spectrum. 
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case  of  Ferdinande's  work  it  is  not  desr,  however,  why  an  extreme  rela¬ 
tivistic  cross  saction,  assumed  valid  at  very  large  endpoint  energies, 
would  ba  used  in  the  energy  range  from  5  to  21  MeV.  Also  unclear  is  the 
reason  for  using  an  lntegrated-over-angle  cross  section  and  an  approxi¬ 
mated  angle  of  photon  emission  bssed  on  calculations  performed  for  50  to 
300  MeV  bremsstrahlung  (the  Steams'  calculations). 

The  spectrum  resulting  from  the  calculations  performed  in  this 
work  are,  finally,  compared  to  tha  experimental  data  of  O'Dell  (79)  in 
Figure  25.  The  solid  line  represents  the  monoenergetlc  endpoint  calcula¬ 
tion,  the  dotted  line  represents  the  calculation  with  a  ±3%  energy  spread 
(79),  and  the  triangles  represent  the  experimental  data  of  O'Dell.  Tha 
points  were  obtained  at  integral  MeV  photon  energies  from  the  smooth 
curve  drawn  by  the  authors  through  their  experimental  points.  The  error 
bars  shown  on  selected  points  were  obtained  from  the  same  curve.  Although 
the  present  calculation  appears  to  underestimate  the  mean  number  of  pho¬ 
tons  in  the  tip,  the  calculated  spectrum  is  within  the  error  bsrs  of  tha 
experimental  points.  The  ±32  correction  shows  a  marked  improvement  over 
the  uncorrected  calculation.  It  ahould  be  noted  that  Ferdinande's  calcu¬ 
lation  comes  closest  to  estimating  the  shape  of  O'Dell's  experimental 
spectrum.  However,  recall  thst  Ferdinands  assumes  an  extreme  relativistic 
cross  section  containing  more  photons  in  the  tip  than  does  the  Schiff 
croaa  section,  used  here. 

Figure  26  shows  the  shapes  of  the  flattened  (photon  I)  and  unflat- 
ter.ed  (photon  II)  25  MeV  bremsstrahlung  beams  of  the  Sagittalre  linear 
accelerator  as  predicted  by  the  thick-target  bremsstrahlung  calculations 
of  this  work.  The  effect  of  the  lead  flattening  filter  (photon  I)  is 
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clearly  shown.  Tht  flattened  beam  has  a  lower  relative  number  of  high- 
energy  photons  and  a  higher  relative  number  of  lower-energy  photons  than 
does  the  unflattened  beam.  It  is  this  "beam  softening"  effect  that  is 
the  suggested  cause  of  the  reduction  of  photonuclear  ratios  discussed  in 
the  previous  section. 

Appendix  F  contains  selected  calculations  of  thick-target  brems- 
strahlung  spectra.  The  endpoint  energies  for  which  the  calculated  photon 
and  energy  fluences  are  shown  correspond  to  the  nominal  energies  of  com¬ 
mercially  available  high-energy  accelerators.  These  energies  are  25,  24, 
22,  20,  18,  and  15  MeV.  For  each  energy,  both  a  "monoenergetic  endpoint 
spectrum"  and  a  spectrum  produced  by  electrons  having  a  ±2Z  energy  spread 
are  shown. 


Electron  Multiple-Scattering 

The  calculation  of  thlck-targct  bremsstrahlung  spectra  required 
prior  knowledge  of  the  angular  distribution  of  electrons  scattered  by 
multiple  collisions  in  the  slabs  of  the  target.  A  frequently  used  approxi¬ 
mation  to  the  electron  multiple-scattering  distribution  is  obtained  from 
the  theory  of  Mollfere  (55).  The  first  term  of  Molifere's  formulation  Is  a 
Gaussian  function  for  which  the  RMS  angle  of  multiple-scattering  is  given 
by  Eq.  (81a).  Hanson  (72)  measured  the  angular  distzr.L  tion  of  t'ectrons 
scattered  in  thin  foils  and  found  that  a  renormalization  of  Molifere's 
first  term  resulting  in  a  Gaussian  with  reduced  width  represented  a  better 
approximation  to  his  experimental  data.  Hanson's  renormalization  consisted 
of  subtracting  the  number  1.2  from  the  Moll&re  parameter  "B"  (see  Eq. 
(81a)].  Similar  renormalizations  have  been  performed  (82).  The  results 
of  the  meas<  rement  of  the  angular  distribution  of  electron  multiple  scat- 
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tering  in  the  target  of  the  Sagittaire  accelerator  demonstrated  that 
renormalization  of  Molifcre’s  first  term  into  a  Gaussian  with  reduced  width 
also  was  required  in  this  experiment. 

The  results  of  the  investigation  into  the  multiple-scattering 

distribution  of  electrons  in  the  Sagittaire  target  are  shown  in  Figure  27. 

The  open  triangles  represent  relative  number  of  counts  as  a  function 

of  angle  from  the  bremsstrahlung  beam  central  axis.  The  smooth  curve  is 

the  Gaussian  fit  to  the  experimental  points.  The  standard  deviation  of 

the  Gaussian  curve  [©■  of  Eq.  (85)]  was  4.27  degrees.  The  effective 

angle  of  multiple -scattering  obtained  from  the  bremsstrahlung  computer 

program  of  Eq.  (85)]  was  7.57  degrees  resulting  in  a  normalization 

factor  F  of  1.773. 
n 

The  Moli&re  B  is  a  function  of  slab  thickness.  For  a  slab  thlck- 
2 

ness  of  .0386  g/cm  (slab  thickness  used  in  the  computer  program) ,  B  - 
7.35  (67).  The  RMS  angle  of  multiple-scattering  is  a  function  of  the 
square  root  of  B.  Because  the  normalization  factor  was  obtained  from  a 
measurement  of  angle,  it  was  pquared  first  and  then  divided  into  B  ■  7.35 
to  yield  a  corrected  (renormalized)  B  of  2.34.  Thus  B  was  decreased  by 
subtracting  5.01  from  the  tabulated  value  (67).  It  should  be  pointed  out 
that,  although  this  renormalization  amounts  to  a  substantial  correction 
to  the  Moliire  B  factor,  the  renormalization  procedure  was  found  to  have 
a  practically  negligible  effect  on  the  shape  of  the  bremsstrahlung  spectra 
produced  in  the  computer  program.  The  correction  was  still  incorporated, 
however,  for  the  sake  of  simulating  actual  experimental  conditions  as 
closely  as  possible. 


Figure  27.  Gaussian  fit  to  the  experimental  points  of  the 
measurement  of  electron  multiple  scattering  xngle. 
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27  25 

The  Al(y,2p)  Na  Cross  Section 

27  25 

The  measured  normalized  yield  of  Al(y,2p)  Na  reactions  as  a 
function  of  bremsstrahlung  endpoint  energy  from  25  to  33  MeV  is  shown  in 
Figure  28.  The  units  on  the  ordinate  are  number  of  reactions  per  unit 
response  of  the  radiation  monitor  system  in  volts.  The  standard  errors 
of  the  data  points  are  less  than  the  size  of  the  points.  The  data  below 
27  MeV  show  significant  scatter  because  of  the  absence  of  Na  counts  in 
that  region  resulting  in  poor  counting  statistics.  Above  27  MeV  the  yield 
Increases  continuously.  Several  ''breaks,"  changes  in  the  rate  of  climb 
of  yield  as  a  function  of  increasing  endpoint  energy,  are  apparent  in  the 
yield  curve  between  27  and  about  31  MeV.  These  breaks  occur  at  approxi¬ 
mately  27.3  and  30.3  MeV.  Some  of  these  breaks  have  been  observed  previ¬ 
ously  (7,57),  but  a  determination  of  cross  section  has  never  been  made  to 
verify  if  in  fact  the  breaks  corresponded  to  resonances  in  the  cross  sec¬ 
tion. 

27  25 

In  order  to  determine  the  cross  section  of  the  Al(y,2p)  Na 

reaction,  the  normalized  yield  data  of  Figure  28  were  converted  to  reduced 

2 

yield  (number  of  reactions  per  target  nucleus  per  cm  of  sample) .  This 
vaa  accomplished  using  Eq.  (15).  The  factor  E(T)  of  Eq.  (13)  is  the  total 
energy  in  the  bremsstrahlung  spectrum  as  a  function  of  endpoint  energy. 
Becsuse,  for  comparison  purposes,  the  reaction  cross  section  was  obtained 
using  both  Schiff  integrated-over-angle  and  thick-target  bremsstrahlung 
spectra,  both  spectra  were  Integrated  as  s  function  of  endpoint  energy 
and  simple  analytical  expressions  were  fit  to  the  integration  results. 

The  total  energy  curves  for  Schiff  and  thick-target  spectra  as  a  function 
of  endpoint  energy  between  15  and  35  MeV  are  shown  in  Figure  29.  The 


Normalized  yield  of 
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curves  were  normalized  at  25  MeV.  As  expected,  because  of  the  greater 
number  of  lower-energy  photons  in  the  thick-target  spectra  and  their  rela¬ 
tively  large  contribution  to  the  total  spectral  energy,  the  thick-target 
total  energy  curve  Increases  more  rapidly  with  Increasing  endpoint  energy 
than  does  the  Schlff  total  energy  curve.  The  reduced  yield  curve  obtained 
using  thick-target  total  energy  thus  increased  more  rapidly  with  increas¬ 
ing  endpoint  energy  than  did  the  reduced  yield  curve  obtained  with  Schiff 
spectral  total  energy.  This  is  shown  in  Figure  30  where  the  Schiff 
reduced  yield  is  represented  by  circles  and  the  thick-target  reduced 
yield  is  represented  by  triangles.  The  reduced  yields  of  Figure  30  were 

normalized  at  29  MeV.  Only  data  points  at  half-MeV  intervals  are  shown. 

27  25 

The  Al(y,2p)  Na  cross  section  as  a  function  of  endpoint  energy 

obtained  with  the  Lesst  Structure  computer  program  (36)  using  thick-target 

** 

bremsstrahlung  spectra  corrected  for  the  tlX  endpoint  energy  spread  is 

2 

shown  in  Figure  31.  The  x  (amount  of  smoothing)  applied  to  the  cross 

section  from  27  to  33  MeV  was  99  when  the  Lagranglan  multiplier  (X)  was 
—8 

(1.0)10  (see  the  Analysis  ar.d  Solution  of  Photonuclear  Yields  section 
of  Chapter  II) .  Evident  from  Figure  31  is  the  structure  in  the  cross  sec¬ 
tion  expected  from  the  breaks  seen  in  the  normalized  yield  curve  (Fig. 

28).  Figure  31  shows  cross  section  resonances  st  about  27.3  and  30.3  MeV. 
Relative  resonance  magnitudes  are  directly  related  to  relative  changes  in 
slope  of  the  yield  curve.  The  undershoot  below  27.0  MeV  is  due  most  prob¬ 
ably  to  statistical  fluctuations  produced  by  the  scatter  of  measured 
yields  in  that  energy  range. 

The  absolute  magnitude  of  the  cross  section  was  obtained  by  nor¬ 
malizing  the  Integrated  cross  section  to  thst  of  Aull  and  Whitehead  (21). 


Figure  30.  Schiff  (o)  and  thick  target  (a)  reduced 
yields  (normalized  at  29  MeV) . 


The  Integrated  cross  section  reported  by  the  above  authors  agreed  with 
previous  measurements  (83).  After  normalization,  the  Integrated  cross 
section  from  25  to  32  MeV  was  1.6  MeV-mb,  and  the  peak  cross  section  was 
.35  mb  at  27.3  MeV.  The  magnitude  of  the  cross  section  resonance  at  30.3 
MeV  was  .31  mb. 

27  25 

There  are  differences  between  the  Al(y,2p)  Ha  cross  section 

obtained  here  using  thick-target  bremsstrahlung  spectra  and  the  published 
27  25 

Al(y,2p)  Ha  cross  section  of  Anil  and  Whitehead  (21),  which  is  shown 
in  Figure  32.  One  difference  is  the  lack  of  structure  in  the  Aull  and 
Whitehead  cross  section.  The  above  authors  made  measurements  of  yield  in 
1.0  MeV  intervals  from  25  to  36  MeV  and  in  2.0  MeV  intervals  from  36  to 
65  MeV.  In  addition,  they  smoothed  the  yield  curve  before  analyzing  it 
for  cross  section  determination.  Finally,  the  cross  section  itself  was 
smoothed,  as  illustrated  in  Figure  32  by  the  difference  between  the  smooth 
curve  and  the  points.  The  smoothing  procedures,  in  combination  with  the 
lack  of  resolution  resulting  from  yield  measurements  at  only  integral  MeV 
values  from  25  to  36  MeV,  may  account  for  the  lack  of  structure  in  that 
energy  range.  In  the  present  work,  yield  measurements  were  made  in  0.25 
MeV  increments  and  no  smoothing  was  applied  to  the  yield  curve.  The 
increased  resolution  of  the  present  measurements  and  the  use  of  Least 
Structure  to  analyze  the  data  (requiring  no  yield  smoothing)  permitted  a 
more  detailed  examination  of  possible  structure  xn  the  cross  section. 

Other  differences  between  the  cross  section  reported  here  and 
that  of  Aull  and  Whitehead  are  the  location  of  the  peak  cross  section  and 
the  cross  section  trend  with  increasing  endpoint  energy.  The  above 
authors*  cross  section  increased  constantly  with  increasing  endpoint 


Figu-e  32.  The  27Al(r,2p)25Na  cross 
section  reported  by  Aull  and  Whitehead. 
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energy  end  peaked  at  32.0  MeV  (21).  Tne  cross  section  reported  here  peaked 
at  about  27.3  MeV,  contained  a  resonance  at  30.3  MeV,  and  thereon  had  an 
essentially  downward  trend  with  Increasing  endpoint  energy.  It  is  suspec¬ 
ted  that  the  difference  in  peak  cross  section  location  is  due  to  the  dif¬ 
ferent  bremsstrahlung  spectra  used  for  cross  section  determination.  Aull 
and  Whitehead  analyzed  their  data  using  Penfold  and  Lelss  matrices  (34) 

which  were  obtained  from  thin-target  (Schiff)  brerasatrahlung  apectra  even 

2 

though  they  report  a  target  thickness  of  .98  g/cm  tungsten.  The  croas 
section  reported  hare  was  obtained  using  thick-target  bremsatrahlung 

27  25 

spectra.  The  Al(y,2p)^  Na  yield  curve  measured  in  this  work  was  also 
analyzed  using  Schiff  spectra  to  demonstrate  tha  dramatic  effect  of  brems¬ 
atrahlung  spectral  shape  on  tha  resultant  crosa  section. 

The  cross  section  as  a  function  of  endpoint  energy  obtained  by 
27  25 

analyzing  the  Al(y,2p)  Na  yield  curve  measured  in  this  study  utilizing 

2 

Schiff  spectra  is  shown  in  Figure  33.  The  x  applied  to  the  cross  section 

-12 

from  27  to  33  MeV  was  58  when  X  “  (1.0)10  .  The  cross  section  now  peaked 

around  32.5  MeV.  This  result  is  in  agreement  with  that  of  Aull  and  White- 
head  (21),  who  obtained  a  peak  cross  section  at  32.0  MeV.  The  differences 
in  the  two  cross  sections  obtained  from  the  same  yield  data  uaing  both 
thin-  and  thick-target  bremsstrahlung  spectra  can  be  due  only  to  differ¬ 
ences  in  thin-  and  thick-target  spectral  shapes. 

Figure  34  shows  the  Schiff  and  thick-target  30  MeV  spectra  used 

27  25 

in  the  Least  Structure  program  to  produce  the  Al(y,2p)  Na  thin-  (Schiff) 
and  thick -target  cross  sections.  The  curves  have  been  normalized  at  27.375 
MeV.  Note  the  much  largar  proportion  of  lower-energy  photons  in  the  thick- 
target  spectrum  as  compared  to  the  relatively  constant  number  of  low-energy 
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Figure  34.  30  MeV  Schiff  and  chick- target  spectre 

used  in  the  cross  section  analyses. 


140 


photons  in  the  Schiff  spectrum.  In  an  analysis  of  photonuclear  yield  to 
produce  a  cross  section,  this  characteristic  shape  of  the  thick-target 
spectrum  will  have  a  profound  effect  on  the  resultant  cross  section.  At 
a  given  endpoint  energy,  the  numbar  of  photonuclear  reactions  induced  in 
a  a ample  are  produced  by  all  photons  in  the  spectrum  of  energy  exceeding 
the  threshold  energy  of  the  reaction.  In  the  case  of  thick-target  brems- 
atrahlung,  because  of  the  large  proportion  of  lower-erergy  photons ,  more 
reactions  are  produced  by  photons  in  the  lower  energy  range  than  are  pro¬ 
duced  by  photons  near  the  tip  of  the  spectrum  (assuming  tha  cross  section 
does  uot  change  appreciably  over  the  same  energy  range).  This  fact  was 

verified  by  examining  the  resolution  function  of  the  thick-target 
27  25 

Al(y,2p)  Na  crosa  section  and  noting  the  influence  of  lower  energies 
on  crosa  sections  at  higher  energies.  Thus,  for  a  given  photonuclear 
reaction,  yields  produced  by  thick-target  bremsstrahlung  increase  more 
rapidly  with  increasing  endpoint  energy  than  will  yields  produced  by  thin- 
target  (Schiff)  bremsstrahlung.  If  a  yield  curve  produced  with  thick- 
target  bremsstrahlung  is  analyzed  for  cross  section  using  Schiff  spectra, 
the  resultant  cross  section  will  show  a  continuous  increase  with  increas¬ 
ing  endpoint  energy.  Thia  ia  due  to  the  combined  effect  of:  1)  an 
Increased  availability  of  photons  possessing  energies  that  exceed  the 
threshold  energy  of  the  reaction  thus  resulting  in  a  more  rapidly  increas¬ 
ing  yield,  and  2)  the  lack  of  a  proportionately  larger  number  of  low- 
energy  photons  in  the  bremsstrahlung  spectrum  used  to  analyze  the  yield. 

The  aforementioned  effect  ia  responsible  for  the  cross  section 
shown  in  Figure  33.  The  cross  section  continuously  increases  because  of 
the  increased  yields  produced  by  lower-energy  photons  in  the  thick -target 
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bremsstrahlung  used  to  produce  the  yields.  The  good  agreement  between 

the  thin-target  cross  section  obtained  here  and  the  Aull  and  Whitehead  • 

cross  section  seems  to  suggest  that  the  above  authors  may  have  made  an 

unfortunate  choice  of  bremsstrahlung  spectrum  in  their  analysis.  Aull 

and  Whitehead  irradiated  their  aluminum  samples  with  bremsstrahlung  pro- 

2 

duced  in  a  0.98  g/cm  tungsten  target  (21)  and  analyzed  their  yield  curve 

using  Schiff  spectra.  The  choice  seems  inappropriate,  particularly  in 

view  of  the  relatively  good  agreement  of  the  thick-target  calculations 

performed  In  this  work  with  the  data  of  O’Dell  (79)  which  was  obtained 

2 

with  a  target  of  less  thickness  (0.735  g/cm  ) . 


CHAPTER  V 


CONCLUSIONS 

Photonuclear  Ratios 

Photonuclear  reactions  are  particularly  well-suited  for  the  study 
of  high-energy  bremsst rah lung  spectra.  Ihelr  aptness  as  a  research  tool 
for  this  type  Investigation  Is  afforded  by  virtue  of  tha  dependance  of 
photonuclear-reactlon  cross  sections  on  photon  energy.  Photonudear-reac- 
tlon  thresholds  and  cross-eectlon  resonances  have  been  used  In  the  energy 
calibration  of  linear  accelerators  (6-8)  and  in  the  determination  of  thin- 
target  x-ray  spectra  (5) .  Rati  os  of  photonuclear  ylalds  In  salect  mater¬ 
ials  have  been  used  as  Indices  of  thin-target  x-rsy  quality  (9,10).  The 
research  reported  bare  utilized  photonuclear  reactions  to  study  the  dis¬ 
tribution  In  energy  of  photons  In  thick-target  x-ray  spectra  and  the  effact 
of  spectral  shape  on  photonuclear  yield  measurements  and  cross  saction 
determinations. 

Photonuclear  yield  ratios  In  K.S1F.  and  In  teflon  hava  been 

Z  D 

obtained  which  serve  as  indices  of  thlck-targat  bremsstrahlung  quality  In 
the  ensrgy  ranga  from  about  14  to  30  MeV.  Many  commercially  available 
linear  accelerators  possess  nominal  energies  in  this  range,  with  endpoint 
energies  of  15,  18,  20,  22,  24,  and  25  MeV  being  most  common.  Photonuclear 
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ratios  published  previously  (9,10)  were  obtained  from  essentially  thin- 
target  x  rays  and  failed  to  provide  unique  indices  from  15  to  25  MeV. 

Ratios  in  I^SiF^  demonstrated  a  sensitivity  (percent  change  of  ratio  per 
MeV)  of  about  262  at  15,  312  at  18,  and  172  at  20  MeV.  Ratios  in  teflon 
had  sensitivities  of  about  75Z  at  22,  412  at  24,  and  272  at  25  MeV. 

The  ratios  obtained  here,  and  hence  the  reported  sensitivities 
as  a  function  of  endpoint  energy  changs,  apply  in  a  strict  sense  only  to 
lightly-filtered  (2.7  mm  Pb)  thick-target  (4  mm  W)  bremsstrahlung  produced 
by  electrons  posssssing  an  energy  spread  of  ±22  of  the  mean  kinetic  ensrgy. 
The  insertion  of  a  thicker  flattening  filter  (2.8  cm  Pb),  on  the  other 
hand,  reduced  the  ratios  by  only  approximately  152  in  the  worst  case  (29 
MeV  teflon  ratio).  The  effect  of  different  targst  thicknesses  on  ratio 
senaitivity  could  not  be  inveatigated.  Nath  and  Schulz  (9) ,  however, 

-  ce  gg  ^ 

report  only  a  102  decreaae  of  their  Cu/  Y  ratio  at  30  MeV  when  the 
target  thickness  was  increased  from  about  .35  to  2.2  radiation  lengths. 
Different  initial  electron* energy  distributions  also  could  affect  ths 
reported  magnitudes  and  sensitivities  of  the  ratios.  The  result  of  elec¬ 
tron-energy  spread  on  the  bremsstrahlung  spectrum  aeems  to  be  a  net 
increase  in  the  number  of  photons  comprising  the  apectrum  tip.  Thus,  it 
ia  suspected  that  an  increased  electron-energy  sprsad  may  increase  the 
magnitude  of  the  ratio  at  the  mean  endpoint  energy.  On  the  other  hand, 
ratio  senaitivity  should  not  be  appreciably  affected. 

Because  ratios  in  K0SiFg  and  teflon  are  independent  of  the  masa 
of  the  sample,  the  flux  of  photons  and  the  efficiency  of  the  counting 
system,  they  should  be  extremely  useful  to  medical  physicists  seeking  to 
augment  their  accelerator  quality  control  program  with  a  simple  yet  sen- 
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sitive  index  of  x-ray  quality.  The  absence  of  strong  lines  in  the  vici¬ 
nity  of  the  .511  MeV  photopeak  suggests  that  a  Nal(Tl)  counting  system 
could  be  used  in  lieu  of  a  semiconductor  system,  provided  that  background 
corrections  are  made  in  a  careful  and  reproducible  manner. 

The  Bremsstrahlung  Spectrum 

The  review  of  the  photonuclear  literature  which  was  conducted  in 
preparation  for  this  work  conveyed  the  fact  that  thick-target  bremsstrah¬ 
lung  spectra  had  never  before  been  utilized  in  the  analysis  of  photonuclear 
yields  for  cross-section  determination.  To  this  purpose  a  computer  program 
was  developed  which  consisted  of  a  numerical -analysis  c:  ’culation  of  thick- 
target  bremsstrahlung  spectra.  The  computer  program  was  capable  of  calcu¬ 
lating  spectra  at  any  endpoint  energy  from  10  to  35  MeV  from  targets  of 
essentially  any  thickness  and  atomic  number.  Provisions  were  made  to 
Include  the  effects  of  endpoint-energy  spread.  Electron-energy  spread  has 
never  been  Incorporated  before  into  calculations  of  this  type.  The  program 
was  checked  by  comparing  the  results  of  a  simulation  made  of  the  21  MeV 
thick -target  spectrum  of  O'Dell  (79)  with  his  experimental  measurements 
and  with  the  results  of  other  calculations  (51,53).  In  general,  the 
results  of  the  simulation  compared  quite  favorably  with  the  published 
measurements  and  calculations,  particularly  when  electron -energy  spread 
was  considered.  The  largest  discrepancies  occurred  at  the  high-energy  end 
(tip)  of  the  spectrum.  The  reasons  for  the  discrepancies  are:  the  inabil¬ 
ity  of  Born  approximation  calculations  to  characterize  the  spectrum  tip, 
the  unavailability  of  calculations  appropriate  for  spectra  of  energies 
less  than  about  30-50  MeV  that  do  not  incorporate  the  Born  approximation, 
and  the  extreme  experimental  difficulties  of  measuring  precisely  the  tip 
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of  bremsstrahlung  spectra.  Even  to  this  day,  further  theoretical  and 
experimental  research  is  needed  badly  in  this  area. 

The  ^Al(y,2p)^Ka  Cross  Section 
27  25 

The  Al(y,2p)  Na  cross  section  reported  here  was  obtained  from 
measurements  of  yield  made  at  0.25  MeV  increments  from  25  to  33  MeV.  This 
cross  section  had  never  been  calculated  with  such  energy  resolution.  The 
increased  resolution  of  this  experiment  allowed  inspection  of  possible 
structure  in  the  cross  section.  The  cross  section  of  Figure  31  shows  two 
resonances,  one  at  27.3  MeV  and  one  at  30.3  MeV.  This  structure  is  unique, 
possibly  indicating  the  existence  of  states  at  these  high  energies. 
Absorption  by  2  fioj  single-particle -like  transitions  from  deep-lying  shells 
may  be  responsible  for  the  cross  section  at  these  energies  (43).  Transi¬ 
tions  from  lp  and  Is  shell  states  may  be  important  in  the  absorption  of 

27 

photons  above  20  MeV  in  nuclei  such  as  Al. 
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If  the  potential  energy  la  independent  of  time,  we  cen  separate  space  and 
time  veriables  (a  technique  for  solving  certain  dlffe.sntiel  equetions 
which  is  used  frequently)  by  setting 

*  ■  <K*,y,*,)c(t)  (A2) 

where  t(x,y,z)  and  c(t)  are  functions  exclusively  of  space  and  of  time 
variables,  respectively.  Substituting  Eq.  (A2)  into  Eq.  (Al)  and  dividing 
by  yields 

#$**-££•  <«> 

Since  the  left-hand  side  of  this  equation  depends  only  on  space  variables 
and  the  right-hand  side  only  on  time,  the  equetion  cannot  be  satisfied 
for  all  points  in  space  at  ell  times  unless  each  side  is  equal  to  the 
same  constant.  Call  this  “separation  constant"  the  total  energy  of  the 
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system,  E.  The  right-hand  side  of  Eq.  (A3)  is  readily  solvable  yielding 
solutions 

;  -  exp[-i(E/ft)t]  (A4) 

where,  for  convenience,  the  constant  of  integration  has  been  set  equal  to 
unity.  The  left-hand  side  of  Eq.  (A3)  can  be  written  (after  rearranging) 
as 

j.2  2 

+  V<|>  -  E<|>  .  (A5) 

This  equation  is  the  time -independent  Schrodinger  equation,  an  equation 
which  will  be  referred  to  often  in  this  °«ction. 

The  basic  assumption  for  any  shell  model  is  that  despite  the 
strong  overall  interaction  between  nucleons  which  provides  the  binding 
energy  of  the  nucleus,  the  motion  of  each  nucleon  is  practically'  indepen¬ 
dent  of  that  of  any  other  nucleon.  Each  nucleon  is  then  assumed  to  move 
in  and  under  the  influence  of  the  same  potential.  Under  this  assumption, 
the  model  is  known  as  the  single -particle  or  independent -par  tide  shell 
model.  In  the  simplest  case,  the  potential  is  assumed  to  be  spherically 
symmetric  in  nature.  Under  this  condition,  Schrodinger 's  time-indepen¬ 
dent  equation,  Eq.  (A5),  may  be  expressed  in  spherical  coordinates  (Fig. 
Al)  r,  ©,  +,  as  (86,87) 

1  3_(r2  3£)  1  3_(sin©  9*) 

r7  3r  3r  r^sin©  3©-  30- 

2 

+  r2^  W  +  &  (E  “  V(r))*  ‘  0  (A6) 

where  now  ^  »  ip(r,0,$)  and  the  potential  V(r)  is  a  function  of  r  exclu¬ 
sively. 
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Equation  (A6)  may  be  readily  separated  Into  three  Independent 
equations,  each  Involving  only  one  coordinate 

-  R(r)-e«(e)*(*)  (A7) 

where  R,  ift,  and  *  depend  only  on  r,  6>,  and  4  respectively.  Carrying  out 

2 

the  separation  and  defining  separation  constants  ,  L(L  +  1),  where 
and  t  are  to  be  defined  shortly,  yields  three  differential  equations 
each  a  function  of  a  single  variable  (87) 


The  solution  of  Eq.  (A8)  along  with  the  condition  that  4  mu3t  be 
single-valued  at  a  given  point  In  space  results  In  the  condition  that  m^ 
be  either  0  or  a  negative  or  positive  Integer.  Equation  (A9)  has  a 
rather  complicated  solution  In  terms  of  associated  Legendre  functions 
which  exist  only  when  the  constant  t  Is  an  Integer  equal  to  or  greater 
than  the  absolute  value  of  m^.  The  above  requirements  can  be  expressed 
as  conditions  on  m^  such  that 

"^1»  •••  *  0,  —  1»£  (All) 

and  on  t  such  that 


l  -  0,  1,  2,  3,  ... 


(A12) 
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The  constants  m^  and  L  are  known  as  the  magnetic  quantum  number 
and  the  orbital  quantum  number,  respectively.  Quantum  numbers  describe 
the  "state”  of  a  nucleon.  Similarly  to  atomic  physics,  for  t  -  0,  1,  2, 

3,  4,  5  the  spectroscopic  notation  s,  p,  d,  f,  g,  h  is  used.  From  Eq. 
(All)  it  is  evident  that  for  a  given  orbital  quantum  number  t%  the  mag¬ 
netic  quantum  number,  m^,  may  take  on  one  of  the  2£  +  1  values.  These 
values  correspond  to  the  2£  +  1  possible  orientations  of  the  angular 
momentum  vector  of  a  nucleon  when  in  the  presence  of  a  magnetic  field. 

A  p-state  it  -  1)  nucleon,  for  example,  may  have  any  one  of  three  magnetic 
sublevels:  m^  ■  -1,  0,  or  1.  A  third  quantum  number,  commonly  denoted 
by  m  or  s,  has  been  called  the  intrinsic  spin.  It  describes  the  direc- 

3 

tlon  of  the  nucleon's  spin  and  may  take  on  the  values  m  -  \  or  m  ■  -h • 

•  s 

The  radial  wave  equation,  Eq.  (A10) ,  can  be  transformed  into  a 
form  which  is  analogous  to  a  simple  one-dimensional  wave  equation  (86). 

To  accomplish  this,  the  so-called  modified  radial  wave  function,  x(r)» 
defined  by 


x(r)  -  rR(r)  (A13) 

la  substituted  into  Eq.  (A10).  The  radial  wave  equation  then  becomes 
simply  (after  rearranging) 


2m 


2d2 

dP 


-  V(r)  + 


Lit  +  Dtt" 

2m  r* 


EX 


(A14) 


Solution  of  this  equation  requires  knowledge  of  the  exact  shape  of  the 
potential  function  V(r). 

The  radial  shape  of  the  nuclear  field  is  well-known  from  scat¬ 
tering  experiments,  but  its  use  in  Eq.  (A14)  makes  the  solution  of  the 
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Schrodinger  wave  equation  extremely  complicated.  Therefore,  simpler  forms 
of  the  potential,  V(r) ,  commonly  are  used.  The  choice  of  the  potential 
form  is  restricted  by  two  rather  stringent  conditions  (88): 

1.  Since  V(r)  represents  the  actions  of  the  other  nucleons  and 
the  nucleon  density  at  the  point  r  ■  0  is  not  essentially  different  from 
that  at  any  other  place  inside  the  nucleus,  V(r)  at  r  *  0  is  not  expected 
to  be  singular.  Furthermore,  since  V(r)  is  spherically  symmetric. 


dr 


0 


(A15) 


2.  The  potential  must  go  to  zero  rather  abruptly  at  the  nuclear 
surface,  that  is  when  r  approaches  the  nuclear  radius 


for  r  ■*  R  . 


-  <A16) 


There  are  two  well-known  potentials  that  satisfy  those  conditions:  the 
square  well  and  the  harmonic  oscillator  potentials. 

The  square  well  potential  has  the  form 


V(r)  - 


r  <_  R 
r  >  R 


(A17) 


The  shape  of  the  square  well  potential  is  presented  in  Figure  A2a.  A 
solution  to  Eq.  (A14)  when  V(r)  is  given  by  Eq.  (A17)  is  sought;  that  is. 


-ft2  d2 

ara?7 


-  V  + 


IU  +  m 
— 


Ex 


(A18) 


When  t  «  0  (s  states),  the  problem  reduces  to  that  of  a  particle  in  a  box 


•  Figure  A2 .  The  square  well  (a) 

and  harmonic  oscillator  (b)  potentials. 

r-e 
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of  sides  R  and  potential  within  of  -V^.  The  energy  eigenvalues  which 
result  are  given  by 


where  n,  the  radial  quantum  number,  is  the  number  of  nodes  of  Bq.  (A13). 
When  t  >  0,  the  situation  is  more  complicated.  Solutions  are  given  by  a 
set  of  functions  known  as  spherical  Bessel  functions  of  order  t  (89*90) 

■  <-1)f  (li  Cr  hf  TT 

where  the  J^(kr)  are  the  spherical  Bessel  functions  of  order  t  and  the 
constant  k  is  given  by 

k2  -  2nE/tt2  ^  (A21) 

The  energy  levels  resulting  from  Eqs.  (A19-A21)  are  shown  in 
Figure  A3a.  Each  level  is  2(2 l  +  1)  orientationally  degenerate.  The 
numbers  in  parentheses  above  each  energy  level  of  Figure  A3a  correspond 
to  the  accumulated  total  of  nucleons  resulting  in  subshell  closures  up 
to  that  level.  These  are  2,  8,  18,  20,  34,  40,  38,  etc.  and,  with  the 
exception  of  2,  8,  and  20,  do  not  correspond  to  the  experimentally 
obtained  "magic  numbers"  2,  8,  20,  28,  SO,  82,  and  126,  corresponding  to 
the  number  of  protons  or  neutrons  in  the  nuclei  of  particularly  stable 
elements . 

The  harmonic  oscillator  potential  has  the  form 
f -Vq(1  -  (r/R2)]  r  <  R 

V(r)  - 


0 


r  >  R 


(A22)  . 
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The  shape  of  the  harmonic  oscillator  potential  is  shown  In  Figure  A2b. 
Solutions  (91)  of  the  radial  wave  equation  (Eq.  (A14)] ,  when  the  poten 
tial  is  given  by  Eq.  (A22) 


~-r  -  V  (1  -  y  ■  Ey 

STd?  V1  W  +  2inr7  X  EX 


result  in  energy  eigenvalues  given  by  (88) 


Ene  "  ^2(n  “  ^  +  ^Rt0  +  ~  Vc 


(A23) 


(A24) 


where 

u  -  (2V  /mr2)*3 
o 


(A25) 


The  energy  levels  resulting  from  Eqs.  (A24-A25)  are  depicted  in 
-Figure  A3b.  Note  that  in  addition  to  the  2(2£  +  1)  orientational  degen¬ 
eracy,  except  for  the  first  two  levels,  this  level  scheme  has  "accidental" 
degeneracy  between  levels  of  different  n  and  t.  Note  also  that  after  the 
first  two  levels,  each  subsequent  oscillator  shell  begins  with  l  larger 
by  one  unit  and  contains  all  t  values  of  the  same  parity  (a  symmetry 
operator  having  to  do  with  wave  function  evenness/oddness)  down  to  the 
lowest  one.  As  in  Figure  A3a,  the  numbers  in  parentheses  (2,  8,  20,  40, 
70,  112)  represent  the  total  occupation  number  of  all  lower  levels  up  to 
the  one  undar  consideration  (inclusive).  With  the  exception  of  the  first 
three  occupation  numbers  (2,  8,  20),  again  the  magic  numbers  (2,  8,  20, 

28,  50,  82,  126)  are  not  reproduced. 

So  far  the  classification  of  levels  has  been  basad  only  on  the 
orbital  momentum  quantum  number  £,  under  the  assumption  that  spin-orbit 
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Interactions  play  only  an  unimportant  role.  This  Is  so  In  the  atomic 
case,  where  this  Interaction  causes  only  a  alight  splitting  of  the  levels 
j  ■  Z  ± \  which  Is  small  compared  to  the  distance  between  levels  of  dif¬ 
ferent  Z.  In  the  nuclear  realm,  however,  experimental  evidence  has 
revealed  the  presence  of  a  strong  Interaction  between  orbital  angular 
momentum  and  Intrinsic  spin  angular  momentum  which  may  cause  a  rather 
large  splitting  of  energy  levels  that  may  no  longer  be  Insignificant  when 
compared  to  the  distances  between  levels  of  different  Z. 

A  spin-orbit  coupling  model  has  been  auggested  In  which  there  Is 
a  strong  spin-orbit  interaction,  proportional  to  L*S,  in  addition  to  the 
potential  V (r) ,  acting  on  a  nucleon  In  the  nucleus .  The  magnitude  of 
this  spln-orblt  Interaction  depends  on,  in  addition  to  the  magnitude  of 
the  orbital  angular  momentum  vector,  the  relative  orientation  of  the  spin 
and  the  orbital  angular  momentum  vectors.  The  spln-orblt  Interaction  in 
this  case  (In  contrast  to  the  atomic  caae)  Is  Inverted,  which  means  that 
a  nucleon  in  a  j  *  Z  +  %  state  has  less  energy  than  a  nucleon  In  a  j  ■ 

Z  -  H  state.  Since  the  energy  of  a  level  depends  strongly  on  the  align¬ 
ment  of  spin  and  orbit,  It  Is  no  longer  justifiable  to  characterize  a 
level  simply  by  Its  Z  value.  The  total  angular  momentum  j  must  be  speci¬ 
fied  also. 

The  level  sequence,  a  harmonic  oscillator  potential  [Eq.  (A22)] 
c  dified  by  strong  spln-orblt  coupling,  Is  shown  In  Figure  A4  (88). 

Owing  to  the  spln-orblt  coupling,  the  nuclear  shells  are  not  the  same  as 
the  oscillator  shells.  It  Is  seen  that  In  the  three  loweat  oscillator 
levels  groups  (Is,  lp,  2ald)  the  grouping  Is  not  appreciably  changed. 
Hence,  the  lower  shell  occupation  numbers  (2,  8,  20)  are  not  affected  by 


Figure  A4.  Rarmonic  oscillator  energy 
levels  with  spin-orbit  coupling.  Ref.  (88). 
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the  spin-orbit  splitting.  However,  the  level  lf7/2  of  the  next  oscillator 
level  group  is  appreciably  lowered  and  is  brought  Into  Isolation  frost  the 
other  levels  which  indicates  shell  closure  at  occupation  number  28.  Still 
mora  drastically,  the  lg9/2  level  from  the  oscillator  group  is  brought 
down  close  to  the  praceding  oscillator  group  producing  a  wide  gap  after 
occupation  number  50.  Here  the  spin-orbit  splitting  clearly  becomes  tha 
dominant  feature  in  tha  level  arrangement.  Thus,  tha  spin-orbit  coupling 
reproduces  all  magic  numbers  exactly  .is  can  be  evidenced  from  Figure  A4. 

Although  the  spin-orbit  sheil  model  of  nuclear  structura  pro¬ 
vides  a  sound  explanation  for  many  nuclear  phenomena  which  are  experimen¬ 
tally  observed,  the  simple  modal  given  above  is  not  sufficient.  For 
exanple,  tha  model  cannot  explain  why  any  even  number  of  identical  nuc¬ 
leons  couples  to  zero  ground  state  spin.  Evidently  there  must  be  a 
nucleon-nucleon  interaction  which  favors  the  pairing  of  nucleons* with 
opposing  angular  momenta.  An  attractive  interaction  between  nucleons 
must  therefore  be  added  to  the  single-nucleon  spin-orbit  interaction, 
which  gives  rise  to  a  pairing  anergy.  From  theoretical  considerations, 
it  appears  likely  that  the  magnitude  of  tbe  pairing  energy  Increases  with 
the  t  value  of  the  pair  and  that  for  this  reason  tha  hlgb-spln  states 
predicted  by  tha  spin-orbit  model  are  not  found  in  odd-A  ground  state 
spins  (85). 

Another  feature  not  included  in  the  simple  shell  model  is  the 
distorting  effect  of  the  outermost  nucleons  on  the  other  nucleons  in  a 
nucleus.  Suppose  a  single  nucleon  is  added  to  a  closed  shell  nucleus. 

This  nucleon  will  usually  have  a  high  t  value  and  therefore  its  wave 
function  will  peak  close  to  tha  nuclear  radius.  In  terms  of  a  classical 
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picture  the  nucleon  circles  around  a  closed  shell  nuclear  core.  The 
attraction  between  the  nucleon  and  the  core  distorts  the  core.  1£  more 
nucleons  are  added  outside  of  the  core,  a  point  will  be  reached  when  it 
Is  permanently  deformed  with  an  accompanying  effect  on  the  orbits.  The 
deformation  may  result  in  a  nucleus  of  ellipsoidal  shape.  A  quantum- 
mechanical  body  which  has  an  axla  of  symmetry  such  as  an  ellipsoid  of 
revolution  can  undergo  rotations  about  an  axis  perpendicular  to  the  sym¬ 
metry  axis.  The  energy  spectrum  of  such  a  "rotator"  Is  quite  character¬ 
istic  (85) .  There  ire  many  cases  of  even-even  nuclei  In  which  this  spec¬ 
trum  has  been  found.  This  leads  to  the  conclusion  that  deformed  nuclei 
indeed  exist. 

In  any  finite  potential  well  there  are  bound  energy  levels, 
corresponding  to  states  of  energy  E  <  0  [assuming  V(r  -*■•*)  ■  0] ,  and 
unbound  or  virtual  levels  of  states  of  energy  E  >  0.  Virtual  states 
result  from  the  fact  that  the  deBroglle  wave  of  a  nucleon  Is  reflected  at 
the  edge  of  the  potential  well  even  though  Its  total  energy  may  exceed 
the  potential  energy  In  the  well.  Thus,  approximately-standing  waves  are 
formed  within  the  potential  well,  the  amplitude  of  the  wave  function  can 
be  very  large  within  the  well,  and  a  virtual  state  occurs.  A  nucleon 
may  vacate  a  virtual  state  after  aoaa  mean  life  t,  causing  it  to  have  a 
width  for  particle  emission  r  given  by 

r  -  */t  (A26) 

A  virtual  state  may  also  decay  to  a  lower  state  by  gasma-ray  emission. 
Bound  levels,  on  the  other  hand,  may  decay  only  by  gamma  emission. 

According  to  the  single  particle  shell  model,  a  given  nucleus 
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(Z,N)  consists  of  7.  protons  and  N  neutrons  each  pieced  into  bound  levels 
of  increasing  energy  in  accordance  with  the  Pauli  exclusion  principle. 

In  the  ground  stete  of  the  nucleus,  all  nucleous  are  in  their  lowest 
energy  states.  The  simplest  excited  stetes  of  the  nucleus  are  formed  by 
raising  the  outermost  nucleon  to  a  higher  stete.  The  corresponding  exci¬ 
tation  spectrum  of  the  nucleus  is  called  the  single  perticle  level  spec¬ 
trum  with  levels  as  shown  in  Figure  A4.  However,  upon  inspection  of 
actual  level  spectra,  many  more  levels  than  expected  from  the  single- 
particle  spectrum  are  observed.  Present  theories  indicate  that  the  single¬ 
particle  spectrum  is  ectually  "dissolved"  among  many  levels  eech  of  which 
consists  of  complicated  excitations  of  more  then  one  particle  (85) . 
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In  the  quantum  mechanical  treatment  used  to  explain  nuclear 
structure  (Appendix  A),  the  simplifying  assumption  that  the  potential 
energy  term,  V,  of  the  Schrodinger  equation  was  constant  in  time  vas  made. 
This  assumption  was  necessary  to  allow  separation  of  the  complete  equation 
into  time-dependent  and  time-independent  equations.  Use  was  made  of  the 
time-independent  equation,  Eq.  (A5) ,  and  time -invariant  expressions  for 
nuclear  potentials  to  obtain  energy  states  that  were  necessarily  constant 
in  time.  No  explanations  concerning  transitions  of  any  kind  from  one 
state  to  another  were  provided.  In  order  to  describe  transitions  and 
allowed  states,  quantum  mechanical  perturbation  theory  is  employed. 

The  cases  in  which  the  Schrodinger  equation,  Eq.  (Al) ,  can  be 
solved  exactly  are  extremely  rare.  Perturbation  theory  is  the  most  impor¬ 
tant  method  for  obtaining  approximate  solutions  to  Schrodinger *s  equation. 
The  method  is  applicable  to  any  situation  in  which  the  Hamiltonian  H  (see 
below)  describing  the  system  under  consideration  is  not  much  different 
from  the  Hamiltonian  H^,  describing  some  similar  system  simple  enough  so 
that  solutions  to  the  Schrodinger  equation  are  known  exactly  (89) . 

The  Hamiltonian,  or  total  energy  operator,  H  given  by 


H 


+  V 


(Bl) 


is  simply  the  operator  on  the  left-hand  side  of  Eq.  (Al) .  When  a  system 
is  perturbed  (such  as  in  the  emission  or  absorption  of  electromagnetic 
radiation),  its  total  Hamiltonian,  H,  is  affected  and 

H  -  H  +  H’  (32) 

o 


where  Hq  is  that  part  of  the  Hamiltonian  that  permits  calculation  of  the 
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properties  of  the  allowed  states  of  the  system  and  H'  is  the  effect  of 
the  perturbation  on  the  Hamiltonian. 

The  operator  H*  may  be  thought  of  as  "inducing  transitions  from 
one  allowed  state  to  another"  (84) .  If  the  initial  state  of  the  system 
(prior  to  transition)  is  labeled  ^  and  the  final  state  of  the  system 
(after  the  transition)  is  labeled  i|if,  then  a  measure  of  how  much  the  wave 
function  H'4^  is  like  the  wave  function  ^  is  given  by  the  overlap  inte¬ 
gral 

M  -  iH'i^  *=  J \f*Ef^±dT  (B3) 

where  dr  is  an  element  of  volume  and  the  integral  is  over  all  space.  The 
quantity  M  is  called  the  "transition  matrix  element."  The  matrix  element 
is  an  important  quantity  since  its  square  is  proportional  to  the  prob- 
ability  that  the  perturbation  will  transform  the  state  described  by  ^  to 
the  state  described  by  The  matrix  element  contains  the  wave  functions 
of  the  states  between  which  the  transition  occurs.  It  is  therefore  the 
quantity  that  determines  the  effect  of  the  spins,  parities,  and  detailed 
structures  of  the  two  states  on  the  transition  probability. 

To  calculate  the  probability  for  a  transition,  an  expression  for 
H*  is  needed.  H*  is  the  energy  of  the  interaction  between  the  electric 
and  magnetic  fields  associated  with  the  incoming  photon  and  those  associ¬ 
ated  with  the  target  nucleus.  These  electric  and  magnetic  fields  may  be 
expressed  as  an  expansion  in  terms  of  their  component  electric  and  mag¬ 
netic  multipoles.  This  expansion  may  be  better  illustrated  by  reference 
to  Figure  Bl.  The  electrostatic  potential  V  et  the  point  P  due  to  a  pro¬ 
ton  of  charge  e  at  the  point  (r,4>)  is  (84) 
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Figure  Bl.  Definition  of  quantities  relevant  to 
the  discussion  of  nuclear  electric  multlpoles. 
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V  -  ~~  e/  [d[1  -  (2r/d) cos©-  +  rW]*5]  .  (B4) 

o  1 

This  equation  is  usually  expressed  In  terms  of  the  Legendre  polynomials 
P^CcosQ*)  using  the  Identity 
00 

L  (f)Pn(c<ia^  “  i/U  -  (2r/D)cos-&  +  rZ/D2]^  (B5) 


in  which  case  the  expression  for  V  becomes 


4..r 


[  h  ® 


(cosS) 


This  equation  is  the  electric  multipole  expansion.  To  calculate  the  poten¬ 
tial  at  P  due  to  protons  in  the  nucleus,  it  is  necessary  to  sum  the 
potential  at  P  due  to  each  separately.  After  summing  the  potential  at  P 
due  to  all  Z  protons  in  a  nucleus,  the  first  term  of  Eq.  (B6)  becomes 
Ze/D,  the  electrostatic  potential  that  would  result  by  a  total  charge  Ze 
located  at  the  nuclear  center  of  mass.  The  second  term  of  the  expansion 
becomes  proportional  to 


I 


1 


ercosO  » 


(B7) 


This  quantity  ia  the  electric  dipole  moment.  The  third  term  of  the  expan¬ 
sion  gives  the  electric  quadrupole  moment,  the  fourth  the  electric  octu- 
pole,  and  so  cn.  A  similar  expansion  can  be  made  for  the  magnetic  multi- 
poles. 

For  each  term  of  the  multipole  expansion  of  H',  a  matrix  element 


can  be  calculated  such  that  the  square  of  its  absolute  value  Is  propor¬ 
tional  to  the  probability  that  the  transition  occurred  by  interaction  of 
the  photon  field  with  that  particular  electric  or  magnetic  multipole 
moment  of  the  nuclear  field.  The  total  probability  of  a  transition  per 
unit  time  depends  on  the  squares  of  the  values  of  all  the  electric  and 
magnetic  multipole  transition  matrix  elements  between  the  initial  and 
final  nuclear  states. 

There  are  csrtain  simple  rules  that  determine  which  matrix 
elements  may  contribute  to  the  transition  in  the  particular  case.  The 
electric  L-pole  (EL)  matrix  elements  and  the  magnetic  L-pole  (ML)  matrix 
elements  are  always  zero  unless  the  spins  of  the  Initial  and  final  states 
differ  vectorially  by  L  units  of  angular  momentum.  In  other  words,  elec¬ 
tric  dipole  (El)  or  magnetic  dipole  (Ml)  transitions  may  occur  only 
between  states  of  spin  differing  by  one  unit  of  angular  momentum,  quad- 
rupole  transitions  (E2  and  M2)  may  occur  only  between  states  differing  by 
two  units  of  angular  momentum,  etc.  Thus,  if  and  represents  the 
spins  of  the  initial  and  final  states,  respectively,  and  2**  Indicates  the 
pole  order  (L  -  1,  dipole;  L  ■  2,  quadrupole;  L  ■  3,  octupole;  and  so  on), 
the  first  transition  "selection  rule"  is 

[ Ji  -  Jf]  <  L  <  4  Jf  .  (B8) 

The  second  transition  selection  rule  is  that  of  conservation  of 
parity.  An  electric  multipole  has  even  parity  when  L  is  even  and  odd 
parity  when  L  is  odd.  Magnetic  multipole  radiation  has  odd  parity  when 
L  is  even  and  even  parity  when  L  is  odd.  In  general,  parity  of  electric 
multipoles  is  it  ■  (-1)^,  and  parity  for  magnetic  multipoles  is  v  m  -(-1)^ 
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whore  +1  means  even  parity  and  -1  means  odd  parity  (92).  If  the  parity 
of  the  Initial  state  Is  deaoted  by  and  the  parity  of  the  final  state 
Is  given  by  ir^,  then  If  a  particular  radiation  multipole  has  a  parity 
change  ir^,  conservation  of  parity  would  require  that 

W1  "  1TfTT'-  *  (B9) 

This  Is  the  second  transition  selection  rule.  Transitions  are  forbidden 
between  two  levels  which  both  have  spin  zero.  Table  B1  (92)  summarizes 
the  transition  selection  rules. 

The  transition  rate,  X  between  an  initial  and  a  final  state 
Is  given  by  "Fermi's  Golden  Sule  #2" 


In  this  expression  M  is  the  transition  matrix  element  and  dN/dE  is  the 
number  of  final  states  per  unit  energy.  Blatt  and  Welsskopf  (93)  have 
estimated  electric  multipole  and  magnetic  multipole  transition  rates 
(Xg^  and  X^)  based  on  the  independent  particle  shell  model  of  the  nucleus. 
The  expressions  they  obtained  are  (85,92) 

e2  /R\2L 

xel  - 2™  r s  (I)  (B11) 


where  o  Is  the  radiation  frequency,  e  is  the  electron  charge,  R  Is  a 
linear  dimension  characterizing  nuclear  size,  K  is  the  radiation's  reduced 
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Table  B1 


Transition  Selection  Rules 


Classification 

Symbol 

Change  in  Angular 

Momentum 

Parity  Change 

Electric  dipole 

El 

1 

Yes 

Magnetic  dipole 

Ml 

1 

Mo 

Electric  quadrupole 

E2 

2 

Mo 

Magnetic  quadrupole 

M2 

2 

Yes 

Electric  octupole 

E3 

3 

Yes  * 

Magnetic  octupole  ■■ 

M3 

3 

Mo 

Electric  2^  -  pole 

EL 

L 

Mo  (L  even) 
Yes  (L  odd) 

Magnetic  2L  -  pole 

ML 

L 

Yes  (L  even) 
Mo  (L  odd) 
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wavelength  (<  ■  \/Zv)t  c  Is  the  velocity  of  light,  M  is  the  nucleon  mass 
and  S  is  a  statistical  factor  that  depends  on  L  and  la  shown  in  Table  B2 
(92). 

Inspection  of  Eqa.  (BllfB12)  reveals  the  relative  proportions 

each  multipole  may  contribute  to  a  given  transition.  In  general,  the 

transition  rate  of  the  multipole  of  lowest  order  usually  exceeds  that  of 

*)  A 

all  the  other  multipoles  by  a  factor  of  at  least  l(j  to  l(r  (85) .  Addi¬ 
tionally,  for  a  given  value  of  L,  electric  multipole  radiation  is  more 
probable  than  magnetic  multipole  radiation  by  a  factor  of  approximately 
4.4A2/3  (92). 
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Table  B2 

Value*  of  the  statistical  factor  S 
aa  a  function  of  multlpole  order  L. 


L 


S 


1 

2 


(2.5)  Kf1 
(4.8)  10“3 

(6.3)  10”5 

(5.3)  10“7 
(3.1)  10~9 
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The  linear  accelerator  used  In  this  study  was  a  Saglttalre/Therac 
40  electron  linear  accelerator  manufactured  by  CGR-MeV  In  France  (94). 

This  accelerator  Is  a  disk— loaded,  traveling-wsve  clinical  machine.  The 
description  preaented  In  this  Appendix  pertains  to  traveling-wave  accel¬ 
erators  like  the  Saglttalre.  The  discussion,  however,  has  been  kept  of 
sufficient  generality  to  apply  to  any  traveling-wave  clinical  machina. 

An  electron  linear  accelerator  Is  a  device  that  accelerates  elec¬ 
trons  by  virtue  of  an  axial  electric  field.  When  the  energy  at  Injection 
is  negligible,  the  final  electron  kinetic  energy  obtained  la  given  by 

Te  "  1  (eVdl)  (C1> 

o 

where  Tft  ■  electron  kinetic  energy,  e  -  electron  charge,  E^  -  axial  elec¬ 
tric  field,  dl  •  lucre  meat  In  length,  and  L  ■  accelerating  tube  length. 

Tha  required  axlel  electric  field  le  provided  by  electromagnetic  weves  of 
microweve  frequency,  commonly  3000  MHz. 

Evolution  of  the. modern  electron  linear  accelerator  has  been  made 
possible  by  the  development  of  high-power  red io frequency  generatore  and 
apecially  designed  microwave  cavities  and  waveguides.  Because  of  the 
Important  role  of  microwave  technology  in  llneer  accelerators,  a  descrip¬ 
tion  of  microweve  cavltlee,  weveguidaa,  and  power  eources  will  be  provided 
In  the  succeeding  peregrephs  (95-97).  Also  Included  le  a  discussion  of 
the  prlnclplee  of  electron  acceleration  by  microvavea  with  particular 
emphasis  pieced  on  traveling -wave  high-energy  electron  linear  accelera¬ 
tors.  Components  characteristic  of  medical  accelerstora  alao  will  be 


described. 
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Microwave  cavities  and  guides  support  certain  specific  electro¬ 
magnetic  wave  configurations  and  are  used  in  microwave  power  generation, 
microwave  transport,  and  for  electron  acceleration.  Energetic  electrons 
can  generate  microwave  power  by  exciting  the  cavities  of  microwave  power 
amplifiers  such  as  the  klystron.  On  the  other  hand,  a  linear  array  of 
cavities  can  accelerate  electrons  by  transferring  electromagnetic  energy 
from  the  cavity  fields  to  an  electron  beam.  Microwave  cavity  geometry 
and  dimensions  determine  the  particular  application  of  the  cavity.  Rec¬ 
tangular  waveguides  are  used  primarily  for  high-power  microwave  energy 
transport  while  cylindrical  cavities  are  used  in  microwave  power  amplifi¬ 
cation  and  electron  acceleration. 

The  rectangular  waveguide  is  the  principal  means  by  which  micro- 
wave  energy  is  transported  long  distances.  Rectangular  waveguides  are 
used  to  transmit  microwaves  from  the  power  amplifier  to  the  electron 
acceleration  section  of  a  linear  accelerator.  This  type  of  waveguide  is 
operated  in  what  is  known  as  the  transverse  electric  (TE)  mode.  In  the 
TE  mode,  the  electromagnetic  wave  has  no  electric  field  component  in  the 
direction  of  travel.  Although  microwave  guides  in  general  will  transport 
vavea  of  an  infinite  number  of  resonant  frequencies,  guide  dimensions  are 
designed  for  transport  in  the  most  convenient  frequencies.  In  many  cases 
this  is  the  TCj.0  mode.  The  power  transported  is  given  by  the  vector  cross 
product  of  the  transverse  electric  field  and  the  magnetic  field  Intensity. 

Cylindrical  cavities,  in  electron  linear  accelerators,  are  used 
in  both  microwave  power  amplification  and  in  electron  acceleration.  Simi¬ 
larly  to  the  rectangular  cavity,  the  lowest  frequency,  and  most  dominant, 
modes  often  are  used.  Cylindrical  cavities  are  generally  operated  in  what 


This  mode,  characterized 


is  called  the  transverse  magnetic  mode,  TM^. 
by  an  axial  electric  field  and  an  azimuthal  magnetic  Induction,  is  parti¬ 
cularly  suitable  for  transferring  power  to  electron  beams.  In  addition, 
microwave  power  amplification  is  usually  accomplished  in  high-energy 
accelerators  by  the  klystron,  a  device  utilizing  TM^-raode  cylindrical 
cavities . 

The  klystron,  in  its  simplest  form,  consists  of  two  interconnected 
cylindrical  cavities,  the  first  of  which  is  called  a  buncher  and  the 
aecond  a  catcher.  Microwaves  from  a  radiofrequency  oscillator  are  trans¬ 
ported  to  the  buncher  cavity  of  the  klystron.  Electrons,  thermionically 
emitted  from  a  hot  filament,  enter  the  buncher  cavity  and  are  grouped  or 
"bunched"  by  the  alternating  axial  electric  field.  Bunched  electrons  are 
accelerated  further  by  a  >3C  field  and  travel  through  a  drift  tube  to  the 
catcher  cavity.  The  catcher  cavity,  resonant  at  the  arrival  frequency  of 
the  electron  bunches,  extracts  a  large  portion  of  the  electron  kinetic 
energy  converting  it  to  electromagnetic  energy.  As  a  result,  mlcrowavea 
of  increased  amplitudea  are  produced.  Even  greater  power  amplification 
can  be  accomplished  by  klystrons  of  multiple  cavities. 

The  cylindrical  waveguides  used  for  electron  acceleration  may  be 
thought  of  as  a  linear  array  of  many  cylindrical  cavities.  This  is  accom¬ 
plished  by  Introducing  circular  '‘washer-shaped"  disks  called  irises. 
Iris-loading  of  the  cylindrical  waveguide  is  necessary  to  reduce  the  phase 
velocity  of  the  traveling  wave  in  the  waveguide  to  less  than  the  velocity 
of  light.  This  permits  the  coupling  of  the  traveling-wave  energy  to  the 
electrons  under  acceleration.  The  phase  velocity  of  the  traveling  wave 
varies  linearly  with  the  spacing  between  the  irises,  thus,  close  iris 
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spacing  means  smaller  phase  velocity.  On  the  other  hand,  because  the 
axial  electric  field  In  the  waveguide  accelerates  the  electrons.  Iris 
spacing  must  be  increased  along  the  length  of  the  waveguide  so  that  the 
microwave  phase  velocity  increases  in  synchrony  with  the  electron  velocity 
thus  establishing  phase  stability.  Eventually,  as  the  electron  velocity 
approaches  the  velocity  of  light.  Iris  spacing  becomes  constant.  Conse¬ 
quently,  electrons  injected  into  an  iris-loaded  cylindrical  waveguide  of 
proper  dimensions  ere  accelerated  to  kinetic  energies  which  depend  upon 
the  magnitude  of  the  axial  electric  field  and  the  length  of  the  acceler¬ 
ating  guide  as  shown  In  Eq.  (Cl). 

The  remaining  paragraphs  describe  the  operetlon  and  the  major 
components  of  e  typical  traveling-wave  medical  lineer  accelerator.  This 
description  is  best  provided  by  reference  to  Figure  Cl  which  depicts, 
schematically,  the  basic  components.  The  accelerator  shown  on  Figure  Cl 
and  subsaquently  described  possesses  two  acceleretlon  sections  and  a 
phasa-shifter  component  which  are  not  typical  of  smaller  linear  accelera¬ 
tors.  The  basic  accelerator  may  ba  divided  into  three  major  components: 
the  microwave  power  supply  component,  the  electron  injector  and  accelera¬ 
tor  section,  and  the  gantry  and  treatment  head  components. 

Microwave  pcerar  Is  provided  by  an  oscillator,  wh'ch  supplies  3000 
iGz  (in  S-band  accelerators)  low-power  microwaves,  and  a  microveve  ampli¬ 
fier,  io  tills  case  a  klystron,  ruisas  iron  a  100  Hz  oscllletor,  which 
serves  as  a  "traffic  director"  in  the  accelerator,  control  a  modulator 
which  in  turn  controls  a  hlgh-power  switching  device  celled  e  thyretron. 
The  thyretron  transfers  high-voltege  DC  power  to  the  klystron.  At  the 
same  time  the  3003  MHz  oscillator  delivers  microwaves  to  the  klystron  for 
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Figure  Cl.  Block  dlagraa  of  linear  accelerator 


amplification.  The  amplified  microwaves  from  the  klystron  are  tranaported 
to  the  acceleration  section  by  means  of  rectangular  waveguides  filled  with 
a  dielectric  gas  such  as  SF^. 

The  acceleration  section  consists  of  iris-loaded  cylindrical  wave¬ 
guides.  Electrons  are  injected  into  the  acceleration  guides  by  an  elec¬ 
tron  gun  which  delivers  electron'  pulses  et  repetition  rates  also  controlled 
by  the  100  Hz  oscillator.  The  electrons  are  then  accelerated  by  the 
amplified  3000  MHz  microwaves  traveling  through  the  waveguide.  The 
increase  in  electron  momentum  is  given  by  (97) 

3t  "  cEozXo(ar)co8*  -  (C2) 

where  dp/dr  is  the  time  rete  of  change  of  the  electron  momentum,  e  is  the 
electron  charge,  E  is  the  amplitude  function  for  the  axial  electric 
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field,  Iq (ar)  is  e  modified  Bessel  function  of  the  argument  or,  and  +  is 
t*e  phase  angle.  The  veriable  r  describes  the  redlel  dimension  in  the 
ecce  levs  ting  cylindrical  waveguide.  The  final  electron  kinetic  energy, 
as  can  ba.sean  from  the  above  expression,  is  e  function  of  the  time  the 
electron  experiences  the  Influence  of  the  electric  field  (and  hence  total 
waveguide  length),  the  magnitude  of  the  axial  electric  field,  and  of  the 
value  of  the  eoa+  term  which  has  to  do  with  the  phase  coupling  of  the 
traveling  waves  in  the  guide. 

In  the  accelerator  shown  in  Figure  Cl,  electron  endpoint  energy 
is  controlled  primarily  by  varying  the  relative  phase  between  traveling 
waves  in  the  two  accelerating  sections.  This  is  accomplished  by  changing 
the  phase  velocity  of  the  t reveling  wave  in  thet  part  of  the  rectangular 
waveguide  which  feeds  the  first  acceleration  section.  The  component 
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called  the  phase-shifter  performs  this  task.  The  phase-shifter  consists 
of  a  quartz  rod  and  a  mechanism  capable  of  adjusting  continuously  the 
penetration  distance  of  the  rod  into  the  rectangular  waveguide  feeding 
section  1.  The  degree  of  penetrstion  determines  the  dielectric  constant 
of  that  section  of  the  waveguide  thus  establishing  the  phase  velocity  of 
the  wave  in  that  region  and  shifting  the  relative  phase  of  the  microwaves 
in  the  first  section  with  respect  to  the  second  section.  Adjusting  the 
phase-shifter  changes  the  phase  angle  $  for  the  electrons  in  the  second 
acceleration  section  and  changes  the  electron  momentum  gain.  In  addition t 
magnetic  induction  coils  are  wrapped  around  both  accelerating  sections  in 
order  to  focus  the  electrons  into  an  electron  beam  along  the  centerline 
of  the  sections. 

Accelerated  electrons  enter  the  gantry  and  therapy  head  eection 
of  the  linsc  when  the  beam  undergoes  the  necessery  modifications  required 
to  produce  an  electron  or  bremsstrahlung  beam.  The  gantry  contains  a 
series  of  bending  magnets  and  quadrupoles  thet  transport  and  eteer  the 
electron  beam  and  analyze  the  energy.  The  gantry  also  serves  to  support 
the  therapy  heed  which  contains  components  such  as  the  verlable  collima¬ 
tors,  the  bremsstrahlung  target  and  beam-flattening  filter. 

Upon  entering  the  gantry,  the  electron  beam  energy  Is  analyzed 
by  means  of  a  bending  electromagnet  and  energy  elit  system  (94).  For  e 
given  magnetic  field,  electrons  possessing  cettaiu  dutcreie  energies  are 
deflected  through  some  definite  angle.  Only  electrons  with  energies 
within  a  certein  percent ege  of  a  mean  energy  are  trensmitted  through  the 
energy  elit  opening.  Microveve  power  is  controlled  by  a  eervo  energy 
slit  eystem  such  thet  the  electron  beam  energy  is  always  adjusted  to  the 
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energy  appropriate  for  transmission  through  the  slit.  Electron  endpoint 
energy  selection  determines  the  magnetic  field  Intensity  by  establishing 
the  magnet  deviation  current.  The  resultant  electron  beam  energy,  T,  Is 
given  by  (59) 

<T2  +  ImA)**  -  Ber  (C3) 

Where  T  la  the  electron  kinetic  energy,  m  Is  the  electron  mass,  e  is  the 
electron  charge,  c  la  the  speed  of  light,  B  Is  the  magnetic  Induction, 
and  r  la  the  electron  orbit  radius.  After  energy  analysis,  the  electron 
beam  Is  redirected  by  means  of  additional  bending  magnets  to  a  direction 
perpendicular  to  the  ecceleration  waveguide.  The  beam  then  enters  the 
therapy  heed. 

The  llnac  therapy  head  contains  the  elements  that  perform  the 
final  beam  modifications  prior  to  Its  exit  from  the  accelerator.'  These 
ere  the  fixed  and  variable  beam  collimators  which  restrict  the  beam  size, 
e  removable  target  which  allots  bremsstrahlung  beam  production,  a  beam 
flattening  filter  which  creates  a  uniform  beam  Intensity  In  e  plane  per¬ 
pendicular  to  the  beam  axis,  and  an  Ionisation  chamber  system  whose 
purpose  Is  to  monitor  beam  position  and  Intensity. 

The  Ionization  chamber  system  consists  of  two  pancake-type  Ioni¬ 
zation  chambers  each  divided  diametrically  In  two.  The  chambers  are 
stacked  one  on  top  of  the  other  such  that  the  beam  traversas  both  of 
them.  The  chashers  also  are  rotated  90°  with  respect  to  each  other  such 
that  each  of  the  four  halves  can  be  used  to  monitor  the  Intensity  of  the 
beam  In  one  of  four  perpendicular  directions.  The  signals  from  each  of 
the  four  halves  ere  used  as  Input  to  servo  controls  for  beam-steering 


magnetic  quadrupoles,  located  in  the  gantry.  These  set  to  center  the 
beam.  In  addition  to  other  uses,  the  loni2stlon  chambers'  current  is 
compared  to  a  standard  current  supply  and  the  electron  gun  current  is 
automatically  adjusted  so  that  the  beam  Intensity  remains  constant. 
Finally,  the  loni2stion  chambers  serve  to  monitor  the  total  Integrated 
radiation  dose  delivered. 
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APPENDIX  D 


Semiconductor  Detector  Spectrometer  Syat 
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Gamma-ray  spectrometers  are  radiation  measurement  systems  used 
to  determine  the  energy  spectrum  and  fluence  of  gamma-  or  x-ray-emlttlng 
radioisotopes.  Semiconductor  detector  systems  are  particularly  well- 
suited  for  this  purpose.  They  possess  the  advantageous  characteristics 
of  high  sensitivity  and  superior  energy  resolution  which  allow  accurata 
spectroscopic  measurement. 

In  this  appendix,  the  genaral  characteristics  of  semiconductor 
materials  and  the  application  of  these  materials  to  radiation  detection 
are  reviewed  (63,98-103).  The  components  and  basic  principles  of  pulse- 
height  spectrometry  systems  also  are  described.  Particular  emphasis  is 
placed  on  Intrinsic  germanium  semiconductor  detectors  and  on  multichannel 
pulse-height  analysis  systems. 

In  many  radiation  detection  applications,  such  as  gamma-ray 
spectroscopy,  the  use  of  a  solid  detection  medium  Is  preferred  over  a 
gas-filled  chamber  because  of  the  Increased  sensitivity.  The  use  of 
scintillation  detectors  for  spectroscopy  Is  fairly  common.  Their  energy 
resolution,  however.  Is  rather  limited  primarily  due  to  the  large  contri¬ 
bution  that  statistical  fluctuations  make  to  the  Information-carrying 
pulse.  The  only  way  to  reduce  the  Influence  of  statistical  fluctuations 
on  the  energy  resolution  of  a  detector  la  to  Increase  the  number  of  Infor¬ 
mation  carriers  Involved  In  pulse  production. 

The  use  of  semiconductor  detectors  results  in  a  large  number  of 
Information  carriers  per  detection  event  and  a  very  small  amount  of  sta¬ 
tistical  fluctuation  In  the  procasslng  of  the  charge  pulse.  Thus,  much 
better  energy  resolution  Is  achieved.  For  this  reason,  semiconductor 
detector  systems  are  best  suited  for  such  radiation  measurement  appllca- 
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tlOQB  as  garama-ray  spectroscopy  where  accurate  determination  of  the 
photon  energy  Is  required. 

Electrons  in  an  Individual  atom  can  assume  only  certain  discrete 
energy  levels.  When  atoms  ere  brought  into  close  proximity,  such  as  in 
a  crystal  structure,  these  discrete  energy  levels  tend  to  "smear"  into 
energy  bands.  The  most  energetic  electrons  which  ere  still  bound  to 
specific  lattice  sites  within  the  crystal  occupy  what  is  known  as  the 
valence  band.  Electrons  which  are  free  to  migrate  through  the  crystal, 
and  thus  contribute  to  the  material's  electrical  conductivity,  occupy 
even  higher  energy  states  grouped  into  the  conduction  band.  The  valance 
and  conduction  bands  are  separated  by  e  gap  of  normally  forbidden  energy 
states. 

The  size  of  this  band  gap  determines  whether  a  material  is  an 
insulator  or  a  semiconductor.  In  semiconductor*  the  band  gap  energy  is 
of  the  order  of  several  electron-volts  (eV),  whereas  insulators  have  band 
gaps  greeter  than  5  eV  (63).  In  the  absence  of  thermal  energy,  the 
valence  band  is  completely  filled  and  the  conduction  bend  is  completely 
empty.  Because  of  the  relatively  small  band  gap  energy  of  semiconductors, 
in  the  presence  of  induced  excitations,  some  valence  band  electrons  may 
acquire  sufficient  energy  to  cross  the  band  gap  to  the  conduction  band. 

An  electron,  now  in  the  conduction  band,  and  its  associated  vacancy  in 
the  valence  band  form  what  is  known  as  an  electron-hole  pair.  The  pre¬ 
sence  of  electron-hole  pairs,  produced  by  ionizing  radiation,  changes  the 
electrical  conductivity  of  the  material.  Thus  electron-hole  pairs  con¬ 
stitute  the  beslc  information  carriers  of  semiconductor  detectors. 

The  presence  of  impurities  effects  the  conductivity  of  aemicon- 
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ductors.  In  a  high  purity  or  Intrinsic  semiconductor,  the  number  of 
electrons  In  the  conduction  band  equals  the  number  of  holes  In  the  valence 
band.  In  relatively  pure  semiconductors  such  as  germanium  and  silicon, 
however,  even  minute  quantities  of  Impurities  can  determine  their  elec¬ 
trical  properties.  Impurities  may  be  electron  donors  or  electron  eccep- 
tors  depending  on  their  valence  electron  number.  In  the  case  of  an  elec¬ 
tron  donor  Impurity,  excess  loosely-bcund  electrons  remain  after  the 
impuiity  has  formed  covalent  bonds.  A  semiconductor  containing  this  type 
of  Impurity  is  termed  n-type  and  its  electrical  conductivity  Is  determined 
by  the  flow  of  electrons.  If  the  iivcrity  Is  an  electron  acceptor,  empty 
covalent  bond  sites  result  in  the  formation  of  holes.  In  these  materials, 
the  holes  determine  the  electrical  conductivity  and  the  semiconductor  is 
said  to  be  p-type. 

If  donor  and  acceptor  impurities  are  present  in  a  semiconductor 
in  equal  concent ret ions,  the  material  is  said  to  be  compensated.  Compen¬ 
sated  semiconductor  mate  riels  have  some  of  the  properties  of  Intrinsic 
semiconductors,  i.e.  equal  amounts  of  electrons  and  holes.  Compensation 
may  be  accomplished  by  the  lithium  Ion  drifting  process  in  tdiich  donor 
atoms  (lithium)  are  added  to  basically  p-type  semiconductor  materials 
such  as  gerrnaaim  or  silicon.  The  widely  used  Ge(Ll)  and  S1(L1)  semicon¬ 
ductor  detectors  ere  made  in  this  way. 

Directly-ionizing  radiation  deposits  energy  in  semiconductor 
materials  resulting  in  the  production  of  electron-hole  peirs.  The  number 
of  electron-hole  pairs  or  charge  carriers  produced  Is  determined  by  the 
amount  of  energy  deposited.  The  total  nuafcer  of  charge  carriers  In  the 
semiconductor,  however,  will  not  be  determined  exclusively  by  the  energy 
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deposited  but  will  also  Include  ’’background"  charge  carriers  due  to 
thermal  excitations  and  to  the  presence  of  even  minute  amounts  of  impuri¬ 
ties.  Background  carrier  concentration  can  be  reduced  by  applying  an 
electric  field  to  the  semiconductor  detector  region.  Care  must  be  taken 
when  applying  the  field  in  order  to  be  sure  that  cherge  carriers  removed 
at  one  end  of  the  semiconductor  detector  region  are  not  reintroduced  into 
the  region  from  the  other  end.  For  this  purpose,  speciel  "blocking  con¬ 
tacts”  are  required.  The  most  suiteble  blocking  contact  is  the  junctlt.t 
of  p-type  and  n-type  semiconductor  materials .  This  junction  is  referred 
to  as  the  p-n  junction. 

When  a  p-type  semiconductor  is  joined  to  an  n-type  semiconductor, 

a  concentration  gradient  is  created  in  the  conduction  band  at  the  junction 

due  to  the  increased  density  of  conduction  electrons  in  the  n-type  mater- 
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ial.  The  gradient  Insures  thet  the  electrons  will  diffuse  ecross  the 
junction  to  the  p-type  material.  Bare  the  electrons  will  combine  with 
the  holes.  A  similar  and  symmetrical  argument  can  be  made  for  the  migra¬ 
tion  of  holee  from  the  p-type  to  the  n-type  semiconductor  material.  This 
charge  migretlan  reeulte  in  a  net  negetlve  spece  charge  in  the  p-type 
side  of  the  junction  and  e  net  poeitive  epaee  cherge  in  the  n-type  side. 

The  residual  space  charge  produces  an  electric  field  which  oppoees  further 
diffusion  across  the  junction  thus  establishing  e  steady-*stete  charge 
distribution.  The  region  of  the  cherge  imbalance  is  called  the  depletion 
region  because  the  concentration  of  electrons  and  holee  is  greatly  reduced. 

If  an  electric  potent iel  1s  applied  ecross  a  p-r.  junction  euch 
thet  the  p-type  side  is  made  negative  and  the  n-type  eide  is  made  positive, 
the  junction  is  eeid  to  be  reverse  biased.  Under  this  condition,  the 
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applied  potential  enhances  the  junction's  natural  space  charge  distribu¬ 
tion  and  thus  increases  the  potential  difference  across  the  junction. 

This  ultimately  results  in  a  wider  depletion  region  within  which  charge 
carriere  that  nay  be  produced  through  the  action  of  ionizing  radiation 
can  be  efficiently  collected.-  With  the  potential  difference  the  junction 
also  functions  as  e  blocking  contact  because  the  re  introduction  of  charge 
carriexs  from  outside  the  depletion  region  is  prevented  when  charge  car¬ 
riers  produced  within  the  region  ere  collected.  The  depletion  region  of 
e  reverse  biased  p-n  junction  stay,  therefore,  act  as  the  sensitive  volume 
of  a  semiconductor  radiation  detector. 

The  width  of  the  depletion  region  of  e  semiconductor  detector  is 
proportional  to  the  material's  electrical  resistivity  if  the  material  is 
free  from  impurities.  The  reduced  concentration  of  impurity  in  high- 
purity  germanium  (HPGe)  allows  the  formation  of  relatively  large  depletion 
regions.  These  enlarged  depletion  regions  Increase  the  active  volume  of 
the  semiconductor  detector  thus  making  it  suitable  for  the  detection  of 
penetrating  ionising  radiation  such  as  gamma  rays. 

The  semiconductor  properties  of  HTGe  tend  to  be  p-type  due  to 
the  presence  of  small  amounts  of  residual  acceptor  impurities  or  to  accep¬ 
tor  centers  associated  with  the  Ge  crystal  Itself.  The  basic  conflgure- 

4.  + 

tion  of  an  HPGe  detector  is  sometimes  referred  to  as  n  -  r  -  p  where 
n  ,  p  refers  to  high  concentration  s-  end  p-type  scsi conductor  materials 
respectively  end  p  refers  to  the  bulk  HPGe  materiel  Itself. 

HPGe  detectors  ere  often  fabricated  in  whet  le  known  as  coaxial 
geometry  in  order  to  achieve  the  largest  active  volumes.  The  outer  sur¬ 
face  of  the  crystal  (n+-type)  la  formed  by  lithium  diffusion  onto  the 
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crystal  surface.  The  detector’s  depletion  region  is  formed  by  reverae- 
blaslng  this  n+-p  junction.  The  inner  (p+-type)  contect  is  often  formed 
by  removing  the  inner  core  of  the  Ge  crystel,  etching  the  inner  surfeee 
and  then  evaporating  either  gold  or  aluminum  for  electrical  contect  (63). 
If  the  reverse-biasing  voltage  is  high  enough,  the  entire  p-region  of  the 
detector  becomes  totally  depleted  thus  resulting  in  an  active  volume 
almost  equal  to  the  size  of  the  crystal  itself. 

Germanium  detectors  ere  operated  at  cryogenic  temperatures.  This 
is  necassery  to  minimize  leakage  current.  Current  due  to  thsrmally-gener- 
ated  charge  carriers  may  constitute  the  largest  source  of  leakage  current. 
Germanium  detectors  are  specially  susceptible  to  leakage  current  because 
of  the  lower  band  gap  energy,  0.67  eV.  The  current  observed  at  the  edges 
of  a  reverse-biased  junction  when  a  high  voltage  is  applied  constitutes 
another  source  of  leakage.  This  current,  called  surface  leakage',  ie 
negligibly  small  provided  the  detector  has  good  encapsulation. 

The  charge  carriers  produced  in  an  ionization  or  excitation  event 
in  a  semiconductor  detector  produce  a  voltage  pulse  whoee  amplitude  is 
directly  proportional  to  the  amount  of  energy  deposited  in  the  detector. 

An  electronic  eystem  that  ia  capable  of  discriminating  and  recording  the 
distribution  of  aaplltudee  of  a  eeriee  of  pulses  produced  by  a  detector 
is  called  a  pulse-height  spectrometer.  The  pulse-height  epectrometer 
system  generally  consists  of  a  detector,  preamplifier,  an  amplifier,  and 
a  pulse-height  analyzer  frequently  a  multichannel  analyzer.  Theae  compo¬ 
nents  amplify,  shape,  analyze  and  record  pulses  according  to  their  ampli¬ 
tude. 

The  cherge  involved  in  an  individual  detector  output  pulse  is 
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so  small  that  some  Initial  amplification  Is  required  prior  to  processing 
and  analysis.  This  function  Is  accomplished  by  the  preamplifier  which  Is 
usually  located  In  close  proximity  to  the  detector  Itself.  The  preampli¬ 
fier  also  serves  as  an  Impedance  matcher  between  the  detector  and  tha 
other  components.  It  presents  a  high  impedance  to  the  detector  to  mini- 
miss  loedlng  whila  providing  low  impedance  output  to  the  signal  cable  and 
components  which  follow  (101) . 

The  preamplifier  output  pulsa  generally  has  small  amplitude  and 
rapid  rlss  and  long  fall  times.  Linear  amplifiers  are  used  to  convert 
such  a  low  amplitude  pulse  to  one  of  sufficient  amplitude  and  proper  shape 
to  drive  the  pulss  analysis  components  electronics.  Semiconductor  detac- 
tor  pulses  are  generally  shaped,  using  RC  circuit  elements.  Into  Gaussian- 
lilce  pulses  thet  have  rapid  rise  and  rapid  fell  times  (100) . 

The  multichannel  analyzer  (MCA)  Is  he  component  of  the  pulse- 
height  spectrometer  that  discriminates  and  sorts  the  pulses  eccording  to 
their  amplitude.  The  operetion  of  the  MCA  Involves  conversion  of  pulse 
amplitude,  an  analog  signal,  to  an  equivalent  digital  number.  The  conver¬ 
sion  is  accomplished  by  an  analog-to-dlgital  converter  (ADC).  Most  ADC's 
use  the  h’ilktnsoo-type  llneer  r*mp  converter  to  digitize  pulss  height 
although  other  annlog-to-digital  conversion  methods  exist. 

In  the  Wilkins on- type  converter,  e  cepecitor  charges  llneerly  In 
time  to  a  level  which  depends  on  the  amplitude  uf  the  input  pulse.  A 
compere  tor  circuit  comperes  the  magnitude  of  the  eccumulated  cepadtor 
charga  with  the  amplitude  of  the  Input  pulse.  During  the  time  the  caped- 
tor  is  cherging,  clock  pulses  from  e  high-frequency  oscillator  (about  50 
to  100  MHz)  are  counted.  The  number  of  pulses  accumulated  is  proportional 
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to  the  charging  time  of  the  capacitor  which  in  turn  is  proportional  to 
the  input  pulse's  amplitude. 

The  pulse  amplitude  values  from  the  ADC  fell  into  intervals  of 
pulse  height,  or  channels.  The  pulse-height  dlsir;.lnition  is  stored  in  e 
computer-type  memory  which  has  as  many  addressable  locations  as  the  maxi¬ 
mum  number  of  channels.  The  number  of  memory  locations  is  usually  made 
a  power  of  two  with  memories  of  256  to  4096  channels  being  most  common. 
The  distribution  in  frequency  of  pulses  observed  ss  e  function  of  pulse 
height  is  called  the  pulse-height  spectrum.  The  pulse-height  spectnm 
can  be  displayed  or  recorded  on  peripheral  output  devices  such  as  a 
cathode  ray  tube  or  teletypewriter,  can  be  transferred  to  soma  computer- 
accessible  media  such  as  paper  type  or  magnetic  tape,  or  can  be  processed 
by  microprocessor-type  programs  for  spectrum  analysis. 

Intrinsic  germanium  and  multichannel  analyzer  pulse-height  spec¬ 
trometer  systems  are  particularly  well-suited  for  gamma-ray  detection  and 
measurement  applications  where  precise  gamma-ray  energy  determination  is 
necessary.  Intrinsic  germanium  possesses  the  sensitivity  and  energy 
resolution  required  for  such  detailed  analyses.  Modern  multichannel 
analyzer  systems  permit  efficient  data  collection  and  manipulation  so 
that  pulse-height  spectrum  examination  is  facilitated.  Such  systems 
enable  precise  definition  of  the  nature  and  amount  of  radioactive  isotope 
in  a  sample. 
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Thick—' Target  Breuss trahlung  Computer  Program 
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This  appendix  presents  a  listing  of  the  bremsstrahlung  computer 
program  described  In  the  Bremsstrahlung  Spectrum  section  of  Chapter  III. 
The  program  Is  written  In  Data  General  FORTRAN  5.  It  consists  of  a  main 
program  (BREM.FR)  and  seven  subroutines  (INT.FR,  DFS.FR,  AT.FR,  ANG.FR, 
STP.FR,  TRN.FR,  and  PBAT.FR).  Each  is  briefly  documented  by  comment 
statements  appearing  at  the  beginning  of  the  computer  code.  Bremsstrah¬ 
lung  spectra,  contained  In  the  two-dimensional  array  NIJ(I ,J) ,  are  written 
to  the  test  file,  BRMTST.OUT,  to  allow  Inspection  of  the  spectra,  and  to 
the  file  BRM.OUT  which  Is  later  used  as  Input  for  the  Least  Structure 
cross  section  computer  program. 
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c  *************************************************** 

C  ***  PROGRAM  8RCM.FR  *** 

C 

C  PROGRAM  TO  GENERATE  BREM  SPECTRA  FROM  THICK  TARGETS 

C  FOR  PKOTONUCLEAR  CROSS-SECTION  DETERMINATION. 

C  DATA  GENERAL  FORTRAN  5  CODE  IS  USED. 

C  ' 

C  KARL  PRADO,  MARCH  19B3 

c  ******************************************,>******** 

c 

REAL  NIJ(35,35) ,  B (50) ,  K (SO) ,  NK(SO) , 

.  NPCTIJ (35,35) ,  NPCTK(SO) 

DIMENSION  TITLE (5) 

DATA  E/50*0.0/,  K/50*0.0/,  NK/50*0.0/, 

.  NPCTK/50*0 .0/ 

C 

DO  11  1-1,35 
DO  11  J-1,35 
NIJ(£,J)-0.0 
NPCTIJ (I,J)-0.0 
11  CONTINUE 
C 

C  OPEN  INPUT  AND  OUTPUT  FILES 

C 

OPEN  5,  'BRM.IN' 

OPEN  6,  * BRM.OUT' 

C 

C  TEST  FILE  BRMTST.OUT 

C 

OPEN  7,  'BRMTST.OUT' 

C 

C  READ  THE  INPUT  FILE 

C  1.  TITLE 

C  2,  EMIN/DE 

C  3.  NPTS 

C 

READ  (5,190)  TITLE 
190  FORMAT  (5A4) 

READ  (5,290)  EMIN,DE 
290  FORMAT  (F7.4,1X,F€.4) 

READ  (5,390)  NPTS 
390  FORMAT  (13) 

C 

C  DEFFINE  TARGET  PARAMETERS 

C 

TTHICK-7 .72 
DELTAT-0 .0366 
NS LABS- 200 
XFLAT-3.0S5 
C 

C  ENDPOINT  AND  PHOTON  ENERGY  BINS  SET-UP 
C 

XI1-0.0 

DO  10  I-l.NPTS 
E  (I ) -EMIN+XI 1*DE 
K  ( I) ■£ ( I ) -0 . 5*DE 


* 


non  non  non  non  non  non  non  non 
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xn-xmi.o 

10  CONTINUE 

•*>**  START  ENDPOINT  ENERGY  DO  LOOP  **’ 

DO  1  1-1 ,NPTS 
ENDPT-E(I) 

DO  222  IJK-1,50 
222  NPCTK ( I JK) "0 . 0 

**•*  START  PHOTON  ENERGY  DO  LOOP  **** 

DO  2  J-1,I 
CAY-K(J) 

***  21  ENDPOINT  SPREAD  LOOP  *** 

DO  111  I PC T* 1,3 

IF( XPCT. EQ.l) ENDPT- (E (I) 02*E (I) ) 
IF(IPCT.EQ.2)ENDPT-E(I) 

IF (IPCT.EQ. 3) ENDPT-(E (I) +.02*E (I) ) 
lF(IPCT.EQ.l)XPCT-.25 
IF (IPCT.EQ. 2)XPCT«. 50 
IF(IPCT.EQ.3)XPCT».25 

***  START  SLAB  SUMMATION  DO  LOOP  *** 

TEAVE-ENDPT 

XNSUM-0.0 

SWFRST-1.0 

PTHETA-0.0 

XI2-1.0 

DO  3  ISLA8-1 ,N5LA8S  * 

SET  UP  SLAB  DISTANCES 

DISTI- (XI2*OELTAT) - (0 . 5*DELTAT) 

DISTX-TTHICK-DISTI 

XI2-XI2+1.0 

TEAVE  FOR  THE  SLAB 

CALL  STP  (TEAVE* DTDS) 

IF  (SWFRST.EQ.1.0)  GD  TD  71 

TEAVE  FOR  ISLAB  >  1 

DEETEE- (DELTAT*DTDS) /COS (THETA) 

IF  (DEETEE.GT. TEAVE)  DEETEE-TEAVE 

TEAVE-TEAVE-DEETEE 

GO  TD  77 

TEAVE  FOR  ISLAB  *  1 

71  DEETEE- (0 . 5*DELTAT*DTDS) 

IF  (DEETEE.GT. TEAVE)  DEETEE-TEAVE 


on  oonnn  -non  non  non  nno  ooo  non  non  o  non 
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TEAVE“TEAVE -DEBTEE 
SWFRST-0.0 
77  CONTINUE 

CALCULATE  XNT  ft  RMS  ANGLE 

CALL  ANG  (TEA VE , DELTAT , PTHETA, THETA) 

PTHETA-THETA 

AVNO-6.02217E23 
?T0P-AVN0*DELTAT 
BBOT»183.B5*COS (THETA) 

XNT-TTOP/BBOT 

INTRINSIC  SPECTRUM 

CALL  INT  (TEAVE, CAY, DE, THETA, XNSK) 

XNTK-XNSK*XNT 

PHOTON  ATTENUATION  BY  REMAINING  SLABS 

CALL  AT  (CAY ( DISTX , ATPACT ) 

XNTK-XNTK*ATPACT  , 

ELECTRON  TRANSMISSION  FACTOR 

call  trn  (endpt,disti, trans) 

XNTK-XNTK* TRANS 

SUM  SLAB  CONTRIBUTIONS 

XNSUM-XNSUM+XNTK 

***  END  SLAB  summation  DO  LOOP  *** 

3  CONTINUE 

MULTIPLY  BY  SAMPLE  SOLID  ANGLE 
XNPHOT-XNSUM*0. 3288365 

PHOTON  ATTENUATION  BY  FLATTENING  FILTER 

CALL  PBAT  (CAY , XPLAT , PBATEN ) 

XNPHOT-XNPHOT*PBATEN 

WRITE  PHOTON  NUMBER  (XNPHOT)  TO  THE  ELEMENT  (I,J) 

of  The  array  nij,  and  to  the  element  (j)  of  the 

ARRAY  NK. 

NIJ (I  ,J) -XNPHOT*XPCT 

NK(J)-XNPHOT*XPCT 

NPCTIJ ( I , J)"NPCTIJ (I,J)+NIJ(I,J) 

NPCTK(J)-NPCTK(J)+NK(J) 

***  END  2%  LOOP  *** 


non  nnnnn  nnn  n n n  nnn 


111  CONTINUE 


***  END  PHOTON  ENERGY  DO  LOOP  *** 

2  CONTINUE 

WRITE  PHOTON  SPECTRA  TO  TEST  PILE 
WRITE  (7,195)  ENDPT 

195  FORMAT  (5 (/), IX, ' ENDPOINT  ENERGY  -  *,F7.4,/) 
WRITE  (7,295) 

295  FORMAT  (/, IX, ' PHOTON' ,10X, * PHOTONf ,/, IX, 
'ENERGY' ,10X,'NUMBER' ,//) 

DO  44  L»1,I 

44  WRITE  (7,395)  K(L),  NPCTK(L) 

395  FORMAT  (1X,F7.4,4X,E14.6) 

****  ehD  ENDPOINT  ENERGY  DO  LOOP  *** 

1  CONTINUE 

WRITE  BREM  SPECTRA  TO  FILE  BRM.OUT.  FILE 
3RM.0UT  CONTAINS  THE  BREM  SPECTRA  TO  BE 
USED  TO  OBTAIN  THE  PHOTONUCLEAR  CROSS-SECTION 

DO  55  I»1,NPTS 

55  WRITE  (6,495)  (NPCTIJ (I ,J) ,  J«1,I) 

495  FORMAT  (E14.6) 

CLOSE  FILES  AND  END 

CLOSE  5 
CLOSE  6 
CLOSE  7 
CALL  EXIT 
END 


non 
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C  ***  SUBROUTINE  INT.FR  *** 

C 

C  INTEGRATES  THE  DIFFERENTIAL-IN-ANGLE  SCHIFF 

C  SPECTRUM  (FROM  SUB.  DFS.FR)  OVER  ANGLES 
C  UP  TO  THAT  DETERMINED  BY  MULTIPLE  SCATTERING 
C  IN  THE  SLAB  (FROM  ANG.FR) . 

C 

SUBROUTINE  INT  (T,CAY ,DE , RADIAN, SPCTJ 
C 

SPCT-0.0 
DRAD-0. 0174532 
XN-0.0 

SUM-OVER-ANGLES  .DO  LOOP 

DO  2  1-1,44 
XRAD-XN  *DRAD 

IP (XRAD.GT. RADIAN)  GO  TO  300 
CALL  DPS  (T,CAY,DE ,XRADfTMP) 

SPCT-SPCT+TMP 
XN-XN+1 . 

‘  2  CONTINUE 

C 

300  RETURN 
END 


non  o  o  o  nnooo  n 
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**  SUBROUTINE  DFS.PR  ** 

SUBROUTINE  FOR  CALCULATING  SREM.  SPECTRA  BASED 
ON  SCHIPF’S  DIFFERENTIAL-IN-ANGLE  CROSS  SECTION. 
FROMt  LENT  AND  DICKINSON ,  UCRL-50442. 

SUBROUTINE  DFS  (T,CAY,DK, ANGL,SPECTR) 

U-0.511 

2-74.0 

PI-3.14159 

E-T+U 

ROSQ-5.794E-28 
SIGMAO-ROSQ*Z* ( Z+l . ) 

PHI- (E*ANGL) /U 

IF  (CAY.GT.T)  GO  TO  9 
FACTOR-1. +<PHI**2.) 

DIV1-111 .* FACTOR 
DIV2-Z**. 33333 
SUM1- (DIV2/DIV1 ) **2 . 

PROD1-CAY/ (E-CAY)  v 

PROD2-U/ (2.*E) 

SUM2-(PROD2*PRODl)**2.  . 

ONEOVR-SUM1+SUH2 
EMPHI-1 ./ONEOVR 
ENLN-ALOG(ENPHI) 

XKE-CAY/E 

51- l.-XKE 

52-  <  4 . * ( PHI **2 . ) ) /FACTOR* *2 . 

83-l.+(l.-XKE)**2. 

SQRBRK-S3-(S2*S1) 

BRC1-SQRBRK*EMLN 

BRC2- (2 .-XKE) **2. 

Bl- (16.*PHI**2.) /FACTOR**2 . 

BRC3«B1*S1 

BRACE-BRC3-BRC2+BRC 1 
FACT1-1 ./( FACTOR* *2 . ) 

FI- (E/U) **2. 

FACT2- (F1*DK) /CAY 
FACT3- (2.*SIGMAO) /PI 
C0NST-FACT3*FACT2*FACT1 
SPECTR-CON ST* BRACE 
GO  TO  1 
9  SPECTR-0.0 
1  RETURN 

END  OF  SUBROUTINE 
END 


nnnnnnnnonnr.nnnn 
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**  SUBROUTINE  ANG . FR  ** 


SUBROUTINE  ANG  CALCULATES  RMS  ANGLE  (THETA)  FOR 
ELECTRONS  IN  TUNGSTEN.  BICHSEL  IN  ATTIX  RAD.  DOSIMETRY  I 
Z-74,  A-183.85,  2MBC2-1.022,  B-7.35  (INTERPOLATION  TABLE 
XIII,  P.221,  ATTIX  VOL.  I,  FOR  DELTAT  0.0386  GM/CM2) 

MODIFIED  ANGLE  (MODTHETA)  t 

1.  B-B-1.2  (HANSOM  A.O.,  BT  AL,  PHYS  REV  84,  634  (1951)1 
THEREFORE  B-6.15 

2.  MODTHETA* (THETA**2.+PTHETA**2. ) ** .5 
PTHETA*PREVIOUS  SLAB'S  ANGLE 

3.  THETA  LIMITED  TO  0.734  RAD  (42  DEG.) 

(FERDINANDS  B.,  ET  AL,  NUC  INST  METH  91,  135  (1971)) 

4.  WILL  USE  B-B-5.01-2.34  BASED  ON  FWHM  DATA 

SUBROUTINE  ANG  (T ,D,PPHI,PHI) 

XNUMER-0 .157*74 . *75 .*D*2 . 34 
XDENOH*183.85*T* (T+l .022) 

PHISQR-XNUMER/XDENOM 

PBIl-SQRT(PHISQR) 

PBI*SQRT(PHI1**2.'»,PPHI**2.) 

IF (PHI. GT. 0.7340)  PHI-0.7340 

RETURN 

END 
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**  SUBROUTINE  AT.FR  ** 

SUBROUTINE  AT  CALCULATES  THE  PHOTON  ATTENUATION  FACTOR 
EXP- (MU*X)  MU  TUNGSTEN  (STORM  &  ISRAEL,  NUC  DAT  TABLES 
A7 ,  565,  (1970)).  X  (DISTX  IN  MAIN)  IS  IN  GM/CM2  AND  IS 
EQUAL  TO  TTHICK  -  DISTI. 

******  dote  DISCONTINUITY  AT  ABOUT  KK-1.0.  ****** 

AT  KK-1  MU  IS  FOR  1,  AT  KK<1  MU  IS  FOR  0.5,  AT  KK>1  MU 
IS  FROM  CURVE. 

SUBROUTINE  AT  (KK,X, ATNATN) 

REAL  KK,  MU 

FIND  MU (KK) 

IF ( KR.LT. 1.0)  MU-0.1363 
IF(KK.EQ.l.O)  MU-0.0655 

IF (KR.GT. 1.0. AND.KK.LT .4.0)  MU-0 .0493* (ALOGIO (KK) )**2.0 
•  -0.0591* (ALOGIO (KK) ) +0.0577 
IF(KK.GE. 4.0. AND.KK.LT. 20.0)  MU-0.0359* (ALOGIO (KK) ) **2. 
.  -0. 0431* (ALOGIO (KK)) +0.0529 
IF(KK.GE.20.0)  MU-O. 047* (ALOGIO (KK) ) -0.0037 

CALCULATE  ATTENUATION  FACTOR 

EXPORT— 1.0*MU*X 
ATNATN-EXP (EXPORT) 

RETURN 

END 
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**  SUBROUTINE  STP.PR  *• 

SUB.  STP  CALCULATES  STOPPING  POWER  OP  ELECTRONS  IN 
TUNGSTEN.  BERGER  &  SELTZER  NASA  SP-3012 

SUBROUTINE  STP  (EE,DTOS) 

XF(EE.LT.l.O)  DTDS-1.175 

IF (EE. GE. 1.0. AND. EE. LE. 15.0)  DTDS-1.25* (ALOGIO (EE) ) **2. 
.  +1.175 

IF(IE.GT.15.0.AND.EE.LE.30 .0)  DTDS-7. 0433* (ALOGIO (BE) ) 

.  -5'.  391 

IF(BB.GT.3C.O)  DTDS-11 .91* ( ALOG1U (EE) ) -12. 59 

RETURN 

END 
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C 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


**  SUBROUTINE  TRN.FR  ** 

SUB.  TRN  CALCULATES  A  TRANSMISSION  FACTOR  FOR  AN 
ELECTRON  BEAM  OF  ENERGY  ENDPT  THROUGH  A  DISTANCE 
EQUAL  TO  DXSTI.  FROM:  EBERT,  PHYS  REV  183,  422, 
(1969) . 

SUBROUTINE  TRN  (ENDPT, DISTI , TRANS) 
BETA1«387.0*ENDPT 

BETA2-74.0* (1. +0. 000075*74. 0* (ENDPT**2.0) ) 

BETA3-BETA1/BETA2 

BETA*  (BETA3)  **0  .‘i5 

ALPHA* (1 .0- (1 . u/BETA) ) ** ( 1.0-BETA) 

FAC1-(125. 0/(74. 0+112.0) ) *ENDPT*0. 565 
FAC2- (175. 0/(74. 0+162.0)) *0.423 
REX-FAC1-FAC2 

EXPNT* (-1 .0*ALPHA) * ( (DISTI/REX) **BETA) 

TRANS -EXP(EXPNT) 

RETURN 

END 
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***  SUBROUTINE  PBAT.FR  *** 

SUBROUTINE  PBAT  CALCULATES  MU  (CM2/G)  FOR  LEAD 
FOR  THE  PURPOSE  OF  FLATTENING  FILTER  ATTENUATION 
(STORM  AND  ISRAEL,  NUC  DAT  TABLES  A7,  565  (1970)) 

'  SUBROUTINE  PBAT  (XE,XD,XF) 

IF(XE.LT.l.O)  PBMU-.1590 
IP(XE.EQ.l.O)  PBMU-.0703 

IF(XE.GT. 1.0. AND.XE.LT. 3.0)  PBMU-.0746* (ALOGlO (XE) ) **2. 

.  -.0812* (ALOGlO (XE)) +.0637 
IP(XE .GE.3.0.AND.XE.LE.4.0)  PBMU-.0416 

IF(XE.GT.4.0.AND.XE.LE.15.0)  PBMU*. 0328* (ALOGlO (XE) ) **2. 
.  -.0357* (ALOGlO (XE)) +.0513 
IF(XE.GT.IS.O)  PBMU*.04977* (ALOGlO (XE))-. 0038 

PI— 1.0*XD*PBMU 
XF-EXP(Pl) 

RETURN 

END 


210 


APPENDIX  F 

Bremsatrahlung  Spectra  at  25,  24,  22,  20,  18, 
and  15  MeV  Endpoint  Energies 


«  A'A"  a'*  .^V'  -'w  a  *  ■  '  •'  Aj  1  A  A  A  A  A* .  A  A*  . 
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This  appendix  presents  thick-target  bremsstrahlung  spectra,  with 

and  without  the  ±22  endpoint  energy  spread,  at  selected  endpoint  energies, 

obtained  from  the  thick-target  bremsstrahlung  computer  program.  The 

2 

spectra  were  produced  assuming  a  7.72  g/cm  tungsten  target,  a  radiation 

2 

cone  of  0.329  steradians,  and  a  beam  filtration  of  31.78  g/cm  (2.8  cm) 
of  lead.  Although  these  conditions  apply  strictly  to  bremsstrahlung  pro¬ 
duced  by  the  Sagittaire  accelerator,  spectral  shapes  are  considered 
fairly  representative  of  flattened  clinical  beams  generated  at  the  above 
potentials. 
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25.0  MEV  SPECTRUM 


ENERGY 

BIN 

PHOTON 

ENERGY 

PHOTON 

PLUENCE 

ENERGY 

FLUENCE 

1.00 

.50 

1.198E  00 

5.988E-01 

2.00 

1.50 

1.780E  01 

2.670E  01 

3.00 

2.50 

1.355E  01 

3.387E  01 

4.00 

3.50 

9.606E  00 

3.362E  01 

5.00 

4.50 

6.907E  00 

3.108E  01 

6.00 

5.50 

5.127E  00 

2.820E  01 

7.03 

6.50 

3.898E  00 

2.534E  01 

8.00 

'  7.50 

3.019E  00 

2.264E  01 

9.00 

8.50 

2.370E  00 

2.015E  01 

10.00 

9.50 

1.879E  00 

1.785E  01 

11.00 

10.50 

1.500E  00 

1.575E  01 

12.00 

11.50 

1.201E  00 

1.381E  01 

13.00 

12.50 

9.629E-01 

1.204E  01 

14.00 

13.50 

7.705E-01 

1.040E  01 

15.00 

14.50 

6. 135E-01 

8f 896E  00 

16.00 

15.50 

4.829E-01 

7.485E  00 

17.00 

16.50 

3.760B-01 

6.204E  00 

18.00 

17.50 

2.888E-01 

5.055E  00 

19.00 

18.50 

2.176E-01 

4.026E  00 

20.00 

19.50 

1.590E-01 

3.101E  00 

21.00 

20.50 

1.111E-01 

2.277E  00 

22.00 

21.50 

7.179E-02 

1.544E  00 

23.00 

22.50 

4.098E-02 

9.221E-01 

24.00 

23.50 

1.780E-02 

4 . 183E-01 

25.00 

24.50 

3.081E-J3 

7.549E-02 

25.0  MEV  SPECTRUM  (+/“  2%) 


ENERGY 

BIN 

PHOTON 

ENERGY 

PHOTON 

FLUENCE 

ENERGY 

FLUENCE 

1.00 

.50 

1.198E  00 

5.990E- 

-01 

2.00 

1.50 

1.781E  01 

2.671E 

01 

3.00 

2.50 

1.355E  01 

3.388E 

01 

4.00 

3.50 

9.609E  00 

3.363E 

01 

5.00 

4.50 

6.909E  00 

3.109E 

01 

6.00 

5.50 

5.128E  00 

2.821E 

01 

7.00 

6.50 

3.900E  00 

2.535E 

01 

8.00 

7.50 

3.020E  00 

2.265E 

01 

9.00 

8.50 

2.371E  00 

2.015E 

01 

10.00 

9.50 

1.880E  00 

1.786E 

01 

11.00 

10.50 

1.500E  00 

1.575E 

01 

12.00 

11.50 

1.202E  00 

1.382E 

01 

13.00 

12.50 

9.633E-01 

1.204E 

01 

14.00 

13.50 

7.708E-01 

1.041E 

01 

15.00 

14.50 

6.138E-01 

8.900E 

00 

16.00 

15.50 

4.832E-01 

7.489B 

00 

17.00 

16.50 

3.761E-01 

6.206E 

00 

18.00 

17.50 

2.891E-01 

5.059E 

00 

19.00 

18.50 

2. 178E-01 

4.028E 

00 

20.00 

19.50 

1 . 592E-01 

3.104E 

00 

21.00 

20.50 

1 . 113E-01 

2.281E 

00 

22.00 

21.50 

7 . 2CSE-02 

1.549E 

00 

23.00 

22.50 

4 . 129E-02 

9.291E- 

■01 

24.00 

23.50 

1.821E-02 

4.278E- 

•01 

25.00 

24.50 

3.806E-03 

9.326E-02 

24.0  MEV  SPECTRUM 


ENERGY 

BIN 

PHOTON 

ENERGY 

PHOTON 

FLUENCE 

ENERGY 

FLUENCE 

1.00 

.50 

1.10SE  00 

5.528E-01 

2.00 

1.50 

1.647E  01 

2.470E  01 

3.00 

2.50 

1.250E  01 

3.125E  01 

4.00 

3.50 

8.845E  00 

3.096E  01 

..  5.00 

4.50 

6.345E  QO 

2.855E  01 

6.00 

5.50 

4.698E  00 

2.584E  01 

7.00 

6.50 

3.563E  00 

2. il6E  01 

8.00 

7.50 

2.750E  00 

2.063E  01 

9.00 

8.50 

2.151E  00 

1.828E  01 

10.00 

9.50 

1.698E  00 

1.613E  01 

11.00 

10.50 

- 

1.348E  00 

1.416E  01 

12.00 

11.50 

1.073E  00 

1.233E  01 

13.00 

12.50 

8.531E-01 

1.066E  01 

14.00 

13.50 

• 

6.759E-01 

9.124E  00 

15.00 

14.50 

5.320E-01 

7.714E  00 

16.00 

15.50 

4.118E-01 

6.384E  00 

17.00 

16.50 

3.143E-01 

5.186E  00 

18.00 

17.50 

2.35IE-01 

4.114E  00 

19.00 

18.50 

1.7C9E-01 

3.162E  00 

20.00 

19.50 

1.187E-01 

2.314E  00 

21.00 

20.50 

7.658E-02 

1.570E  00 

22.00 

21.50 

4.317E-02 

9.282E-01 

23.00 

22.50 

1.868E-02 

4 . 202E-01 

24.00 

23.50 

3.196E-03 

7.511E-02 

24.0  MEV  SPECTRUM  (+/-  2%) 


ENERGY 

PHOTON 

PHOTON 

ENERGY 

BIN 

ENERGY 

FLUENCE 

FLUENCE 

1.00 

.50 

1.106E  00 

5. 529E-01 

2.00 

1.50 

1.647E  01 

2.470E  01 

3.00 

2.50 

1.250E  01 

3.126E  01 

4.00 

3.50 

8.847E  00 

3.09SE  01 

5.00 

4.50 

6.346E  00 

2.856E  01 

6.00 

5.50 

4.699E  00 

2.584E  01 

7.00 

6.50 

3.564E  00 

2.316E  01 

8.00 

7.50 

2.751E  00 

2.063E  01 

9.00 

8.50 

2.152E  00 

1.829E  01 

10.00 

9.50 

1.699E  00 

1.614E  01 

11.00 

10.50 

1.349E  00 

1.41SE  01 

12.00 

11.50 

1.073E  00 

1.234E  01 

13.00 

12.50 

8 . 534E-01 

1.067E  01 

14.00 

13.50 

6.761E-01 

9.127E  00 

15.00' 

14.50 

5 . 321E-01 

7.716E  00 

16.00 

15.50 

4 . 120E-01 

6.387E  00 

17.00 

16.50 

3. 145E-01 

5.189E  00 

18.00 

17.50 

2.353E-01 

4.117E  00 

19.00 

18.50 

1.711E-01 

3.166E  00 

20.00 

19.50 

1 . 189E-01 

2.318E  00 

21.00 

20.50 

7.681E-02 

1.575E  00 

22.00 

21.50 

4.358E702 

9.369E-01 

23.00 

22.50 

1.907E-02 

4 . 290E-O1 

24.00 

23.50 

3.902E-03 

9 . 169E-02 
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22.0  MEV  SPECTRUM 


ENERGY 

BIN 

PHOTON 

ENERGY 

PHOTON 

FLUENCE 

ENERGY 

FLUENCE 

1.00 

.50 

9 . 300E-01 

4.650E-01 

2.00 

1.50 

1.392E  01 

2.088E 

01 

3.00 

2.50 

1.052E  01 

2.630E 

01 

4.00 

3.50 

7.408E  00 

2.593E 

01 

5.00 

4.50 

5.286E  00 

2.379E 

01 

6.00 

5.50 

3.891E  00 

2.140E 

01 

7.00 

6.50 

2.932E  00 

1.906E 

01 

8.00 

7.50 

2.246E  00 

1.685E 

01 

9. do 

8.50 

1.741E  00 

1.480E 

01 

10.00 

9.50 

1.360E  00 

1.292E 

01 

11.00 

10.50 

1.065E  00 

1.118E 

01 

12.00 

11.50 

8.332E-01 

9.582E 

00 

13.00 

12.50 

6.483E-01 

8.103E 

00 

14.00 

13.50 

4.997E-01 

6.746E 

00 

15.00 

14.50 

3.797E-01 

5.506E 

00 

16.00 

15.50 

2 .805E-01 

4.348E 

00 

17.00 

16.50 

2.011E-01 

3.319E 

00 

18.00 

17.50 

1.378E-01 

2.412E 

00 

19.00 

18.50 

8.805E-02 

1.629E 

00 

20.00 

19.50 

4.944E-02 

9.642E-01 

21.00 

20.50 

2 . 115E-02 

4.336E-01 

22.00 

21.50 

3.453E-03 

7.424E- 

•02 

22.0  MEV  SPECTRUM  (+/-  2%) 


ENERGY 

DIN 

PHOTON 

ENERGY 

PHOTON 

FLUENCE 

ENERGY 

FLUENCE 

1.00 

.50  • 

9.301E-01 

4.650E-01 

2.00 

1.50 

1.392E  01 

2.088E 

01 

3.00 

2.50 

1.052E  01 

2.631E 

01 

4.00 

3.50 

7.409E  00 

2.593E 

01 

5.00 

4.50 

5.287E  00 

2.379E 

01 

6.00 

5.50 

3.892E  00 

2.140E 

01 

7.00 

6.50 

2.932E  00 

1.906E 

01 

8.00 

7.50  ' 

2.247E  00 

1.685E 

01 

9.00 

8.50 

1.742E  00 

1.480E 

01 

10.00 

9.50 

1.360E  00 

1.292E 

01 

11.00 

10.50 

1.065E  00 

1.118E 

01 

12.00 

11.50 

8. 334E-01 

9.584E 

00 

13.00 

12.50 

.  6. 486E-01 

8.108E 

00 

14.00 

13.50 

5.001E-01 

6.751E 

00 

15.00 

14.50 

3.798E-01 

5.507E 

00 

16.00 

15.50 

2.807E-01 

4.351E 

00 

17.00 

16.50 

2.012E-01 

3.320E 

00 

18.00 

17.50 

1. 380E-01 

2.415E 

00 

19.00 

18.50 

8.816E-02 

1.631E 

00 

20.00 

19.50 

4.962E-02 

9.676E- 

•01 

21.00 

20.50 

2.151E-02 

4.410E- 

•01 

22.00 

21.50 

4.117E-03 

8.851E-02 

20.0  NEV  SPECTRUM 


ENERGY 

BIN 

PHOTON 

ENERGY 

PHOTON 

FLUENCE 

ENERGY 

FLUENCE 

1.00 

.50 

7 . 672E-01 

3.836E- 

•01 

2.00 

1.50 

1.157E  01 

1.735E 

01 

3.00 

2.50 

8.700E  00 

2.175E 

01 

4.00 

3.50 

6.089E  00 

2.131E 

01 

5.00 

4.50 

4.317E  00 

1.943E 

01 

6,00 

5.50 

3.155E  00 

1.735E 

01 

7.00 

6.50 

2.357E  00 

1.532E 

01 

8.00 

7.50 

1.788E  00 

1.34IE 

01 

9.00 

8.50 

1.369E  00 

1.163E 

01 

10.00 

9.50 

1.052E  00 

9.992E 

00 

11.00 

10.50 

.  8.072E-01 

8.476E 

00 

12.00 

11.50 

6.158E-01 

7.082E 

00 

13.00 

12.50 

4.636E-01 

5.795E 

00 

14.00 

13.50 

3.412B-01 

4.606E 

00 

15.00 

14.50 

2.439E-01 

3.536E 

00 

16.00 

15.50 

1.650E-01 

2.558E 

00 

17.00 

16.50 

1.039E-01 

1.714E 

00 

18.00 

17.50 

5.750E-02 

1.006E 

00 

19.00 

18.50 

2.439E-02 

4.512E-01 

20.00 

19.50 

3.769E-03 

7.350E-02 
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18.0  MEV  SPECTRUM 


221 


18.0  MEV  SPECTRUM  (+/-  2%) 


ENERGY 

BIN 

PHOTON 

ENERGY 

PHOTON 

FLUENCE 

ENERGY 

FLUENCE 

1.00 

.50 

6.190E-01 

3.095E-01 

2.00 

.1.50 

9.406E  00 

1.41 IE  01 

3.00 

2.50 

7.030E  00 

1.758E  01 

4.00 

3.50 

4.885E  00 

1.710E  01 

5.00 

4.50 

3.434E  00 

1.545E  01 

6.00 

5.50 

2.486E  00 

1.367E  01 

7.00 

6.50 

1.836E  00 

1.193E  01 

8.00 

7.50 

1.372E  00 

1.029E  01 

9.00 

8.50 

1.031E  00 

8.767E  00 

10.00 

9.50 

7 . 737E-01 

7.350E  00 

11.00 

10.50 

5.755E-01 

6.043E  00 

12.00 

llvSO 

4.204E-01 

4.835E  00 

13.00 

12.50 

2.984E-01 

3.730E  00 

14.00 

13.50 

2.022E-01 

2.729E  00 

15.00 

14.50 

1 . 270E-01 

1.841E  00 

16.00 

15.50 

6.965E-02 

1.080E  00 

17.00 

16.50 

2.927E-02 

4.830E-01 

18.00  . 

17.50 

5.312E-03 

9.297E-02 

15.0  MEV  SPECTRUM 


ENERGY 

BIN 

PHOTON 

ENERGY 

PHOTON 

FLUENCE 

ENERGY 

FLUENCE 

1.00 

.50 

4.2I8E-01 

2.109E-01 

2.00 

1.50 

6.473E  00 

9.709E  00 

3.00 

2.50 

4.771E  00 

1.193E  01 

4.00 

3.50 

5.260E  CO 

1.1415  01 

5.00 

4.50 

2  248E  00 

1.012E  01 

6.00 

5. 50 

1.589E  00 

8.738E  00 

7.00 

6.50 

1.139E  00 

7.400E  00 

8.00 

7.50 

8 . 175E-01 

6.131E  00 

9.00 

8.50 

5.814E-01 

4.942E  00 

10.00 

9.50 

4 . 046E-01 

3.844E  00 

11.00 

10.50 

2.7C3B-CI 

2.838E  00 

12.00 

11.50 

1.681E-01 

1.934E  GO 

13.00 

12.50 

9. 253E-02 

1.1576  CO 

14.00 

13.50 

3.890E-02 

5.251E-01 

15.00 

14.50 

6.128E-03 

8 .886E-02 
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15.0  MEV  SPECTRUM  (+/-  2%) 


ENERGY 

BIN 

PHOTON 

ENERGY 

PHOTON 

PLUENCE 

ENERGY 

PLUENCE 

1.00 

.50 

4 .219E-01 

2 . 109E-01 

2.00 

1.50 

6.474E  00 

9.711E  00 

3.00 

2.50 

4.772E  00 

1.193E  01 

4.00 

3.50 

3.261E  00 

1.141E  01 

5.00 

4.50 

2.249E  00 

1.012E  01 

6.00 

5.50 

1.589E  00 

8.742E  00 

7.00 

6.50 

1.139E  00 

7.402E  00 

8.00 

7.50 

8 . 176E-01 

6.132E  00 

9.00 

8.50 

5.818E-01 

4.945E  00 

10.00 

9.50 

4.050E-01 

3.848E  00 

11.00 

1C. 50 

2.706S-01 

2.B42E  00 

12.00 

11.50 

1.686E-01 

1.938E  00 

13.00 

12.50 

9. 262E-02 

1.158B  00 

14.00 

13.50 

3 .911E-02 

5.280E-01 

15.00 

14.50 

6.831E-03 

9.905E-02 

